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Abstract 



This thesis is concerned with two topics in Hamiltonian lattice gauge theory (LGT): im- 
provement and the apphcation of analytic techniques. On the topic of improvement, we 
develop a direct method for improving lattice Hamiltonians for gluons, in which linear 
combinations of gauge invariant lattice operators are chosen to cancel the lowest order 
discretisation errors. A number of improved Hamiltonians are derived and the level of 
improvement tested for the simple cases of U(l) in 3+1 dimensions and SU(2) in 2+1 di- 
mensions. In each case the improved calculations are closer to the continuum limit, when 
working at a given lattice spacing, than their unimproved counterparts. 

On the topic of analytic methods, we extend the techniques that have been used in 
2+1 dimensional SU(2) variational calculations for many years, to the general case of 
SU(Ar). For this purpose a number of group integrals are calculated in terms of Toeplitz 
determinants. As generating functions these group integrals allow the calculation of all 
matrix elements appearing in variational glucball mass calculations in 2+1 dimensions. 
Making use of these analytic techniques, the lowest five glueball masses, or massgaps, in 
various symmetry sectors are calculated, with improved and unimproved Hamiltonians, 
for 2 < A?^ < 5. For simplicity, a minimisation basis containing only rectangular states 
is used. Large scaling regions are obtained for each of the five lowest mass states in the 
J^^ = 0"*""^ and sectors. The results are comparable to competing calculations. How- 
ever, the enumeration of the states obtained requires further work. 

The large N limit of the glueball mass spectrum is also explored. Making use of an- 
alytic variational techniques, glueball mass calculations with N as large as 25 are carried 
out on a desktop computer and the results extrapolated accurately to the N ^ oo limit. 
Evidence of leading order finite N corrections is obtained for a number of excited 



states. An interesting empirical observation is made which demonstrates a close similarity 
between the glueball mass spectrum and the simple harmonic oscillator spectrum in 2+1 
dimensions. 

We finish with a feasibility study of applying analytic Hamiltonian techniques in 3+1 
dimensional calculations. After developing analytic techniques specific to 3+1 dimensions, 
we study the variational glueball mass spectrum for SU(iV) on a single cube, with a min- 
imisation basis containing only two states. Memory limitations do not allow calculations 
with N > 8. We observe promising signs of an approach to asymptotic scaling in the 1"' 
sector as A/" is increased, warranting larger volume studies with additional states in the 
minimisation basis. 
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Chapter 1 



Introduction 

1.1 Historical Background 

Quantum chromodynamics (QCD) is the accepted theory of high energy strong inter- 
actions. It is the gauge theory describing the interactions of coloured particles, quarks 
and gluons, which are commonly accepted to be the fundamental degrees of freedom of 
hadronic matter. 

The non-abelian nature of QCD results in a collection of non-trivial complications 
when compared to the archetypal gaiige theory quantum electrodynamics (QED). Most 
serious is the phenomenon of asymptotic freedom which dictates that the QCD coupling 
constant, , is not in fact constant. Rather it depends on the energy scale of the process in 
question. The calculation of the one-loop (3 function [3, 4] demonstrated that the QCD cou- 
pling constant is small only for high energy processes; a fact observed in experiment. This 
discovery was instrumental in the establishment of QCD as a theory of strong interactions. 

Asymptotic scaling poses a serious problem to quantitative tests of QCD. The tra- 
ditional perturbative techniques used in high precision tests of QED rely on expansions 

about a small coupling constant. Consequently such techniques arc only applicable to high 
energy QCD processes. Perturbative QCD has had a great deal of success explaining such 
processes. 

The first quantitative demonstration of the validity of perturbative QCD was in the 
violation of Bjorken scaling in deep-inelastic lepton-hadron scattering (for example see 
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Ref. [5]). More recently high order QCD calculations of jet rates in hadron colliders have 
shown impressive agreement with experimental data. Less striking tests of perturbative 
QCD have been made in the analysis of r lepton decays, heavy-quarkonium decay and 
also in e+e~ collisions. For a recent review see Ref. [6]. 

For low energy processes the QCD coupling constant is large. Consequently traditional 
perturbative techniques can not be used to test the validity of QCD in this regime. A great 
deal of interesting physics appears in this nonperturbative regime. One would like to test 
whether the same theory that successfully explains a broad range of high energy processes 
also gives rise to quark confinement, provides a mechanism for chiral symmetry breaking, 
gives rise to the correct values for the hadronic mass spectrum and explains other low- 
energy phenomena. Without the use of non-perturbative techniques the situation is dire. 
One cannot calculate such mental quantities as the masses of mesons and baryons even 
if given the coupling constant and the constituent quark masses. The need to probe the 
low-energy regime of QCD has given rise to many QCD motivated models [7, 8, 9, 10]. 
However, only one first principles approach to nonperturbative QCD is currently known 
and that is lattice gauge theory (LOT). 

LGT was originally developed in the Lagrangian formulation by Wilson in 1975 [11]. 
A year later the corresponding Hamiltonian formulation was published by Kogut and 
Susskind [12]. The basic approximation of LGT is the replacement of continuous and 
infinite space-time by a discrete lattice of sites and links. The essential difference between 
the Hamiltonian and Lagrangian formulations is that in the former time is left continuous. 
Using LGT one can presumably calculate the full range of hadronic properties from mass 
spectra and wave functions to decay constants and matrix elements. The hadronic matrix 
elements appearing in the weak decays of hadrons are of particular interest. Phenomenol- 
ogists are waiting for lattice estimates to such matrix elements in order to improve the 
theoretical understanding of the CKM matrix [13, 14]. Another area of recent interest has 
been in the study of QCD in extreme conditions. A knowledge of the phase diagram of 
QCD is relevant to the understanding of high density quark systems such as the early uni- 
verse, neutron stars and heavy ion collisions. The scope of LGT is summarised in Ref. [15]. 

In practice LGT cannot, at this stage, predict the complete range of hadronic phenom- 
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ena. The major stumbling block is the requirement of time on the largest supercomputers 
in existence to calculate basic quantities with accuracy. In the early days of LGT strong 
coupling expansions in the Hamiltonian approach provided compelling evidence for the 
confinement of quarks by showing that the potential energy between two infinitely heavy 
quarks grows linearly with their separation [12, 16]. This was without the complication of 
dynamical quarks and in the strong coupling limit which is far from the real world physics 
of the continuum. The inclusion of two flavours of quarks and the extension to physical 
couplings came much later and required the numerical techniques of Lagrangian LGT. The 
qualitative result of the strong coupling expansions did not however change (for example 
see Ref. [17] and references within). 

Most LGT calculations to date have been performed using Monte Carlo techniques in 
the Lagrangian approach. Until recently most simulations used the so called quenched 
approximation in which the vacuum polarization effects of the QCD vacuum are explicitly 
removed. The calculation of certain hadron masses are within 10% of the experimentally 
observed results using the quenched approximation [18]. Calculations utilising full QCD 
with two flavours of quarks have been performed and show significant improvements over 
the quenched approximation [19]. At this stage the advancement of LGT in the Lagrangian 
approach is restricted by the available computational resources. 

There are some areas of research where progress in the Lagrangian formulation has 
been significantly slow. Some examples are in calculations involving excited states and 
the determination of the QCD phase diagram. Such slow progress and the need for sig- 
nificant computational resources in Lagrangian LGT suggests that alternative methods 
should be pursued in parallel. A viable alternative is Hamiltonian LGT, which may have 
a number of advantages. It has been suggested that Hamiltonian LGT could more readily 
handle finite density QCD [20]. Indeed, the problems encountered in finite density QCD 
in the Lagrangian formulation has prompted a return to the strong coupling expansions 
of early Hamiltonian LGT [21]. Additional areas in which Hamiltonian LGT may have 
an advantage over the Lagrangian approach are in calculations involving excited sites, 
time dependent quantities and wave functions. An appealing aspect of the Hamiltonian 
approach is that it reduces LGT to a many body problem. As such a host of analytic 
techniques from many body theory are applicable. Techniques such as the coupled cluster 



4 



Introduction 



method [22, 23, 24, 25], plaquette expansion [26, 27, 28, 29] and variational methods [1, 30] 
have had success in recent years in calculations of glucball masses. These calculations how- 
ever have been performed without the inclusion of quarks, often in less dimensions that 
the physical 3+1 and often using lower dimensional gauge groups than the physical case 
of SU(3). 

The self interactions of gluons in QCD allow the existence of massive states containing 
only gluon degrees of freedom [31, 32]. Although QCD provides a theory of strong inter- 
actions we have little understanding of the physical states of the theory. The constituent 
quark model provides an excellent phenomenological description of meson and baryon 
spectroscopy. QCD inspired models however predict the existence of other hadrons. These 
hadrons include not only the glueballs, which are often described as bound states of two 
and three gluons, but other more exotic states as well. So called hybrid particles, which 
have both quark and excited gluon degrees of freedom, have been predicted to exist. Until 
the properties of the physical states of QCD can be reliably calculated, and these calcula- 
tions stand up to experimental tests, we cannot claim to understand the nonpcrturbative 
aspects of QCD. Prom both theoretical and experimental sides there have been significant 
developments in the search for glueballs in recent years. 

On the theoretical side calculations of the glueball mass spectrum have been dominated 
by LGT calculations in the Lagrangian approach. The calculations of the last decade are 
nearing agreement at least for the lowest mass glueballs. The emerging picture is that the 
scalar glueball with quantum numbers J^*^ = 0+"'" is the lightest with a mass in the range 
1.4-1.8 GeV. Interestingly the next lowest mass states in the PC = ++ sector have spin 
ordering 2, 3, 1, whereas the ordering in the PC = — h sector is 1, 3, 2, 0. The status of 
LGT glueball calculations in the Lagrangian approach is reviewed in Refs. [33] and [34]. 

On the experimental side progress has been hampered by numerous technical difficul- 
ties. Perhaps most significant is the fact that glueballs can mix strongly with nearby qq 
resonances. To search for glueballs one needs to probe glue rich processes. Recent searches 
have taken place in radiative J/ ip decay, pp annihilations and central production in pp col- 
lisions. A review of search techniques is presented in Ref. [35]. Perhaps the most promising 
glueball candidates are the /o(1500) and /o(1710) resonances (the subscript denotes the 
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total angular momentum and the argument the mass in MeV) . The interpretation of these 
resonances has been the subject of much debate, particularly with regard to mixing with qq 
states. Some authors argue that the /o(1500) has a dominant glueball contribution while 
the glueball content of /o(1710) is minimal [36, 37]. Other authors argue the opposite 
case [38, 39]. Improved data from non-appearance searches in 77 collisions could possibly 
resolve this issue [40]. The basis for these studies lies in the fact that the glueball to two 
photon widths are small compared to those of mesons since gluons do not couple directly 
with photons. Such studies will also be important in the identification of other glueball 
candidates such as /j(2220) [41]. 

The focus of most of the lattice community, as far as glueballs are concerned, has been 
the physically interesting case of SU(3) glueballs in 3+1 dimensions. In recent years the 
study of glueballs belonging to higher dimensional gauge groups in both 2+1 and 3+1 
dimensions has received intense interest from the string theory community. The reason 

for this is that the SU(A^) glueball mass spectrum in the large N limit provides an ideal 
laboratory for quantitative tests of Maldacena's conjecture and its later extensions. We 
defer the discussion of this topic to Section 7.2. 

While accurate studies of the SU(3) glueball spectrum in 3+1 dimensions may be some 
way off in Hamiltonian LGT, studies of the SU(iV) glueball spectrum in 2+1 dimensions 
are definitely feasible. In fact it appears that the Hamiltonian approach has a distinct 
advantage over the Lagrangian approach in that excited states are easily accessed. Cur- 
rent Lagrangian calculations have extracted the three lowest mass states in some sectors 
in the large N limit. With Hamiltonian techniques it is feasible to even triple this number. 
The accurate simulation of large mass excited states is of immediate interest to QCD in- 
spired models of glueballs such as the bag model [7] , potential models [8] and the flux-tube 
model [9] and its extensions [10]. The last of these in fact shows remarkable agreement 
with the lowest excited state masses calculated on the lattice in 2+1 dimensions for all 
but the 0"+ state [10]. 

The results so far obtained in Lagrangian LGT have been made possible by an ex- 
tensive improvement programme initiated by Symanzik nearly twenty years ago [42, 43]. 
The formulation of QCD on a lattice introduces errors which can be expressed in powers 
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of the lattice spaeing a and the QCD couphng and start at ©(a^). When performing a 
lattice calculation in the Lagrangian approach the most important factor in determining 
the computational cost is the lattice spacing. In fact it has been predicted that the cost of 
a full QCD simulation varies as a"^-^ [44]. Thus it is very important to be able to work on 
coarse lattices. This of course is a move away from the a ^ limit in which the real world 
physics of the continuum is retrieved. This realisation has motivated the development of 
improved actions (for a review see Ref. [45]). These actions aim to systematically remove 
the errors caused by placing QCD on a lattice. In practice these actions allow one be closer 
to the continuum limit when working with a given lattice spacing. The additional compu- 
tational cost of working with the more complicated improved actions is by far outweighed 
by the benefits of working on coarse lattices. The development of improved actions has 
been a crucial factor in the progress of LGT in the last 15 years. It has allowed routine 
calculations to be performed on work stations and has brought the numerical solution of 
full QCD tantalisingly close. 

In contrast the improvement of Hamiltonian LGT has hardly begun. Perhaps the most 
comprehensive study into the derivation of improved Hamiltonians to date is due to Luo, 
Guo, Kroger and Schiitte [46]. Their study started with an improved Lagrangian and 
derived improved Hamiltonians with the standard transfer matrix and Legendre transfor- 
mation techniques. Using these techniques it became clear that the construction of an 
improved kinetic Hamiltonian with a finite number of terms is a non-trivial problem. Us- 
ing a more direct method improved Hamiltonians can be constructed with relative ease. In 
this thesis we propose a method for the direct improvement of Hamiltonian LGT which is 
in the spirit of the original Kogut-Susskind derivation. This topic is the focus of Chapter 3. 

We now move on to a description of the outline of the thesis. 

1.2 Outline of the Thesis 

As a form of introduction, in Chapter 2, the fundamental approximations of LGT are 
detailed. The basic steps involved in a Hamiltonian lattice calculation are also discussed 
as well as the important topic of how the real world physics of the continuum is retrieved. 
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Chapters 3 and 4 arc concerned primarily with the topic of improvement. In Chapter 3 
we develop a technique which allows the straightforward construction of lattice Hamiltoni- 
ans for gluons in d+1 dimensions. A variety of improved Hamiltonians are derived using 
this method. These improved Hamiltonians are devoid of the lowest order discretisation 
errors present in the unimproved Kogut-Susskind Hamiltonian. In Chapter 4 we present 
some simple checks on improvement for the cases of U(l) in 3+1 dimensions and SU(2) in 
2+1 dimensions. In this way we test the level of improvement achieved by a selection of 
the Hamiltonians derived in Chapter 3. 

In Chapter 5 we move on to the topic of analytic methods. Here we introduce tech- 
niques which allow analytic lattice calculations of SU(iV) glueball masses in 2+1 dimen- 
sions for general N. Such techniques extend the SU(2) methods of Arisue [1] to the general 
case of SU(Ar). These techniques are applied in Chapter 6 in the calculation of glueball 
masses in 2+1 dimensions with 2 < N < 5. Kogut-Susskind, improved and tadpole im- 
proved Hamiltonians, all of which are derived in Chapter 3, are used and their results 
compared. 

In Chapter 7 we move on to an analytic study of large A'' physics on the lattice. Here 
we study the large N glueball mass spectrum in 2+1 dimensions using the Kogut-Susskind 
Hamiltonian. Using the analytic techniques of Chapter 5, calculations with N as large as 
25 are possible on a desktop computer, allowing an accurate extrapolation to the N ^ oo 
limit. Chapters 6 and 7 contain the key results of this thesis. 

In Chapter 8 we move on to explore the viability of adapting the analytic techniques, 
Tiscd with success in Chapters 6 and 7, to the physically interesting case of 3+1 dimen- 
sions. Here, the difficulties of a move to 3+1 dimensions are discussed. Calculations of 
glueball masses within a simplistic, one cube, model are presented as a first step in testing 
the worth of analytic variational techniques in 3+1 dimensions. 

Chapter 9 contains a summary of our findings and our conclusions. 



Introduction 



Chapter 2 



Gauge Theory on a Lattice 



2.1 Introduction 

In this chapter wc describe the basic steps involved in Hamiltonian LGT calculations. 
Our intention is to introduce the fundamental approximations of LGT and the topics 
of importance to recovering the real world physics of the continuum. We start with an 
explanation of how quarks and gluons are put on a lattice in Section 2.2. In this section 
the important difference in the building blocks of LGT and continuum gauge theory are 
discussed. We then move on to explain the steps involved in performing a calculation in 
Hamiltonian LGT. We finish in Section 2.6 with a discussion of the important topic of 
how the real world physics of the continuum is extracted from a calculation on the lattice. 



2.2 The Lattice 

The fundamental simplification of LGT is the replacement of continuous and infinite space- 
time with a finite and discrete lattice of sites and links. A label, I = {x,i), is attached 
to each link. Here x is the starting point of link I and i its direction. Quark fields are 
essentially treated in the same way on the lattice as they are in continuum QCD. In both 
theories they are represented by Grassmann variables. However, in LGT the quark fields 
are restricted to lattice sites. As in the continuum, the fermion fields can be integrated 
out of the partition function exactly giving a highly non-local fermion determinant, det M. 
The calculation of this determinant is the key factor in the computational expense of full 
QCD simulations in Lagrangian LGT. The quenched approximation, in which the vacuum 
polarisation effects of the QCD vacuum are explicitly removed, gives detM a constant 
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Figure 2.1: The link operator. 



value, reducing both the computational expense and the connection to real world physics 
of LGT. It is in this approximation that the bulk of Lagrangian LGT calculations have 
been carried out. 



In LGT gluons are incorporated differently to continuum gauge theory in order to 
maintain manifest gauge invariance on the lattice. In continuum SU(A^) gauge theory 
gluons arc represented by the SU(iV) valued fields A^{x). In LGT gluons are represented 
by Wilson lines (or link operators) along the links joining lattice sites. We denote the link 
operator on link I hy Ui. The link operator on the same link but taken in the opposite 
direction is given by Ui . At times it will be convenient to use the more detailed notation, 
Ui{x), for the link operator on link I = {x,i). It is conventional to represent the link 
operator by a directed line segment joining the sites x and x + ai, as shown in Fig. 2.1. 
Here i is a unit vector in the i direction. The details of this construction and precisely how 
it relates to the continuum theory will be introduced in Section 3.3.2. The fundamental 
differences between continuum gauge theory and LGT are caricatured in Fig. 2.2. As in 
Fig. 2.2 it is conventional to denote the lattice spacing by a. 



Having discussed the fundamental approximation of LGT and briefly introduced its 
building blocks we now move on to a description of the steps involved in a complete 
Hamiltonian LGT calculation. The first step is choosing an appropriate lattice. 



2.3 Defining a Lattice 
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(a) Continuum QCD 
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Figure 2.2: A cartoon comparing continuum QCD with lattice QCD. 



2.3 Defining a Lattice 



In Hamiltonian LGT one works on a spatial lattice. The time coordinate remains contin- 
uous. Contrastingly, the Lagrangian formulation discretises all dimensions and works in 
Euclidean space. An obvious advantage of the Hamiltonian approach is that one works 
with a lattice with one less dimension. The cartoon in Fig. 2.2 illustrates a cubic lattice 
in three dimensions. LGT is not, however, restricted to 3 dimensions and the lattice itself 
is not restricted to a cubic arrangement. 



For reasons of simplicity, it has been conventional to study LGTs built on cubic lattices. 
Not a great deal of effort has been devoted to exploring the effect of using lattices with 
different structures [47]. One could however imagine that working on a non-cubic lattice 
could be advantageous in some situations. For example, if the degrees of freedom were 
arranged into products of link operators (which is often convenient) around the elementary 
faces of the lattice (for example squares of side length a in the cubic case), working on 
non-cubic lattices could be advantageous when the number of faces and number of links 
need to be related in a precise way to incorporate constraints such as Gauss' law [24]. In 
this thesis we work on a (i dimensional spatial lattice with d = 2 and 3. 



When performing numerical LGT calculations the primary restriction is computer 
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memory. The value of every link operator on the lattice needs to be stored in mem- 
ory at a given time. Since the number of links on the lattice increases quickly as the 
volume increases, one is in practice limited to lattices with small volumes. The majority 
of work in LGT to date has involved small volume numerical simulations in the Lagrangian 
approach. Typical calculations on supercomputers with speeds of the order of hundreds 
of gigaflops^ with two flavours of dynamical quarks are performed on 24^ x 48 lattices. It 
has been estimated that with access to a 10 teraflop machine 36^ x 96 lattices will be able 
to be used for equivalent calculations with similar effort [44]. 

Variational studies in Hamiltonian LGT in 3+1 dimensions (which go beyond the 
strong coupling limit), to our knowledge, have only dealt with small lattices; typically 
with sizes smaller than 6^ [48, 49]. We discuss the viability of using analytic Hamiltonian 
techniques in 3+1 dimensions in Chapter 8. 

Having defined a lattice we now move on to the important step of constructing a 
Hamiltonian. 

2.4 Constructing a Hamiltonian 

When constructing a Hamiltonian on the lattice there are three important points to con- 
sider. Firstly, the lattice Hamiltonian must reduce to its continuum counterpart in the 
limit of vanishing lattice spacing. Prom here on we refer to this limit as the continuum 
limit. Secondly, we must incorporate as many symmetries of the continuum theory as pos- 
sible. Finally, the lattice Hamiltonian must not be prohibitively expensive to use. These 
points are closely related and often an attempt to address one point is detrimental to 
another. For example, if one attempts to build a Hamiltonian that is closer to its con- 
tinuum counterpart for a given lattice spacing, one with a weak dependence on a, one is 
left unavoidably with a more computationally intensive Hamiltonian. Another example 
relates to restoring rotational symmetry. On the lattice, continuous rotational symmetry 
is broken down to a discrete analogue which depends on the particular lattice in use. One 
can attempt to construct Hamiltonians with improved rotational symmetry but again the 
resulting Hamiltonian is more computationally expensive than the simplest versions. 

^ 1 flop = 1 floating point operation per second 
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One symmetry that is essential to incorporate in the lattice Hamiltonian is local gauge 
invariance. One may ask why it is important to uphold local gauge invariance when other 
symmetries such as rotational invariance and Lorentz invariance are explicitly broken. 
The reason lies in the fact that with local gauge invariance upheld the bare gluon mass 
is zero and the couplings appearing at multiple gluon and gluon-quark vertices are all 
equal. Without local gauge invariance these couplings would need to be tuned indepen- 
dently in order to retrieve real world physics. This is a complication that would become 
prohibitively expensive when processes with loop diagrams were taking into consideration. 

In building lattice Hamiltonians we have a great deal of freedom. Our only restrictions 
are gauge invariance and the need to obtain the correct continuum limit. Our approach 
is to start with gauge invariant lattice operators and form linear combinations such that 
the correct continuum limit is obtained. In this way Hamiltonians, which couple links 

with various separations, can be easily constructed. This process is discussed in detail in 
Sections 3.3 and 3.4. 

Having constructed a Hamiltonian the next step in a Hamiltonian LGT calculation is 
to choose an appropriate wave function. 

2.5 Choosing a Wave Function 

In this section we discuss the issues involved in choosing a wave function. We start with 
a discussion of lattice symmetries and finish by mentioning some practical concerns. 

When performing a LGT calculation one is usually interested in studying a particle 
with certain symmetry properties. For instance in this thesis we study the masses of 
glueballs in 2+1 and 3+1 dimensions. These glueballs can be classified according to their 
total angular momentum, J, charge conjugation, C, and parity, P eigenvalues. In order to 
study these particles we need to construct suitable wave functions. The restrictions on the 
wave function are firstly, that it be gauge invariant, and secondly, possess the symmetry 
properties of the continuum state we wish to study. As we will see in Section 3.3.3, the 
former forces us to construct our state from Wilson loops, the latter requires us to choose 
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particular combinations of Wilson loops with certain symmetry properties. 

Suppose for example we work with a Hamiltonian symmetric under parity transforma- 
tions^, P, charge conjugation, (7, and rotations by 7r/2. Then the Hamiltonian commutes 
with C, P and the lattice spin operator J. Consequently the eigenstates of the Hamil- 
tonian can be labelled by their respective (7, P, and J eigenvalues. C has the effect of 
replacing the lattice coupling, which we denote by e, by — e. This is equivalent to the 
replacement Ui uj. To define the parity operator and explore rotations on the lattice 
we need to be more careful. We start with a description of parity. 

With three spatial dimensions, the effect of P is to replace lattice sites x by —x. 
However in two spatial dimensions, with one less axis to rotate about, the same prescription 
(of course with the z component taken to be zero) is equivalent to a pure rotation. In 2-1-1 
dimensions we can define P as a reflection in one axis. For example, in this thesis we will 
take 

P|V'(x,y)) = |V^(x,-y)). (2.1) 
Next we discuss the more difficult topic of rotational symmetry. 

The rotational symmetry properties of lattice operators are restricted to the particular 
rotational symmetries of the lattice. In the 3-1-1 dimensional continuum rotational symme- 
try is continuous and the spin of a state is characterised by the irreducible representations 
of the rotation group SO (3). The situation on the lattice is quite different. On a cubic 
lattice the 80(3) rotational symmetry of the continuum is broken down to the cubic or 
octahedral group O. Since the lattice Hamiltonian is symmetric under cubic rotations its 
eigenstates must be classified in terms of the irreducible representations of O, not of S0(3). 
The classification of lattice states in terms of the irreducible representations of O was de- 
veloped by Berg and Billoire [50]. An analogous situation occurs in 2-1-1 dimensions. As 
we will see in Section 6.5.2 the result is an ambiguity in spin identification on the lattice. 
For instance on a 2 dimensional square lattice one can only determine the continuum spin 
of a state modulo 4 with standard techniques [51]. One can determine the spins of the 



^We use the notation, X, to represent an operator and, X to represent an eigenvalue. 
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lowest mass states when using smeared links within the Lagrangian approach [51]. This 
approach, however, does not allow the spins of excited states to be determined. Techniques 
are available to counter this problem [52, 53] but currently their use is at an exploratory 
stage. 

In 3+1 dimensions the situation is less problematic. Since there is an additional axis to 
rotate about, one can count the multiplicity of states to help determine continuum spins. 
Additional methods, specific to 3+1 dimensions, have been developed to distinguish spin 
and spin 4 states with a careful choice of lattice operators [54]. 

A further consideration in the construction of appropriate wave functions is transla- 
tional invariance. On the lattice, states with particular momenta are constructed in the 
usual way via Fourier transforms. In this thesis we will be concerned with the calculation 
of masses, for which translationally invariant zero momenta states suffice. 

The last issue to discuss is one of a practical nature. A benefit of the Hamiltonian 
approach is that analytic techniques are available for use. This is only the case if one 
makes use of wave functions which result in expectation values for which analytic results 
are available. If the expectation values are not able to be handled analytically one must 
resort to Monte Carlo techniques for their calculation [49] . Unsurprisingly wave functions 
which allow analytic calculations have less appealing continuum limits than those which 
must be handled numerically [49]. 

Having discussed the issues involved in choosing a trial state we now move on to discuss 
the important topic of how continuum physics is extracted from a lattice calculation. 

2.6 Extracting Continuum Physics 

In this section we explain how continuum results are extracted from a lattice calculation 
in 2+1 and 3+1 dimensions. This will serve as an introduction to the form in which the 
results of later chapters are presented. In the context of this thesis the primary difference 
between 2+1 and 3+1 dimensions is the dimension of the coupling constant. As we shall 
see in Section 3.3.5 the coupling constant, e^, is dimensionless in 3+1 dimensions but has 
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the units of mass in 2+1 dimensions. We start with a general discussion of how continuum 
results are extracted from a LGT calculation. We then specialise to the cases of 3+1 and 
2+1 dimensions. 

When performing a calculation on the lattice one is always interested in eventually 
extracting continuum physics. The process of unmasking continuum physics however is 
highly nontrivial. Naively one would expect continuum results to be obtained simply in 
the limit as a — > 0. However, if this naive limit is taken in a lattice calculation of a 
given physical quantity either zero or infinity is obtained depending on the dimension of 
the quantity under consideration. To achieve physical results one must tune the coupling 
constant in such a way that a finite result is obtained for all physical quantities. 

To illustrate the process consider the calculation of a physical quantity, Q, on the 
lattice with dimension in units of the lattice spacing d{Q) . As an example Q could be the 
mass of a physical state, in which case d{Q) = — 1. In Hamiltonian LGT one typically 
calculates Q for a number of couplings, arriving at 

Q = a'^(^^f{e). (2.2) 

Here / is a dimensionless function which may depend on the coupling, e, either directly 
or via the parameters of the wave functions used in the calculation of Q. It is clear from 
Eq. (2.2) that in the naive limit, a ^ , Q approaches cither infinity or zero (depending 
on the sign of d{Q)). It is clear that for a finite result one must let e vary with a in such a 
way that a sensible result is obtained for Q as a — > 0. It is only in this limit that the real 
world physics of the continuum is recovered. The nontrivial aspect of this limit is imme- 
diately clear. The dependence of e on a must be chosen so that all physical observables 
have a well defined a ^ limit. The complexity lies in the fact that for each physical 
observable Q there is a different /(e) in Eq. (2.2). The coupling must be tuned in such 
a way that as it approaches its a — limit all physical observables tend to constant values. 

The lattice is a regularisation scheme which reduces the continuum theory to a theory 
with a finite number of degrees of freedom. The process of tuning the coupling constant 
described above is the renormalisation process which precisely removes the dependence of 
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the calculation on the regularisation scheme. Precisely how the coupling depends on a is 
dictated by the Callan-Symanzik /3 function 

m = a% (2.3) 

One can make specific predictions about the form of /3 in weakly coupled QCD. Including 
up to two loop contributions it has been shown [55, 56] that for 3+1 dimensional SU(A/") 
gauge theory with n/ flavours of quarks, 
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^ 3-"^- iV +C?(e')-(2.4) 



A four loop calculation of /3(e) is also available [57]. In this thesis we will always be 
concerned with the quarkless case n/ = 0. Making use of Eq. (2.4) one can integrate the 
one loop part of Eq. (2.3) directly to obtain for the quarkless case 

where ^ is the inverse 't Hooft coupling ^ = {Ne^)~^. The integration constant, A^,, with 
dimensions of mass is commonly called the lattice scale parameter. It can be related per- 
turbatively, at the one [58, 59] and two loop level [60], to the scale parameters obtained 
from other regularisation schemes. Such results are important when comparing results 
from LGT to results obtained with other regularisation schemes. It is important to note 
that by inverting Eq. (2.5) it is clear that we must have e — *■ as a — *■ 0. Since they 
are different regularisation schemes, the scale parameters of Hamiltonian and Lagrangian 
LGT also differ. Results for their ratios are available [59, 61] and are necessary for the 
comparison of results obtained in the two formulations in the continuum limit. 



When attempting to extract continuum physics from a Hamiltonian LGT calculation 
one generally looks for a range of couplings for which Eq. (2.5) applies. If such a range 
of couplings, called a scaling window, exists all observables with dimensions of mass must 
scale in the same way as A^; all masses must be proportional to A^ in a scaling window. 
The window does not, in practice, extend all the way to e = 0. Generally the approxima- 
tions used to perform the calculation break down before then. One generally assumes that 
if the approximations used in the calculations were improved the scaling window would 
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extend closer to e = 0. 

To demonstrate how continuum physics is extracted in practice consider the calculation 
of a mass on the lattice in pure SU(iV) gauge theory. We calculate the mass, M(e) for a 
range of couplings. To observe asymptotic scaling we require 

aM(e) = M*akL (2.6) 

for some range of couplings. Here M* is the continuum limit mass in units of A^. Making 
use of Eq. (2.5) we see that 

\og{aM{e)r'"''') = log(M*) + ^ log - (2-7) 

Thus in a scaling window a plot of log(aM(e)^^^"'^/"'^^"'^) versus ^ yields a straight line with 
gradient — 247r^/ll. The observation of such behaviour allows the constant of propor- 
tionality, M*, to be extracted. We refer to the type of scaling where all masses scale 
as Eq. (2.5) as "asymptotic scaling". If such behaviour is observed one can be confident 
that the lattice spacing is small enough that continuum physics is being revealed and an 
accurate value for the continuum limit mass is obtained. 

A less general type of scaling can be seen in the ratios of masses. Since all masses 
scale as Eq. (2.6) for small enough lattice spacing, the ratios of masses become constant 
in a scaling window. However, a ratio of masses may happen to be constant on a much 
larger scaling window, one where the lattice scale parameter does not obey Eq. (2.5) for 
its entire range of couplings. In such a region the non-zero errors, which depend on a and 
e, cancel to some extent and the ratio reproduces the correct continuum ratio. We refer 
to this more general situation as "scaling" which is not to be confused with asymptotic 
scaling. 

The final step in obtaining a result in familiar units is to fix the value of A^^ in physical 
units (MeV for example). This is most commonly done by by comparing the string ten- 
sion -y/a calculated on the lattice in units of A^ with the physical value in units of MeV 
obtained from the spectroscopy of heavy quarkonia. The string tension is the coefficient of 
the linearly rising potential energy as a quark-antiquark pair are separated. We introduce 
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the string tension in Section 4.2. A commonly used value for SU(3) in 3+1 dimensions is 
^0^ = 440 ±38 MeV [62]. 

The above discussion treats the case of 3+1 dimensions. To renormalise the results of 
an LGT calculation on a lattice with different numbers of dimensions requires a different 
treatment. For the case of 2+1 dimensions, with which we are concerned for the bulk of 
this thesis, the situation is particularly simple. For this case the coupling constant, e^, has 
dimensions of mass, in contrast to the 3+1 dimensional case for which it is dimensionless. 
The coupling constant can therefore be used to explicitly set the scale of the lattice theory. 
This has become standard practice in 2+1 dimensional LGT, both in the Lagrangian and 
Hamiltonian approaches. Consequently, in the limit of vanishing lattice spacing all masses 
calculated on the lattice must approach a constant value when measured in units of the 
coupling, 

M(e) = M*e^ (2.8) 

Here M* is the continuum limit mass in units of e^. In 2+1 dimensional LGT, the Hamilto- 
nian and Lagrangian formulations again present different regularisations of the continuum 
gauge theory. To compare the results of the two formulations one needs to acknowledge 
that the coupling constants of the two formulations differ. Their ratio can be calculated 
perturbatively [61] for small couplings in pure SU(A'') LGT and is close to unity. 

2.7 Conclusion 

In this chapter we have introduced the basic approximations made in LGT. We have 
discussed briefly the steps involved in a Hamiltonian LGT calculation, introducing in par- 
ticular the important steps of choosing wave functions and extracting continuum physics. 
Our intention was to provide a method-independent introduction to the process of a Hamil- 
tonian LGT calculation in order to introduce the notation in which the results of later 
chapters are presented. In the next chapter we consider in detail the process of constructing 
and improving lattice Hamiltonians. 
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Chapter 3 



Constructing and Improving 
Lattice Hamiltonians 



3.1 Introduction 



In this chapter we develop a technique which allows the straightforward construction of 
Hamiltonians for gluons in d + 1 dimensions. We start with a historical background in 
Section 3.2 to motivate the need for additional Hamiltonians beyond the original Kogut- 
Susskind Hamiltonian. In Section 3.3 wc derive the simplest Hamiltonian for gluons on 
the lattice, the Kogut-Susskind Hamiltonian, focussing on the errors occurring when it 
is compared to its continuum counterpart. These errors can be divided into two classes, 
classical and quantum. The origin of each class of error is discussed in Section 3.3.2. 
We then turn our attention to the removal of classical errors in Section 3.4, deriving the 
simplest classically improved lattice Hamiltonian by adding appropriately weighted gauge 
invariant lattice operators. We refer to this technique as direct improvement. In Section 3.5 
we discuss the topic of tadpole improvement, a technique for correcting a particular type 
of quantum error, which has revolutionised Lagrangian LGT in the last 15 years. In 
particular we discuss its implementation in the Hamiltonian formulation. We finish the 
chapter in Section 3.6 with a derivation of an improved lattice Hamiltonian which couples 
distant lattice sites. This highlights the ease with which improved Hamiltonians can be 
constructed with the direct method of improvement. 
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3.2 Historical Background 

The Hamiltonian formulation of LGT was developed by Kogut and Susskind [12]. Their 
Hamiltonian was derived by demanding the correct continuum limit be obtained in the 
limit of vanishing lattice spacing. This was the style in which the Wilson action was de- 
rived a year earlier. Creutz showed that the Kogut-Susskind Hamiltonian could be derived 
from the Wilson action using the transfer matrix method [63] . Later Kogut demonstrated 
that the same could be done by taking the continuous time limit of the Wilson action and 
performing a canonical Legendre transformation [64]. 

When building a gauge theory on a lattice, the infinite number of degrees of freedom 
of continuum theory is reduced to a large but manageable number. The price to pay, how- 
ever, is that errors are made. In order to obtain an accurate extrapolation to continuum 
physics these errors need careful attention. 

As discussed in Section 1.1 much work in the last fifteen years has been devoted to 
improving lattice actions. The motivation for improvement in the action formulation is 
computational cost. Since the cost of a Lagrangian LGT calculation varies as a~^'^^ [44], 
it is by far more efficient to build improved theories than to work on finer lattices. The 
improvement programme, in particular tadpole improvement, has allowed accurate calcu- 
lations to be performed on relatively coarse lattices on desktop computers and brought 
the simulations full QCD within the reach of today's most powerful computers. 

In contrast, the improvement of lattice Hamiltonians has only recently begun. With 
most computational techniques in the Hamiltonian formulation, as one moves closer to 
the continuum limit calculations require the symbolic manipulation of an increasing num- 
ber of increasingly complicated lattice operators. The motivation for improvement in the 
Hamiltonian approach is that when working at a given lattice spacing (or coupling) one is 
closer to the continuum limit when using an improved Hamiltonian. In this way the effort 
required to reach the continuum limit is reduced, provided the improved Hamiltonian is 
not prohibitively difficult to work with. Perhaps the most extensive treatment to date is 
due to Luo, Guo, Kroger and Schiitte who discussed the classical and tadpole improve- 
ment of Hamiltonian LGT for gluons [46] . In their study it was discovered that deriving 
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an improved Hamiltonian from a classically improved action, whether by transfer matrix 
or Legendre transformation methods, results in a highly non-local kinetic term with an 
infinite number of terms. To derive a kinetic Hamiltonian with only nearest neighbour 
terms it was found necessary to start with an improved action with an infinite number of 
terms, coupling arbitrarily distant lattice sites. 

With the technique of Luo et al the order errors are removed from the Kogut- 
Susskind Hamiltonian. However, generating Hamiltonians with further improvement would 
seem exceedingly difficult. This is because one would need to start from a Liischer-Weisz 
improved action with non-planar terms [65]. For this reason we propose a move to the 
direct approach. That is, to construct improved Hamiltonians directly by adding appro- 
priate gauge invariant terms and fixing their coefficients so that errors are cancelled. 

Tadpole improvement has been studied briefly in the Hamiltonian approach [46, 66]. 
However there has been some disagreement between possible implementations. Other 
forms of quantum improvement have not yet been considered in detail in Hamiltonian 
LGT. 

In this chapter we derive the Kogut-Susskind Hamiltonian for gluons and develop a 
technique for constructing improved Hamiltonians directly. This approach is in the spirit 

of the original Kogut-Susskind derivation and parallels the strategies used for constructing 
improved actions. For simplicity we will not consider quarks at this stage. 

3.3 The Kogut-Susskind Hamiltonian 
3.3.1 Preliminciries 

In this section we derive the Kogut-Susskind Hamiltonian for pure SU(A^) gauge theory in 
d-l-l dimensions. We maintain a continuous time coordinate and work on a (i dimensional 
spatial lattice. 

We start with the formal continuum Hamiltonian operator expressed in terms of the 
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chromoelectric and chromomagnetic fields: 

H = J d'^xTi {E^ + B^) . (3.1) 

The space-like components of the continuum chromomagnetic field are given in terms of 
the space-like components of the field strength tensor, F^, by 

= 2^ijkFjk, (3.2) 

where i = 1,2, . . . ,d is the Dirac index labelling the space-like directions and summation 
over repeated indices is understood. £n...i„ is the totally antisymmetric Levi-Civita tensor 
defined to be 1 if {ii, . . . ,in} is an even permutation of {1,2, . . . ,n}, —1 if it is an odd 
permutation and otherwise (i.e. if an index is repeated). The field strength tensor is 
given in terms of the continuum gluon field, A^, by 

F/^u = d^A^ - d^A^ - ie [A,,, A^] , (3.3) 

where e is the QCD coupling. 



The chromoelectric and chromomagnetic fields can each be expressed in terms of their 
colour components as follows, 

Ei = X^Ef Bi = X'^Bf. (3.4) 

Here {A" : 1 < a < A^^ — 1} is a basis for SU(A^). It is common to use the Gell-Mann 
basis, in which case A" = G"/2, where {G" : 1 < a < A?"^ — 1} is the set of traceless N x N 
Gell-Mann matrices. It is conventional to normalise the basis as follows, 

1v(a«A'') = ^6ab, (3.5) 

so that Eq. (3.4) can be inverted to write the colour components of the chromoelectric and 
chromomagnetic fields as. 



Ef = 2Tr (SjA") and Bf = 2Tr (Sj A") . 



(3.6) 
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A final relation that is of use in this chapter is the commutation relation, 



abc \c 



(3.7) 



where the /"^^ are the totally antisymmetric, real structure constants for SU(A'^). 



3.3.2 Errors in Lattice Gauge Theory 

Before discussing Hamiltonian improvement we must first understand how deviations be- 
tween LGT and its continuum counterpart arise. The deviations can be separated into 
two classes, classical and quantum errors, which will be described in what follows. We 
start with classical errors. 



Rather than being constructed from gluon fields, pure gauge theory on the lattice is 
built from link operators 

Here e is the dimensionful QCD coupling, I labels the link in question and Ai is the lattice 
gluon field on the link I. In what follows we will adopt the convention of writing the lattice 
version of a continuum quantity in a calligraphic typeface. We define the lattice gluon field 
on link I to be the average continuum gluon field along link I, 

Ai = - fdx A. (3.9) 
« J I 

The definitions of Eqs. (3.8) and (3.9) lead to the diagrammatic representation of the link 
operator shown in Fig. 2.1. 



In practice the lattice gluon field is defined at only one point on or nearby a link. 
This leads to interpolation errors in the integral in Eq. (3.9). For example, by choosing to 
evaluate the gluon field at x, the midpoint of the link, we have, 

1 /•«/2 

Ai = - dtAi{x + ti), (3.10) 

a J -a/2 

where i is a unit vector in the direction of link I. The interpolation errors can be obtained 
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by expanding in a Taylor series as follows: 

1 

a J -a 



1 /-a/a 
A, = - I dt 

i/2 



Ai{x) + tdiAi{x) + ^^edfAi^x) + 0(a3) 



= A,{x) + ^^djA{x) + r^dUi{x)... (3.11) 

We sec that the lattice gluon field reduces to its continuum counterpart in the a ^ 
limit, but that they differ by interpolation errors of order a?. From Eq. (3.11) we build 
the sequence of approximations to the lattice gluon field: 

Af\x) = A,{x) 

Af\x) = A{x) + ^a'dfA{x) (3.12) 
Af\x) = A{x) + ^a'dfAix) + ^a'dtAix) 



Having discussed classical errors we now move on to the more difficult topic of quantum 
errors in LGT. Quantum errors arise in LGT in two different contexts. Firstly, the lattice 
acts as an ultraviolet regulator allowing the simulation of only those states with momenta 

less than it /a. The absence of high momentum states results in a deviation between the 
lattice and continuum theories. Secondly, non-physical interactions arise due to the use of 
the link operator in constructing the lattice theory. To demonstrate this we expand the 
link operator in powers of e, 

C/^(x) = 1 + ieaAi,{x) - ^A^(x)A^(x) + . . . , (3.13) 

and note that the interaction of any number of gluons is allowed. Naively, the unphysical 
interactions are suppressed by powers of a. However, when contracted, products of pairs 
of gluon fields produce ultraviolet divergences (cx: 1/a^) which exactly cancel the a depen- 
dence of the expansion. These terms can be uncomfortably large and result in what are 
known as tadpole errors. 



In the last decade the improvement programme has led to a good understanding of 
both classical and quantum errors in quark and gluon actions (See Ref. [45] and refer- 
ences within). In contrast, only the lowest order classical errors have been corrected in 
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the Kogut-Susskind Hamiltonian. Conjectures have been made about the structure of a 
quantum improved Hamiltonian [46], but a perturbative study, which would be needed to 
fix its precise form, has not yet been carried out. 



We describe how tadpole improved lattice Hamiltonians can be constructed in Sec- 
tion 3.5. 



3.3.3 Gauge Inveiriance on the Lattice 



Under a local gauge transformation, A(a;), the link operator, Ui{x), transforms as follows, 

Ui{x) A{x)Ui{x)A^x + at). (3.14) 

Consider the closed path, C, on the lattice consisting of links {xi,ii), {xi+aii,i2), ■ ■ ■ , {xi — 
aijn,im)- The ordered product of link operators along C is given by 

P^{x) = Ui^{xi)Ui2{xi + ai2)...Ui^{xi-aim)- (3.15) 

Making use of Eq. (3.14), under a local gauge transformation Ti{x) we have 

Pf{x) ^ Ui,ixi)Pfix)ul{x^). (3.16) 

Since the trace of a matrix does not depend on the choice of representation (i.e. Tr{AU A'^) = 
TtU), by taking the trace of Pf{x) we obtain a gauge invariant operator. The traces of 
closed loops are known as Wilson loops and arc the building blocks of LGT. The simplest 
Wilson loop on the lattice is called the plaquette. It is the trace of the ordered product 
of link operators around a square of side length a. Using different combinations of Wilson 
loops one can construct lattice theories which have different approaches to the continuum 
limit. 
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3.3.4 Commutation Relations 

In continuum gauge theory, when employing canonical quantisation the following equal 
time commutation relations are postulated, 



Et{x),A]{y) 
'At{x),A){y) 



iSijSabS{x - y) 
Ef{x),EUy) 



0. 



(3.17) 



When constructing a lattice Hamiltonian the building blocks should be related in such a 
way that the correct continuum commutation relations are restored in the limit of zero 
lattice spacing. To construct the corresponding commutation relations on the lattice we 
need to use the lattice chromoelectric field, £, and the link operator, Ui. While we can 

determine the continuum behavior of the link operator with a simple Taylor expansion 
about a, the requirement that the correct continuum commutation relations be restored 
as a — > defines the relationship between the lattice and continuum chromoelectric fields. 

Let us now consider the commutation relations between the lattice chromoelectric field 
and the link operator. 



(3.18) 



Making use of the Campbell-Baker-Hausdorff formula for AT x AT matrices A and B, 



e-^Be^ = B + [B,A] + ^^ [[B, A],A] + ^^ [[[B, A] , A] , A] + . . . , (3.19) 



and Eq. (3.12) gives to lowest order in a 



^1 1 



lea 



tea 



'-'I ' •^m 



(3.20) 



Here m is the link {y,j)- In order to produce the correct continuum limit commutation 
relations, the lattice electric field must be proportional its continuum counterpart to lowest 
order in a; 



si'^^ = a[Ef{x) + 0{a')] 



(3.21) 
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Here / is the link (x, i) and a is a constant which we now determine. Substituting Eq. (3.21) 
in Eq. (3.20) and using the continuum commutation relations Eq. (3.17) gives 



c(0)a TT 



= teaa Ef{x), A'j{y)\^ X'Um + 0{a') 
= -eaaSijS{x - y)\^Um + C(a^) 



(3.22) 



In the last line we have related the discrete and continuous delta functions by 5xy ~ 

a'^5{x — y) in the small a limit. We set a = — a^~^/e so that the commutation relations 
between the lattice chromoelectric field and the link operator do not depend on either the 
coupling or the lattice spacing. This leads to, 



sf'^^ = -"^[Et{x)+0{a')\^ 



(3.23) 



and the conventional commutation relation, 



c(0)a TT 



(3.24) 



Proceeding exactly as above, or simply by pre and post-multiplying both sides of Eq. (3.24) 
we can show that the analogous commutation relation for Um is given by 



(3.25) 



3.3.5 The Kinetic Hamiltonian 



Using Eq. (3.23) we see that the kinetic part of a lattice Hamiltonian with order errors 
can be defined as follows 

^(0) ^ ^^ffff. (3.26) 
I 

By comparison with the kinetic part of the continuum Hamiltonian of Eq. (3.1), it is clear 
that Eq. (3.26) has the correct form in a ^ limit. To show this we substitute Eq. (3.23) 
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in Eq. (3.26) which leads to, 

I x,i 

« J d'^xTr{E^) +0{a'^). (3.27) 

Here, in the small a limit, we have replaced the sum over x hy a d dimensional integral. 
This is the correct continuum limit kinetic Hamiltonian up to order corrections. Since 
the Hamiltonian has dimensions of mass, from Eq. (3.26) we see that in d + 1 dimensions 
the coupling constant, e^, has units of (length)''"^. Thus, in 2+1 dimensions the coupling 
constant has dimensions of mass and in 3+1 dimensions it is dimensionless. We can define 
a dimensionless coupling by 

£,2 = e^a^-'^. (3.28) 
In terms of the dimensionless coupling the kinetic Hamiltonian is 

= (3.29) 



3.3.6 The Potential Term 



We have seen in Section 3.3.3 that gauge invariant operators on the lattice are built from 
the traces of products of link operators around closed loops. The simplest loop, the square 
of side length a is called the plaquette. Consider a plaquette which joins the lattice sites 
X, X + ai, X + a{i + j) and x + aj. We define the plaquette operator on such a plaquette 

by 



Pijix) 



1 - — ReTr 
N 



Ui{x)Uj{x + ai)Ul {x + aj)U\ {x 



(3.30) 



We will often use the equivalent pictorial expression to condense notation: 



(3.31) 
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The directed square denotes the traced ordered product of hnk operators, starting with 
Ui{x), around an elementary square, or plaquette, of the lattice. 



Q := Tr [Ui{x)Uj{x + ia)uj{x + ja)u]ix) 



(3.32) 



The reason for choosing the particular form of Eq. (3.31) for the plaquette operator, Pij, 
becomes obvious once we expand it in powers of a. Making use of the definition of the 
link operator (Eq. (3.8)) and Stokes' theorem we have. 



^ii(x) = 



1 — — ReTrexp I ie (b dx ■ A 
N \ In 



1 - —ReTrexp 
N 



■a/2 

le I ds 

a/2 



a/2 
-a/2 



dtFij{x + si + t3)+0{A'^) 



. (3.33) 



Expanding the exponential in the above and expanding Fy in a Taylor series in a leads to 
Pij{^) = ^^TT[F,,{xf]+^Tr[F,,ix){Df + D])F,,ix)]+Oia'). (3.34) 



2AN 



Here Di = di — ieAi is the gauge covariant derivative. The potential term of a lattice 
Hamiltonian with order corrections is thus given by 



x,i<j x,i<j 



(3.35) 



To see this we make use of Eq. (3.34) to give 



V(°) = a'^ J]TV[F,,-(a;)2]+C»(a2) 

x.i<j 

= J d'^xTr (B^) + 0{a^). 



(3.36) 



The last line follows simply from the definition of the chromomagnetic field in Eq. (3.2). 



Combining the kinetic and potential parts (Eqs. (3.29) and (3.35) respectively) leads 



32 



Constructing and Improving Lattice Hamiltonians 



to the Kogut-Susskind Hamiltonian 




(3.37) 



which has the correct continuum form up to order corrections. 

3.4 Classically Improved Hamiltonians 

3.4.1 Introduction 

In this section we derive a classically improved Hamiltonian directly by adding additional 
gauge invariant terms and fixing their coefficients in order to cancel errors. As a first step 
we aim to correct the classical order errors arising in the Kogut-Susskind Hamiltonian 
of Eq. (3.37). We treat the kinetic and potential terms of the Hamiltonian separately. We 
start with the potential term. 

3.4.2 Improving the Potential Term 

To improve the potential term of the Kogut-Susskind Hamiltonian, we follow the process 
of improving the Wilson action. We introduce the rectangle operator Rij{x) in the i-j 
plane (with the long side in the i direction), 



Here the Wilson loop joins the lattice sites x, x + 2ai, x + a{2i+j) and x + aj. Following 
the procedure of the previous section we can expand the rectangle operator in powers of 
the lattice spacing as follows 



Rijix) = -—TV [Fijixf] + -—Tr [Fijix)i4Df + D])Fijix)] + 0(a«). (3.39) 



We see that the expansion of Rij{x) gives precisely the same terms as the expansion of 
Pij{x) in Eq. (3.34), but weighted differently. This allows us to take a linear combination 
of Rij{x) and Pij{x) and fix the coefficients so that the order errors cancel. We write 




(3.38) 
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the improved potential term as the hnear combination 




(3.40) 



Substituting Eqs. (3.34) and (3.39) into Eq. (3.40) we have 



V(i) = o'^ ^ { (X + 4y) TV [Fij{xf] 



x,i<j 



+^ (X + lOY) Tr [Fij{x)iDf + D^)Fij{x)] } + ©(a^). (3.41) 



In order to obtain the correct continuum hmit we must have X + AY = 1 and to cancel 
the order error we require X + lOF = 0. Solving simultaneously gives X = 5/3 and 
y = — 1/6. This leads to the improved potential term 



In principle, the next lowest order classical errors could be corrected by including addi- 
tional, more complicated Wilson loops in the potential term. This has not been done 
because many additional diagrams are required to cancel the large number of order 
error terms. Since these errors are overwhelmed by order a^g^ quantum errors in the La- 
grangian approach, addressing quantum corrections in the Hamiltonian approach would 
seem to be of more immediate importance. 

3.4.3 Improving the Kinetic Term 

Constructing an improved kinetic Hamiltonian with a finite number of terms has proven to 
be a nontrivial exercise. Luo, Guo, Kroger and Schiitte demonstrated an interesting trade 
off when using either the transfer matrix or Legendre transformation methods to derive 
an improved Hamiltonian [46]. Both techniques require the starting point to be an im- 
proved action. When one starts from an improved action incorporating rectangular terms 
the resulting Hamiltonian has infinitely many terms and couples links which are arbitrar- 
ily far apart. To produce a Hamiltonian which couples only nearest neighbour links, it 
was found necessary to start from a carefully constructed highly non-local improved action. 

Here we demonstrate an alternative approach, similar in nature to the classical improve- 




(3.42) 



34 



Constructing and Improving Lattice Hamiltonians 



ment of the potential term in the previous section. One only needs to include additional 
gauge invariant terms with appropriate continuum behaviour in the kinetic Hamiltonian. 
The coefficients of the additional terms are chosen so that the order errors vanish. 



Perhaps the most important property of the electric field is that it generates group 
transformations. Mathematically, this translates to the lattice electric fields and link op- 
erators satisfying the commutation relations given by Eqs. (3.24) and (3.25). It is desirable 
for these to hold on the lattice for any degree of approximation. Let us consider what hap- 
pens to these commutation relations for the approximation labelled by the superscript (1) 
in Eq. (3.12). Making use of the Campbell-Baker- Hausdorff formula (Eq. (3.19)) we have, 



{l)a, 



(3.43) 



We observe that if the lattice electric field is taken to be related to the continuum electric 
field via Eq. (3.23), an order a? error arises in the commutation relation. To cancel this 
error we set 



-.d-l 



Et{x)--d^Et{x] 



We can take this to order by setting 

„d-l 



,(2)a 



[X) 



(3.44) 



(3.45) 



In this way a sequence of approximations to the lattice electric field in terms of the con- 
tinuum electric field can be constructed. 



Making use of these approximations we can analyse the classical errors arising in the 
kinetic Hamiltonian. To cancel these errors we take the approach of adding new terms 
and fixing their coefficients in order to cancel the order a? error. We have a great deal of 
freedom in choosing additional terms. They are restricted only by gauge invariance and 
the need for an appropriate continuum limit. 



To understand the construction of gauge invariant kinetic terms involving the lattice 
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electric field, it is important to recall that the lattice electric field transforms as follows 
under a local gauge transformation, A{x): 

Si{x) A{x)£i{x)A\x). (3.46) 

The response of the link operator to a gauge transformation is given in Eq. (3.14). Con- 
sequently, the next most complicated gauge invariant kinetic term we can construct (after 
Tr [fi(a;)^]) couples nearest neighbour electric fields: 



Tr 



£i{x)Ui{x)£i{x + cd)UJ{x) . (3.47) 



More complicated gauge invariant terms are easily constructed. One only needs to couple 
electric fields on different links anywhere around a closed loop. Consequently, generat- 
ing Hamiltonians with higher degrees of improvement would be seem to be more readily 
achieved within this approach. 



3.4.4 Continuum Limits 

In this section we examine the small a limit of gauge invariant lattice operators involving 
electric fields for possible use in the kinetic term of the Hamiltonian. 

We start with the simplest possible gauge invariant operator containing lattice electric 
fields, the kinetic term of the Kogut-Susskind Hamiltonian, 



2a 



q2^3 d 



(3.48) 



X,l 



X,l 



Making use of Eq. (3.44) to express the lattice electric fields in terms of their continuum 
counterparts gives 



2a 



x:(o) = 



,2d-2 



2e2 



E 

x,i 



Et{x)Et{x) - -Et{x)djEt{x) + 0{a') 



. (3.49) 



In the small a limit we can replace the sum over x with a d dimensional integral; 



x:(o) 



Et{x)Et{x) - -Et{x)d^Et{x) + 0{a') 



(3.50) 
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Eq. (3.50) quantifies the order discrepancy between the Kogut-Susskind kinetic Hamil- 
tonian and its continuum counterpart. 



We now move on to the next most compUcated term; one with nearest neighbour 
correlations. Let us consider the small a limit of the term 

A/C = ^Tr [fi(a;)C/i(a;)fi(a; + az)C//(a;) 

x,i 

= ^£^{x)S^{x + ai)Tr\x'^Ui{x)\^U^i{x)\. (3.51) 



Assuming an infinite lattice we can shift the summation variable x by any finite number 
of lattice sites without changing the result: 

1 



A/C = -^<^^f(x)^f(x + ai)Tr X''Ui{x)X^uj (x) 



e^{x)£^{x - ai)Tr \\<^Ui{x - ai)\'^uUx - ai) 



A"A^ + iea 



+£^{x)£^{x-ai)Ti A'^A" - iea A'^.A" Aiix - ai) j + 0{e^a^). {3. 52) 

Here we have used the definition of Ui (Eq. (3.8)) and the Campbell-Baker-Hausdorff 
identity (Eq. (3.19)). Making use of the totally antisymmetric SU(iV) structure constants, 
/"'^T, defined in Eq. (3.7), and the relation Ai{x) = Ai{x + at/2) + 0(^0^), we expand A/C 
in a keeping the lowest order terms to give 



x,i 



£^{x + ai) + £^{x - ai) 
+Tr [ea/°^TA^A(a;) + O(ea^) £^{x) £^{x + at) - £^{x - at) . (3.53) 



Here we have made use of the antisymmetry of the structure constants. We now expand 
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the lattice chromoelectric fields in powers of a to obtain 



+ea^rf^''A]{x)S^{x)diS^{x) + 0{ea\ a^) 



(3.54) 



The final step is to express the lattice electric fields in terms of their continuum coun- 
terparts. Let us work to order so that £ = £^^^ + O(a^). Substituting Eq. (3.44) in 
Eq. (3.54) gives 



2e2 



Er{x)Enx) + -a'Enx)dfErix) + 0{ea\ a') 



(3.55) 



It should be pointed out that while the Kogut-Susskind and nearest neighbour kinetic 
terms both produce the correct continuum limit with order classical corrections, the 
nearest neighbour term introduces a new quantum error of order ea^. We briefly discuss 
quantum errors in Section 3.7. 



3.4.5 The Improved Kinetic Term 

Having derived the small a expansions of the two simplest gauge invariant kinetic terms 
we can form a linear combination of them and choose the coefficients so that the correct 
continuum limit is obtained and the order classical error vanishes. Incorporating nearest 
neighbour interactions leads to the simplest improved kinetic Hamiltonian: 

/CW = ^ ^Tr yX£i{x)£i{x) + Y £i{x)Ui{x)£i{x + ai)C//(a;)J . (3.56) 



Substituting the small a expansions from Eqs. (3.50) and (3.55) into this linear combination 
we obtain 



^ + Et{xf+I^-^^] E,{x)dfE,{x)+0{ea\a'^ 



Ei{x)dtEi{x) + 0{ea\a'' 



. (3.57) 



To obtain the correct continuum limit we must set X + y = 1. To cancel the order 
error we require 5Y — X = 0. Solving these equations simultaneously gives X = 5/6 
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and Y = 1/6. Substituting back into Eq. (3.56) results in the order a? improved kinetic 
Hamiltonian, 



a ^ 



^Ei{x)Si{x) + ^Si{x)Ui{x)£i{x + a%)u\{x) 



(3.58) 



This is the result of Luo, Guo, Kroger and Schiitte [46] . We can take the this to order 
o by including next nearest neighbour interactions. A similar calculation using 
from Eq. (3.45) gives 



a ^ 



X,l 



97 1 - + 

—£i{x)£i{x) + -£i{x)Ui{x)£i{x + ai)Uj{x) 



-^£i{x)Uiix)Ui{x + ai)£i{x + 2ai)uJ{x + ai)uj{x) 



(3.59) 



Combining Eqs. (3.58) and (3.42) leads to the simplest classically improved lattice 
Hamiltonian devoid of 0{a^) errors, 



^£i{x)£i{x) + ^£i{x)Ui{x)£i{x + ai)uj{x) 



a ^ 

X,l 



(3.60) 



X,t<J 



3.5 Tadpole Improvement 

Tadpole improvement is an important step in producing more continuum-likc actions in 
LGT. Its implementation for pure glue on both isotropic and anisotropic lattices is well 
known in the Lagrangian approach. When tadpole improvement is included agreement 
between Monte Carlo and perturbative calculations on the lattice is achieved [45] . 



Tadpole improvement in the Lagrangian formulation is carried out at tree-level simply 
by replacing all link operators, U, by the tadpole improved link operator, U/uq, where 
uq is the mean link operator. The scheme was developed a decade ago by Lepage and 
Mackenzie [67]. 



In the Hamiltonian formulation the implementation of tadpole improvement in the 
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potential term follows the procedure in the action approach; all links are divided by the 
mean link. However for the kinetic term two conflicting implementations have been sug- 
gested. The earliest starts from a tadpole improved action and carries factors of uq into 
the Hamiltonian [46]. More recently it was suggested that no tadpole improvement was 
necessary in the kinetic term of the improved Hamiltonian [66]. Here we present our own 
views on the correct implementation. 



In the Hamiltonian approach the question of whether the lattice electric field should be 
rescaled under tadpole improvement arises. This question is easily answered by considering 
the commutation relations between the link operator and electric field: 



We observe that the lattice electric field cannot be rescaled and still maintain the correct 
commutation relations. Thus under tadpole improvement the electric field cannot change. 
We must, however, divide the second of the kinetic terms by a factor of u^. Tadpoles arise 
in this term because the chromoelectric and gluon fields do not commute. 



Including tadpole improvement in Eqs. (3.42) and (3.58) leads to the simplest order 
tadpole improved Hamiltonian: 



[£nx),Uj{y)] = 6ij6^yX''Ui{x). 



(3.61) 



We see that if we divide all link operators by uq we have 



[£^{x), —Uj{y)] = dij6^yX'^—Ui{x). 



(3.62) 




(3.63) 
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Here, the tadpole improved plaquette and rectangle operators are defined as follows 



3.6 Additional Improved Hamiltonians 



In this section we use the direct method to calculate improved Hamiltonians coupling ar- 
bitrarily distant lattice sites in both the kinetic and potential terms. Such Hamiltonians 
may be of use in calculations where vacuum wave functions with long distance correlations 
are used. Such wave functions are used in coupled cluster calculations of glucball masses. 
It is possible that more accurate calculations of high mass excited states may be obtained 
using such extended Hamiltonians. Again we consider the kinetic and potential terms 
separately. We start with the potential term. 



For the potential term we will couple lattice sites that are separated by, at most, m 
links. Consider the extended rectangle operator, 

(3.65) 

' 'rfi 

Here the rectangular Wilson loop is a 1 x m rectangle joining the lattice sites x, x + mai, 
X + a{mi + j) and x + aj. We can follow the procedure of the last section to calculate the 
small a expansion of Rm,ij{x). After some algebra, the procedure of Section 3.3.6 leads to 

Rm,i,{x) = '^:^Tr{F^.) + '^^Tr[F,,im'Df + D])F,,]+0{a'). (3.66) 

We can build an improved potential term as in the previous section by taking an ap- 
propriate linear combination of plaquette and extended rectangle operators. Let such an 
improved potential term be defined by 

2N ( Y 1 

^ x,i<j ^ ^ 



Rm,ij{x) = l-T7Re 
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Substituting the small a expansions from Eqs. (3.66) and (3.34) gives 

' TV [Fij {Df + D])Fij] + 0{a^) \ . (3.68) 



x,i<j 



e^a^ r m?{rin? + 1) 

X H ^i; -Y 



12 

We see that in order to obtain the correct continuum limit and cancel the order correc- 
tions we must have 



X + m^Y 



^^mW + 1)^ = 0. 



Solving these equations simultaneously results in 

-1-1 2 

X= and Y = ^-^ — -. (3.69) 

— 1 m'^{m'^ — 1) 

Substituting these values in Eq. (3.67) leads to the order improved potential term 

= ^ E {^^^^^(^) - ^2(^2 _ 1) [Rrn,j{x) + RmA^)]} ■ (3.70) 

We note that the result for the simplest improved potential term given by Eq. (3.42) is 
obtained by setting m = 2. 

We can construct an extended kinetic term in a similar fashion. In the kinetic term 
we couple lattice sites that are separated by, at most, n links. We start by considering the 
small a behaviour of the extended kinetic operator, 

A/CW = J2^[s,{x)U^^^^^^^,{x + nai)Ul^^^J. (3.71) 

03,1 

Here we have introduced the following notation for the path ordered product of link oper- 
ators joining sites x and x + nai, 

n-l 
1=0 

The small a expansion of AK^n^ can be determined using the procedure of Section 3.4.4. 
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After some algebra we arrive at 



X,l 



6n^ - 1 



(3.73) 



We define the extended improved kinetic Hamiltonian by taking the following linear combi- 
nation of the standard Kogut-Susskind kinetic term and the extended kinetic term A/Cn^^ : 



Again we fix the constants X and Y to produce the correct continuum result and cancel 
the order correction term. Substituting the small a expansions from Eqs. (3.49) and 
(3.73) into Eq. (3.74) gives the following: 



(6n2 - l)Y - X 



12 



a^Et{x)dlEt{x) . (3.75) 



To obtain the correct continuum limit and cancel the order error term we set 



X = 1- and Y = 



(3.76) 



Substituting in Eq. (3.74) we obtain the following extended improved kinetic term 



4^) = ^^TV 



1 - ^ ) £i{x)£i{x) 



£i{x + nai)U_ 



x^x+nai 



(3.77) 



We notice the result for the simplest improved kinetic term is recovered by setting n = 1. 



Combining Eqs. (3.77) and (3.70) leads to the classically improved Hamiltonian cou- 
pling lattice sites, at most, n-links apart in the kinetic term and, at most, m-links apart 
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in the potential term, 



a ^ 



1 - ^) £i{x)Ei{x) 



x,i<j ^ 



m?{m? — 1) 



[Rm,ij{x) + Rm,ji{x)\\ ■ (3-78) 



It should be pointed out here that Luo, Guo, Kroger and Schiitte derived an improved 
Hamiltonian [46] with a kinetic term which is a sum over n of A/Ci 
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Having calculated some examples of classically improved Hamiltonians we now move 
on to briefly discuss the more difficult topic of quantum improvement. 



3.7 Quantum Errors 

The correction of quantum errors has been instrumental in the progress of Lagrangian 
LGT. While classically improved actions are built from linear combinations of traced Wil- 
son loops with the constant coefficients chosen to cancel leading order errors, quantum 
improved actions determine the coefficients to one-loop order giving them a dependence. 
A method for determining these coefficients in weak coupling perturbation theory was de- 
veloped by Liischer and Weisz [65]. Unable to find sufficient conditions in the physical 
theory, Liischer and Weisz considered the situation where two dimensions are compactified 
and the gauge field obeys twisted periodic boundary conditions. This mechanism, famil- 
iar from Kaluza-Klein theories, gives the gluons mass and provides many simple on-shell 
quantities for the study of quantum improved actions. For the purpose of calculating the 
coefficients of a quantum improved action at the one-loop level, Liischer and Weisz used 
the masses of the asymptotic gluon states and some simple scattering amplitudes. The 
result of their SU(3) calculation in 3-1-1 dimensions, with errors starting at 0{g^a^), is 
termed the Liischer- Weisz improved action. 

Sun = 

r 



+ 0.2370C/2 ' '^^ ~ ' 12 ^ 0.025215^ ) £,2 - ^m^^Xg^ 



(3.79) 
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with 




(3.80) 



The classically improved action is recovered from Eq. (3.79) by modifying the coefficients 
so that only the order terms remain. It is important to note that since the coefficient 
of £3 is small, a quantum improved action is obtained, to good approximation, by keeping 
only plaquette and rectangular terms. The tadpole improvement of the Liischer-Weisz 
action at the one-loop level, which we do not discuss here, is straightforward [68] . A con- 
vincing demonstration of improvement using the tadpole improved Liischer-Weisz action 
was first shown in the restoration of rotational symmetry in the qq potential [69] . 

In the Hamiltonian formulation quantum errors are yet to be considered in detail. It 
would seem that their treatment would require considerable effort. Luo, Guo, Kroger and 
Schiitte have suggested a form for a quantum improved Hamiltonian by taking the time-like 
lattice spacing to zero in a quantum improved Liischer-Weisz action, and then performing 
a Legendre transformation [46] . To extend their suggestion to a precise form would require 
new perturbative calculations in the Lagrangian approach in the at <C limit^, possibly 
following the techniques of Lxischer and Weisz described above. The remaining step would 
be to relate the Lagrangian and Hamiltonian lattice couplings. This could be done using 
the background field method, extending the work of Hamer for the unimproved case [61], 
to the improved case. 

3.8 Conclusion 

In this chapter we have introduced a direct method for the construction of improved 
lattice Hamiltonians. This method differs from Legendre transform and transfer matrix 
techniques in that one does not need to start with a lattice Lagrangian. One simply needs 



^Here at is the time-like and O3 the space- Uke lattice spacing. 
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to construct suitable linear combinations of gauge invariant terms and fix their coefficients 
so that the correct continuum limit is obtained to the desired level of approximation. 
The advantage of the direct approach is that it can be easily applied in the construc- 
tion of various improved Hamiltonians, as was demonstrated in Sections 3.4 and 3.6. We 
have also discussed the implementation of tadpole improvement within the Hamiltonian 
formulation and how the construction of quantum improved Hamiltonians may be possible. 

The extension of the direct method to the construction of lattice Hamiltonians for 
quarks and gluons is straightforward. One simply needs to construct appropriate gauge 
invariant quark operators and form linear combinations of them such that the correct con- 
tinuum Hamiltonian is obtained. The construction of quantum improved Hamiltonians 
is of more immediate interest but their calculation seems to require significant additional 
work. 

Having constructed a variety of lattice Hamiltonians we now move on testing the levels 
of improvement that they offer. In the next chapter we introduce some simple tests of 
improvement. In Chapter 6 we continue the discussion of improvement by employing 
classically improved and tadpole improved Hamiltonians in the calculation of glueball 
masses. 
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Chapter 4 

Testing Improvement 



4.1 Introduction 

When performing a calculation on the lattice we are interested in the continuum limit in 
which the lattice spacing is taken to zero. It is only in this limit that physical results 
are obtained. For an improved theory one would expect the presence of lattice artifacts 
to be less pronounced, taking the lattice theory closer to continuum physics for a given 
lattice spacing. Whether or not this is the case for the improved Hamiltonians generated 
in Chapter 3 will become evident in specific calculations. 

In this chapter we perform some relatively simple checks on improvement. The first 
check presented is for the simple case of U(l) gauge theory. In this case we perform a per- 
turbative calculation of the static quark-antiquark potential in the strong coupling regime 
for both the Kogut-Susskind and improved Hamiltonians. One would expect an improved 
Hamiltonian to exhibit a break from the strong coupling limit at a larger coupling, cor- 
responding to a larger lattice spacing, than the unimproved case. We demonstrate in 
Section 4.2 that this is indeed the case. 

The second check presented is a variational calculation of the lattice specific heat for 
the case of SU(2) in 2+1 dimensions. For this case analytic results for the integrals involved 
are known and have been used in variational calculations for nearly 20 years. The peak 
in the lattice specific heat has been shown to correspond to the transition region from 
strong to weak coupling [70] . One would expect the peak of the improved lattice specific 
heat to be appear at a larger coupling (corresponding to a larger lattice spacing) than its 
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unimproved counterpart. We demonstrate in Section 4.3 that this is in fact the case. 



Without resorting to a detailed computation, a relatively straightforward check on im- 
provement can be made in a perturbative calculation of the static quark-antiquark po- 
tential in the strong coupling regime. Wc adopt the simplistic model of Kogut and 
Susskind [12] and consider a source and a sink of colour flux, each infinitely heavy, sepa- 
rated by r lattice sites. The potential existing between the quark and antiquark is described 
in terms of a string of excited links joining the particles. 

Consider the case of U(l) on a 3 dimensional spatial lattice with L links. The Kogut- 
Susskind Hamiltonian is given by 



In the strong coupling limit {g 2> 1) only the kinetic term survives. In this limit the po- 
tential term can be treated as a perturbation whose effect can be handled with standard 
Rayleigh-Schrodinger perturbation theory. Before proceeding with the calculation, we first 
introduce the Fock space of link excitations in which we intend to work. 

We define the lattice vacuum by the direct product of all strong coupling link vacuua: 



Here |0)j denotes the strong coupling vacuum of the i-th link upon which the action of the 
lattice electric field gives Si\0)i = 0. An excited link is defined by acting a link variable on 
the vacuum. For example: 



4.2 The U(l) Static Quark-Antiquark Potential 




(4.1) 



|0) = |0)i®|0)2®---®|0)l. 



(4.2) 



\l)i = Ui\0)i, \2)i = UiUi\0)i, etc. 



(4.3) 



Here the index I is equivalent to the pair of indices x (starting point) and i (direction) 
which label the link in the notation Ui{x). Employing the commutation relations between 
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S and U which, for the gauge group U(l), are given by 

[SuUm]=5imUi, (4.4) 

we sec that excited hnk states arc cigcnstates of the strong couphng Hamiltonian. In 
particular the energy associated with a singly excited link is /2a, and with a chain of r 
links rg^ /2a. 



We now return to the calculation of the static quark-antiquark potential. The lowest 
energy state is given by the direct line of flux joining the source and sink with all links 
singly excited. Its energy is given by 

V-r^. (4.5, 

Making use of Rayleigh-Schrodinger perturbation theory to include the effects of the po- 
tential term, the lowest order shifts in energy start at order g~^. Taking into account 
all possible contributions yields the static quark-antiquark potential in 3-1-1 dimensions 
(relative to the perturbed vacuum), 



2a \ 3g^ 



In 2+1 dimensions the perturbation is half the 3-|-l dimensional result. This is a result of 
there being only two possible orientations for a plaquette sharing a link with a line of flux 
in 2-1-1 dimensions, whereas in 3-1-1 dimensions there are four possible orientations. The 
2-1-1 dimensional result is then 

.(0) _ rg^ (, 1 



We can perform the same calculation with the classically improved Hamiltonian, which 
for U(l) is given by 



2a ^ 



x.i '- 



5 1 

-£i{xf + -£i{x)£i{x + ai) 



aq 
^ x,i<3 
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The improved strong coupling ground state energy of the chain of A'^ excited links is given 
by: 



"2a 



5 1 , 



11 



2a 



(4.9) 



The difference between this and the unimproved result arises because the improvement 
term, 3. ^^i(a;)^i(a3 + ai), gives non-zero contributions only when acting on pairs of 
adjacent excited links. Consequently, in a straight chain of r links there are only r — 1 con- 
tributions from this term. This difference vanishes however in the limit of large separation. 



The calculation of the lowest order perturbation in the improved case is considerably 
more complicated than the unimproved case. Including rectangular Wilson loops in the 
potential term results in a large number of contributions to the order shift. Collecting 
all contributions leads to the improved static quark-antiquark potential (relative to the 
perturbed vacuum) , 

<).!^ri-^UH^. (4.10) 

The 2+1 dimensional result can be derived as above: 

..(1) rg' A l-8825\ 0.59219 

^2+^ ^ 2aV 9' J ag^ ' ^^"^^^ 

We notice a large change in the magnitude of the first order perturbation in going from 
the Kogut-Susskind to the improved lattice Hamiltonian. In the limit of large separations 
we ignore the trailing order g~^ terms which are independent of r. 



The string tension, a, is defined as the coeflacient of the linear part of the quark- 
antiquark potential. It is a physical quantity and is measurable, in principle, in deep 
inelastic scattering experiments. Prom Eqs. (4.7) and (4.10) we obtain the Kogut-Susskind 
and classically improved U(l) string tensions in 3+1 dimensions respectively, 



r(0) = ^ 



2 / 2 



(4-12) 



2a2 V 

Since it is a physical quantity, the string tension must have a finite continuum limit. For 
this to be the case the coupling, g, must depend on the lattice spacing. Precisely how is 
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determined by removing the dependence of the string tension on a by setting 

d , , , , da dq da , . 

-.(<..9(a))=0=^ + ^^. (4.14) 

This is known as the renormahsation group equation. It is usual to describe the dependence 
of the couphng on the lattice spacing with the Callan-Symanzik (3 function, which, for 
Abelian gauge theories, we define by 

The renormahsation group equation then allows us to write the /? function in terms of the 
string tension as follows 

= (416) 

With this result we can use Eqs. (4.12) and (4.13) to obtain the Kogut-Susskind and 
improved (3 functions shown in Fig. 4.2, 

/?W , 15.0604 
-— ^ 1 (4.18) 

We can use this result to check the level of improvement offered by the improved Hamil- 
tonian. To do this we examine the break away from the strong coupling limit, —13/g = 1. 
One would expect a calculation using an improved Hamiltonian to break away from the 
strong coupling limit at a larger lattice spacing (or equivalently a larger coupling) than 
one using an unimproved Hamiltonian. We see from Fig. 4.2 that this is indeed the case. 



To be more precise, suppose the break from the strong coupling limit occurs at a 
coupling g^"^ for the Kogut-Susskind Hamiltonian and g^"^ for the improved Hamiltonian. 
A relationship between g^^ and g^^ can be deduced using Eqs. (4.17) and (4.18): 

.i"«(^^i^)%f«1.24,<»>. (4.19) 
The same result is obtained, to this order, in 2+1 dimensions. 
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Figure 4.1: The Callan-Symanzik (3 function versus coupling for U(l) in 3+1 dimensions. 

Having observed improvement in a simple strong coupling calculation for pure U(l) 

gauge theory wc now move on to a more detailed check on improvement away from the 
strong coupling limit. 

4.3 SU(2) Lattice Specific Heat in 2+1 Dimensions 

4.3.1 Introduction 

In this section we calculate lattice specific heats using various improved SU(2) Hamilto- 
nians in 2+1 dimensions. We start with the classically improved Hamiltonian given in 
Eq. (3.78), which couples links, at most, n lattice sites apart in the kinetic term, and m 

sites apart in the potential term. Following that wc consider tadpole improvement, dis- 
cussing its practical implementation in the Hamiltonian approach. Finally, we calculate 
the lattice specific heat with the order classically improved kinetic term. A discussion 
of the results follows. 

4.3.2 Preliminciries 



Before proceeding with the calculation we need to make some preliminary definitions. 
Consider for the time being SU(iV) pure gauge theory on a lattice with Np plaquettes. 
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The lattice specific heat is defined by 

Cv = (4.20) 

where (3 = N/g^ is the inverse coupling, not to be confused with the Callan-Symanzik 
(3 function of Section 4.2. Prom here on f3 will always denote the inverse coupling. The 
vacuum energy density, cq, is given by 

a (0o|7^|0o) , . 

We compute the vacuum energy density variationally, using the one-plaquette trial state, 

exp I c ^ ReQ j |0). (4.22) 

x,i<j J 

Here |0) is the SU(A'") analogue of the strong couphng vacuum defined in Eq. (4.2). At this 
stage Eq. (4.21) defines the vacuum energy density as a function of the coupling and the 
variational parameter, c. The variational parameter is fixed as a function of the coupling by 
minimising the vacuum energy density at each coupling. Using the commutation relations 
between and Um in Eqs. (3.24) and (3.25) the vacuum energy density can be expressed in 
terms of the expectation values of the Wilson loops that appear in the Hamiltonian under 
consideration. This can be done easily for all and any number of dimensions. The 
difficulty lies in the calculation of the expectation values themselves. There is no difficulty 
however for the special case of SU(2) in 2+1 dimensions for which analytic results are 
available. 

With the necessary preliminary definitions made, we now proceed to the calculation of 
lattice specific heats using a selection of the Hamiltonians derived in Chapter 3. 

4.3.3 The Extended Hamiltonian 

In this section we calculate lattice specific heats for SU(2) LGT in 2+1 dimensions using 
the extended improved Hamiltonian Tinln of Eq. (3.78). Let us denote the lattice specific 
heat calculated with Hnln by 



{n,m) 



(4.23) 
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where the vacuum energy density is given by 



^(n,m) _ g {M'Hn)ni\<t>o) 

° Up (^o|0o) 



(4.24) 



Here |^o) is the one-plaquette trial state of Eq. (4.22). Making use of Eq. (A. 6), after 
some algebra, ep*'"*^ can be expressed in terms of the expectation values of plaquettes and 
extended rectangles as follows: 



1 \ iV2 - 1 2/3 m2 + 1 
1 - I — — — c 



Qv? J 2/3 

4/3 1 



N m2 



+ 



A'' m?{m? — 1) 



r]{D 

'3 



+ 2/3 1 + 



(4.25) 



For the special case of SU(2) in 2+1 dimensions, the expectation values appearing in 
Eq. (4.25) can be evaluated analytically [1, 30], with the result. 



„ \ (<^o| P] |</'o) ^ ^ 



m 



3) -Mo 



/2(4c) 
/i(4c) 



(4.26) 



Here /„ is the n-th order modified Bessel function of the first kind defined for integers n 
by 



In{2x) = 



2^ ^2k+n 



k=0 



k\{k + n) 



(4.27) 



Making use of Eq. (4.26) in Eq. (4.25) gives the SU(2) improved vacuum energy density 
in 2+1 dimensions. 



(n,m) _ 
^0 — 



(3 ^ m? — 1 



J2(4c) 
/i(4c) 



+ 



_4^ / /2(4c) 



m?{m? — 1) \Ii(4c) 



2/3 1 + ^ 



(4.28) 
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The result of Eq. (4.28) should be compared with the unimproved result of Arisue, Kato 
and Fujiwara [1], derived from the Kogut-Susskind Hamiltonian: 



With an analytic expression for the energy density as a function of (3 in hand, the 
lattice specific heat may be calculated using Eq. (4.23). We now consider the effect of 
varying the extent of the correlations in the kinetic term and potential term separately. 
We start with the kinetic term. 

Fig. 4.2 shows a plot of the lattice specific heats Cy''^\ Cy''^^ and Cy^^'^"* (those with 
2x1 rectangles in the potential term), as well as the unimproved result. The location of the 
peak, /3*, indicates the crossover region between the strong and weak coupling regimes [70]. 
The graphs exhibit the behaviour expected from an improved Hamiltonian. It is clear that 
for the improved Hamiltonians the transition from strong to weak coupling occurs at a 
larger coupling (smaller /3), and hence a larger lattice spacing, than the unimproved case. 
The approximate values of /3* are: 0.8850 for the unimproved case, 0.5955 (n = 1), 0.6391 
(n = 2) and 0.6533 (n = 100). This implies that when using an improved Hamiltonian 
one is closer to continuum physics when working at a given lattice spacing. It should be 
noted that the degree of improvement, does not vary greatly with n. This is expected, 
since distant correlations in the kinetic term will not contribute to the energy density 
unless equally distant correlations arc included in the trial state. For such complicated 
trial states, Monte Carlo techniques are necessary to calculate the required expectation 
values. Such calculations will not be presented here. They are made difficult by very 
shallow minima in the vacuum energy density [49] . 

We now examine the dependence of improvement on the potential term. Fig. 4.3 shows 
a plot of the lattice specific heats Cy''^\ Cy '"^^ and Cy'^\ as well as the unimproved result. 
Again, it is clear that for the improved Hamiltonians, the transition from strong to weak 
coupling occurs at a larger coupling (smaller /3), and hence a larger lattice spacing, than 
for the unimproved case. The values of (3* are: 0.5955 (m = 2), 0.7133 (m = 3) and 0.7768 
(m = 5). The degree of improvement is best for m = 2 and becomes worse as m increases. 
This is to be expected since the errors present in Eq. (3.68) are 0{m^a^), which is large for 




(4.29) 
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Figure 4.2: The 2+1 dimensional SU(2) lattice specific heat calculated with various Hamil- 
tonians. Cy'^\ C^'^^ and Cy'^^''^^ are improved results derived from Hamiltonians with 
varying kinetic terms. The unimproved result refers to the calculation of Arisue, Kato and 
Fujiwara [1]. 

large m. This suggests that there may be limits on the utility of improved Hamiltonians 
which couple distant lattice sites. It is clear that the degree of improvement is, by far, 
more sensitive to varying the potential term than varying the kinetic term with the chosen 
trial state. 



Having examined classical improvement in the context of improved Hamiltonians cou- 
pling distant lattice sites, we now move on to consider the effect of tadpole improvement 
on the location of the lattice specific heat peak. 

4.3.4 Tadpole Improvement 

In this section we extend the work of Section 4.3.3 to consider tadpole improved SU(2) 
LGT in 2+1 dimensions. As discussed in Section 3.5, tadpole improvement is an important 
step in the improvement of LGT. Tadpole errors arise when expectation values involving 
products of link operators are taken. When contracted, products of gluon fields, Ai{x), 
produce ultraviolet divergences. Such divergences spoil small a perturbation series, leading 
to inaccurate results in the continuum limit. Lepage and Mackenzie demonstrated that 
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Figure 4.3: The 2+1 dimensional SU(2) lattice specific heat calculated with various Hamil- 
tonians. Cy''^\ Cy'^^ and Cy'^^ are improved results derived from Hamiltonians with 
varying potential terms. 

tadpole errors are largely cancelled by dividing each link operator by the mean link uq [67] . 

In this section we calculate the lattice specific heat using the simplest tadpole improved 
Hamiltonian given by Eq. (4.3.3). We define the SU(iV) mean link in terms of the mean 
plaquette as follows: 

In the Lagrangian formulation, the more difficult choice of fixing to the Landau gauge and 
calculating the mean link directly is often used. We choose to define the mean link in 
terms of the mean plaquette for simplicity. In this case one does not need to perform a 
gauge fixing, which would be difficult in the Hamiltonian approach. 

To incorporate tadpole improvement in Hn}m is straightforward. We follow the pre- 
scription discussed in Section 3.5 in which all links are divided by the mean link and the 
lattice electric field operators remain unchanged. Making use of this prescription leads to 
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the tadpole improved version of Hnju, 



'■n,m 



2a ^ 



1 - 



Q^^^£i{x)U3:^x+nai£i{x + nai)ul_ 
2N (m? + 1 1 



X,l<J 



(4.31) 



where the tadpole improved plaquette operator is given by Eq. (3.64) and the tadpole 
improved extended rectangle operator is defined by, 



Rij{x) 



2m+2 




Re 



(4.32) 



m 



Here the mean link is given by Eq. (4.30). Proceeding as in Section 4.3.3 we obtain the 
tadpole improved energy density, 



^(n,m) 
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, _ M 3 
6n2 J 
4/3 



2/3 + 1 
~ 1 



/i(4c) 



2m+2 




m?{m? - 1) Vi(4c) 



+ 2/3 1 + 



(4.33) 



Incorporating tadpole improvement in a variational calculation is not straightforward. 
This is because the mean plaquette depends on the variational state, which is determined 
by minimising the energy density. The energy density however depends on the mean pla- 
quette. Such interdependence suggests the use of an iterative procedure for the calculation 
of the energy density. The approach we adopt is as follows. For a given (3 and starting 
value of uo we minimise the energy density of Eq. (4.33) to fix the variational state, |(^o)- 
We then calculate a new mean plaquette value using this trial state and substitute it in 
Eq. (4.33) to obtain a new expression for the energy density which is then minimised. This 
process is iterated until convergence is achieved, typically between five and ten iterations. 
The result of the procedure for the simple case of n = 1, m = 2 is shown in Fig. 4.4. For 
the purpose of comparison we also plot the unimproved lattice specific heat and the most 
improved result of Section 4.3.3. We observe a greater level of improvement for the tadpole 
improved case (/3* = 0.5565) than the order classically improved case (/3* = 0.5955). 



We now move on to consider the effect of incorporating order classical improvement 
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0.0 0.5 1.0 1.5 2.0 

Figure 4.4: The tadpole improved 2+1 dimensional SU(2) lattice specific heat compared 
with the standard improved and unimproved results. 

in the kinetic Hamiltonian. 



4.3.5 Classical Order Improvement 

In this section we calculate the lattice specific heat using a Hamiltonian with an or- 
der classically improved kinetic term with next nearest neighbour interactions given 
in Eq. (3.59). For the potential term we use the standard order improved result of 
Eq. (3.42). Ideally we would use an order improved potential term so that all terms 
were improved consistently to the same order. In principle, such a potential term could be 
constructed by adding appropriately chosen additional, more complicated operators to the 
order potential term of Eq. (3.42). This will not be done here because many additional 
operators are needed to cancel the large number of contributions to the order error. 
Consequently in this section we are not able to compare an order classically improved 
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calculation with its order c? equivalent. Instead we examine the effect of incorporat- 
ing easily calculated addtitional improvements in the kinetic term of the Kogut-Susskind 
Hamiltonian. 



Tadpole improvement is easily incorporated in the Hamiltonian described in the pre- 
ceding paragraph by dividing all links by the mean link and leaving the lattice electric field 
operators unchanged as described in Section 3.5. Making use of this prescription leads to 
the Hamiltonian 



7^(2,1) = ^yiv 



97 1 ' t 

-^ei{x)ei{x) -^£i{x)Ui{x)£i{x + ai)Ul{x) 

^ Si{x)Ui{x)Ui(x + ai)£i{x + 2ai)U^{x + ai)U^{x) 



+^ E {l^^ji^) - ^ i^ji^) + Rd^)]] ' (4-34) 

" x,i<j ^ ^ 

where the tadpole improved plaquette and rectangle operators are given by Eq. (3.64). 
With this Hamiltonian we obtain the following tadpole improved SU(2) energy density in 
2-1-1 dimensions: 

_/97c 10/3\l2(4c) P 

The improved and tadpole improved lattice specific heats for this case are shown in Fig. 4.5 
along with the results of Section 4.3.4 for comparison. Without tadpole improvement we 
observe that the case of order improvement in the kinetic term (/3* = 0.5866) displays 
marginal improvement over the order improved case (/?* = 0.5955). Similarly small 
levels of improvement are achieved when tadpole improvement is included (/3* = 0.5475) 
compared to the order classically improved result of Section 4.3.3 (/3* = 0.5565). As 
mentioned in Section 4.3.3, such marginal improvement is expected, since the one-plaquette 
trial state is not particularly sensitive to distant correlations in the kinetic term. A trial 
state more sensitive to distant correlations in the kinetic term would require larger Wilson 
loops in its exponent. 
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0.0 0.5 1.0 1.5 2.0 

Figure 4.5: 0{a'^) classical and tadpole improved 2+1 dimensional SU(2) lattice specific 
heats compared with the standard improved and unimproved results. 

4.3.6 Conclusion 

We have shown that the level of improvement offered by an improved Hamiltonian does 
not depend significantly on the degree of separation of the sites coupled in the kinetic term 
when using the one-plaquette trial state of Eq. (4.22). The onset of weak coupling has 
been found to occur at /3 0.6 for each improved case, compared with P ^ 1 for the unim- 
proved case. The tadpole improved Hamiltonian with order improvement in the kinetic 
term produced a lattice specific heat peak at the largest coupling which suggests that this 
Hamiltonian exhibits the largest degree of improvement of those considered in this chapter. 

The weak dependence on distant correlations in the kinetic Hamiltonian is a result of 
only including plaquettes in the trial state. For such a choice the calculation of all matrix 
elements encountered is analytic but the improvement term does not explicitly contribute. 
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This is because the improvement term must act on an excitation of two neighbouring hnks 
to produce a nonzero result. Such a configuration does not occur with the plaquette trial 
state. If extended Wilson loops were included in the trial state exponent one would expect 
significantly different results for the dependence on separation in the kinetic term. One 
would expect the dependence of the lattice specific heat peak on the separation in the 
kinetic term to be significant until it matched the degree of separation in the trial state. 
Beyond that one would again expect a weak dependence. 

To extend the calculation to more complicated trial states requires a change in tech- 
nique. The required matrix elements must be calculated using Monte Carlo techniques as 
no analytic results are available. 

4.4 Summary 

In this chapter we have performed some simple calculations designed to test the level of 
improvement offered by the improved Hamiltonians constructed in Chapter 3. In both a 
strong coupling calculation of the U(l) string tension in 3+1 dimensions and a variational 
calculation of the SU(2) lattice specific heat in 2+1 dimensions we have observed that 

improved Hamiltonians do in fact provide improvement. By improvement wc mean that 
when working at a given lattice spacing (or coupling) we are closer to the continuum limit 
when working with an improved Hamiltonian. 

In Chapter 6 we return to the discussion of improvement. There we make use of 
the simplest classical and tadpole improved Hamiltonians, as well as the original Kogut- 
Susskind Hamiltonian, in calculations of SU(iV) glueball masses in 2+1 dimensions. Before 

presenting these calculations it is necessary to introduce the analytic techniques upon 
which the calculations rely. These techniques do nothing more than allow the calculation 
of expectation values of Wilson loops in 2+1 dimensions for general SU(A'^). 



Chapter 5 



Analytic Techniques in 
Hamiltonian Lattice Gauge 
Theory 

5.1 Introduction 

In Chapter 4 it was shown that when working with the one-plaquette trial state of 
Eq. (4.22) the vacuum energy density can be expressed in terms of expectation values 
of plaquettes and rectangles. In principle, this allows the variational parameter, c, to 
be fixed as a function of /3 in any number of dimensions and for any gauge group. The 
difficulty lies in the calculation of the plaqucttc and rectangle expectation values. Monte 
Carlo simulations can be used to calculate these expectation values in any number of di- 
mensions and for any gauge group. However, in 2+1 dimensions it is possible to calculate 
such expectation values analytically. 

In this chapter we develop techniques which allow the analytic calculation of various 
expectation values in 2+1 dimensions for SU(iV) gauge theory, extending the long known 
analytic results for SU(2) to SU(A^). Analytic expressions are obtainable in 2+1 dimensions 
because the change of variables from links to plaquettes has unit Jacobian [71]. Conse- 
quently, the plaquettes on a 2 dimensional spatial lattice are independent variables. For 
the case of 3+1 dimensions the same change of variables produces a complicated Jacobian, 
which seems only manageable on special lattices. The problems faced in 3+1 dimensional 
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Hamiltonian LGT are discussed in Chapter 8. 

Our motivation for considering SU(iV) expectation values is in their use for calculating 
general S\J{N) glueball masses and in particular examining their large N limit. This limit 
is of current interest due to a proposed correspondence between certain string theories 
and large N gauge theory. This correspondence was first conjectured by Maldacena for 
supersymmetric gauge theories [72] and later extended by Witten to pure SU(A'') gauge 
theories [73]. The Maldacena conjecture and the role of large N glueball masses are dis- 
cussed further in Section 7.2. 

The outline of this chapter is as follows. In Section 5.2 we consider the special cases of 
SU(2), for which analytic results have been known and used for many years, and SU(3), 
for which we rediscover a useful generating function. In Section 5.3 we calculate three 
generating functions for general SU(A/') which provide analytic results for the expectation 
values that appear in calculations of glueball masses. We finish the chapter in Section 5.4 
with a summary of the results obtained. 



In this section we discuss analytic techniques for use in the calculation of group integrals 
for the special cases of SU(2) and SU(3). Before deriving the generating functions which 

are useful in the calculation of SU(2) and SU(3) expectation values, the concept of group 
integration needs to be introduced. To do this wc consider the expectation value of a 
plaquette with respect to the one plaquette trial state of Eq. (4.22). We can write this 
expectation value as an SU(A/") group integral as follows: 



Here the products over / and /' extend over all links on the lattice, while the products 
over p and p' extend over all plaquettes on the lattice. We introduce the notation Zi(p) 
to denote the trace of plaquette p. For each integral in Eq. (5.1) the integration measure 
is given by the Haar measure (also called the invariant measure and less commonly the 
Hurwitz measure) [74, 63]. For any compact group G, the Haar measure is the unique 
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(5.1) 



5.2 Analytic Techniques for SU(2) and SU(3) 



65 



measure dU on G which is left and right invariant, 

/ dUf{U) = [ dUf{VU) = [ dUfiUV) VI/ G G (5.2) 

Jg Jg Jg 



and normahsed, 

/ dU = 1. (5.3) 
Jg 

In Eq. (5.2) / is an arbitrary function over G. 

In 2+1 dimensions the variables in Eq. (5.1) can be changed from links to plaquettes 
with unit Jacobian [71]. The plaquettes then become independent variables allowing the 
cancellation of all but one group integral in Eq. (5.1). All that remains is 



a 



Here, y is a plaquettc variable with TrV = Zi{0). For the case of SU(2), analytic expres- 
sions for the plaquette expectation value in terms of modified Bessel functions (Eq. (4.26)) 
have been used in variational calculations for almost 20 years. 



Integrals over SU(2) are easily handled using the parameterisation of SU(2) [74], 

(cos - + z sin - cos 6 i sin - sin de~^'^ \ 

2 2 2 VC/ G SU(2), (5.5) 

i sin - sin 6e^'^ cos - — i sin - cos 9 I 

2 2 2 / 

where the parameters p, 6 and (f) lie in the ranges < ^ < tt and < p,(f) < 27r. The Haar 
measure corresponding to this parameterisation is 

dU = —r sin^ ^ sin 9dpd(f)d0. (5.6) 
47r^ 2 

With this information it is straightforward to show, 

/ dUe'^'^^ = -/i(2c). (5.7) 

JSV{2) C 

All SU(2) integrals of interest to us can be obtained from Eq. (5.7) by differentiating 
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appropriately with respect to c. For example 



/ dUTiUe'^^ = -l2(2c), (5.8) 

JSV{2) C 

/ dC/(TrC/)V'^^ = -/i(2c) ^ -|/2(2c), etc... (5.9) 

isU(2) C C 

Here we have used the following property of the modified Bessel function under differen- 
tiation: 

dl Ti 

-^{x) = /nil ± -In{x). (5.10) 



The corresponding SU(3) results, which have not been used, to our knowledge, in 
Hamiltonian LGT, follow simply from a paper of Eriksson, Svartholm and Skagerstam [75], 
in which the SU(3) integral, 

[ d[/e^(^^'+^'^), (5.11) 
isu(3) 

is calculated for arbitrary 3x3 matrices M. Following their analysis, but treating M and 
as independent variables, with M = cl, = dl, and 1 the 3x3 unit matrix, leads 
to an expression for the SU(3) generating function: 

y{c,d) = f dC/e^('^^+<^^') 
^SU(3) 

Here F is a closed contour in the complex plane including the pole at z = but excluding 
the pole at z = cd and Jn is the n-th order Bessel function of the first kind defined, for 
integers n, by 

^ ( _']\k 2k+n 

To evaluate the contour integral in Eq. (5.12) we expand the integrand in power series 
about the pole at z = and use Cauchy's integral theorem to eventually obtain the 
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following convergent series: 



y{c,d) 



>(fc + l)!(fe + 2)!^l 



This generating functional is extremely useful in the context of Hamiltonian LGT. In 

principle, it permits an analytic investigation of 2+1 dimensional pure SU(3) gauge the- 
ory. The calculation of various matrix elements for all couplings, reduces to a exercise in 
differentiation. For example: 



{Z^ip)) - 
{Z^{p)Z^{p)) ^ (jo] 



ydcdd 



, and 



d=c 



1 gm+n-y 



(5.15) 



d=c 



As an application of Eq. (5.14) we can calculate the strong coupling limit (c = 0) of 
{[Zi(p)Y^ [Zi(p)]"). These integrals arise in coupled cluster calculations [76]. By differen- 
tiating Eq. (5.14) appropriately, we observe that non-zero strong coupling matrix elements 
occur only when n + 2m = OmodS. For this case we have. 



Sk+3 



2nlm.l 



^n+m+S {k + l)\{k + 2)\{r^-k)\{ 



m+2n 



-k)\ 



,(5.16) 



where the sum runs over all integers ^^^^ < k < min ( ^ m+2n ^ _ fpj^jg g^j-ong coupling 
result has an equivalent combinatorial interpretation as the number of times the singlet 
representation appears in the direct product of m 3 and n 3 representations [63] . Although 
formulas for the case of m = exist [77, 78], we have not seen a general expression for 
m ^ n published elsewhere. These numbers are well known in combinatorics and are 
related to the Kostka numbers [79] . 
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5.3 Analytic Results for SU(Ar) Integrals 
5.3.1 Introduction 

Much work has been carried out on the topic of integration over the classical compact 
groups. The subject has been studied in great depth in the context of random matri- 
ces and combinatorics. Many analytic results in terms of determinants are available for 
integrals of various functions over unitary, orthogonal and symp lectio groups [80]. Unfor- 
tunately similar results for SU(A^) are not to our knowledge available. The primary use of 
these integrals has been in the study of Ulam's problem concerning the distribution of the 
length of the longest increasing subsequence in permutation groups [81, 82]. Connections 
between random permutations and Young tableaux [83] allow an interesting approach to 
combinatorial problems involving Young tableaux. A problem of particular interest is the 
counting of Young tableaux of bounded height [77] which is closely related to the problem 
of counting singlets in product representations mentioned in Section 5.2. Group integrals 
similar to those needed in this thesis have also appeared in studies of the distributions of 
the eigenvalues of random matrices [84, 85]. 

In the context of LGT not much has work been done in the last 20 years on the subject 
of group integration. The last significant development was due to Creutz who developed a 
diagrammatic technique for calculating specific S\J{N) integrals [86] using link variables. 
This technique allows strong coupling matrix elements to be calculated for SU(A^) [76] 
and has more recently been used in the loop formulation of quantum gravity where spin 
networks are of interest [87, 88, 89]. 

In sections 5.3.2, 5.3.4 and 5.3.3 we extend the results of Section 5.2 to calculate a 
number of SU{N) integrals. As generating functions these integrals allow the evaluation 
of all expectation values appearing in variational calculations of SU{N) glueball masses in 
2-1-1 dimensions. To calculate these generating functions we work with plaquette variables 
and make use of techniques which have become standard practice in the fields of random 
matrices and combinatorics. In Section 5.3.2 we derive a generating function which allows 
the calculation of integrals of the form 




(5.17) 
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The work in Section 5.3.3 generalises the generating function of Section 5.3.2 allowing the 
calculation of more complicated integrals of the form 



f dU 

JSV{N) 



The third integral, which we calculate in Section 5.3.4, while not directly relevant to the 
calculation of massgaps here, is of importance in the study of combinatorics. It allows the 
calculation of integrals of the form 



m 

dU ■ 



SU(W) 



Tr(C/') Tr(C/') . (5.19) 



For each integral considered the approach is the same and proceeds as follows. We 
start with a calculation of a U(Af) integral. For example in Section 5.3.2 we calculate 

Gu(^)(c,d) = / dUe'^^^+'^^\ (5.20) 

This is a generalisation of Gu(7v)(c, c), an integral first calculated by Kogut, Snow and 
Stone [90]. We then make use of a result of Brower, Rossi and Tan [91] to extend the 
\J{N) integral to SU(A^) by building the restriction, det ?7 = 1 for all U G SU(Af), into the 
integration measure. In this way S\J{N) generating functions can be obtained as sums of 
determinants whose entries are modified Bessel functions of the first kind. 



5.3.2 A Simple Integral 

In this section we introduce a useful technique for performing SU(A/") integrals. We start 
with the U(7V) integral of Eq. (5.20) and calculate it using a technique which has become 
standard practice in the study of random matrices and combinatorics. 



Since the Haar measure is left and right invariant (see Eq. (5.2)) we can diagonalise U 
inside the integral as 



/ pi<l>0 



u = v 



e^"^! 



\ 



V 



(5.21) 
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In terms of the set of variables, {(l)k}k=i^ the Haar measure factors as dU = dfi{(/))dV [90]. 
Since the integrand is independent of V, the V integral can be carried out trivially using 
the normalisation of the Haar measure given by Eq. (5.3). 



Making use of the Weyl parameterisation for U(iV) [92], 



N 



1=1 



(5.22) 



where A{(f)) is the Vandermonde determinant, with implicit sums over repeated indices 
understood. 



A(<^) 



1 



p. . . piM^-h)JM^-i2) . . . pi<t>N{N-i]^) 



(5.23) 



we can express the \J{N) generating function as follows 



Gu(Ar)(c, d) 



2n 



27r ■■7o 



d^ 



N 



27r 



exp 



N 



^{ce'^^+de 
.1=1 



|A(0)|2. (5.24) 



Substituting Eq. (5.23) gives. 



Gu(Ar)(c, d) 



N r2 



k=l 



H 1^ ^exp i(jfe-ife)</>fc + ce^'^'=+de-<^'= . (5.25) 



To simplify this further we need an expression for the integral, 

gn{c,d) = / — exp{inx + ce^^ + de *^), 
Jo 27r 

which is easily handled by expanding the integrand in Taylor series in c and d, 



9n{c, d) 



c'd'- f dx 

f^f^wy 2^ 



ix{k—l+n) 



= E 



fe=0 /=0 



°° k^k+r, 



k=0 
C 



A;!(A; + n)! 

/„ (^2Vcd) . 



(5.26) 



(5.27) 



5.3 Analytic Results for SU(Ar) Integrals 



71 



Making use of Eq. (5.27) in Eq. (5.25) gives an expression for Gu(Ar)(c, d) as a Toeplitz 
determinant^ , 

1 ^ 

Gl!(N){c,d) = -^£iii2...jjvS^Ju'2...jjv Jl ^^^-^^(C) c?) 

fc=l 



^ ^ fe=l 

(2\rcdS\ . (5.28) 

V ) \\<iA<N 



Here the quantities inside the determinant are to be interpreted as the (z,_7)-th entry of an 
NxN matrix. Now to calculate the corresponding SU(A/") result the restriction det U = 1, 
which is equivalent to X^j^i <t>k=^ mod 2tt in terms of the 0^ variables, must be built into 
the integration measure. To do this we follow Brower, Rossi and Tan [91] and incorporate 
the following delta function in the integrand of Eq. (5.24): 

(N \ oo / N \ 

^0fc-Omod27r| = ^ 2Tr6 i^cpk - ^Trm] . (5.29) 
k=l ) m=-oo \fc=l / 

This is most conveniently introduced into the integral via its Fourier transform, 

oo / ^ \ 

^ exp I 0fe j . (5.30) 



k=\ 



To obtain the SU(A/") integral from the corresponding U(A?^) result the modification is 
therefore trivial. Including Eq. (5.30) in the integrand of Eq. (5.24) leads to the general 

SU(Af) result, 



Gsu(JV)(c,d) = / dUe"" 
Jsv{N) 

oo 

= XI det[5„+j_i(c,ci)]^<.^^.<^. (5.31) 



cTr;7+dTr;7t 



This expression can be manipulated to factor the d/c dependence out of the determinant 



Toeplitz determinant is defined as a determinant of a matrix wliose (i, ji)-tli entry depends only on 

j -i- 
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as follows, 



det 



Il+j-i (2Vcd 



l<i,j<N 



m 



■■lll2--iN'^Jlj2---jN 



m=l 

- det 



l<i,j<N 



.(5.32) 



Making use of this result in Eq. (5.31) leads to the S\J{N) generating function, 



Gsv(N){c,d) = X! (~) Il+j-i [2Vcd 

l=-oo ^'^'^ 



i<i,j<N 



(5.33) 



For the case of SU(2) wc can show that this reduces to the standard result of Arisuc given 
by Eq. (5.7). To do this we need the recurrence relation for modified Bessel functions of 
the first kind, 



2n 

In-l{x) - In+l{x) = —In{x). 



(5.34) 



Recall that for SU(2) the Mandelstam constraint is TrU = TrU^ , so the case Gsu(2)(c, c) 
can be considered without loss of generality; 



G'su(2)(c,c) = J2 [Iii^cf - Ii-i{^c)Ii+^{2c)] 

l=—oo 

= /o(4c)-/2(4c). 



(5.35) 



Here wc have used the standard addition formula [93] for modified Bessel functions. Em- 
ploying the recurrence relation of Eq. (5.34) gives 



1 



Gsu(2) (c,c) = ^-fi(4c). 



(5.36) 



which is the standard result of Arisue given by Eq. (5.7). 



With an analytic form for S\J{N) in hand we can attempt to find simpler expressions 
for Gsu(N)(c, d) analogous to Eq. (5.36). To our knowledge no general formulas are avail- 
able for the simplification of the determinants appearing in Eq. (5.33). Without such 
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formulas we can resort to the crude method of analysing scries expansions and comparing 
them with known expansions of closed form expressions. Since the determinants appearing 
in the generating function are nothing more than products of modified Bessel functions we 
expect that if a closed form expression for the general SU(Ar) generating function exists, it 
will involve the generalised hypergeometric function. With this approach we have limited 
success. The SU(3) result of Eq. (5.14) is recovered numerically but the SU(4) result does 
simplify analytically. 



When analysing the series expansion of Gsu(4) (c, c) we notice that it takes the form of 
a generalised hypergeometric function. In particular we find the following result: 



GSU(4) (C, C) = 2F3 



3 5 
2' 2 

3,4,5 



;16c^ 



(5.37) 



Here the generalised hypergeometric function is defined by 



ai, 02, . . . ,ap 

61, 62, ■■ -,6(7 



00 

= 1: 



(ai)fc(«2)fc • • • (Qp)fc 



(5.38) 



where (x)^ = x{x + 1) . . . {x + k — 1) is the rising factorial or Pochhammer symbol. In 
addition to Eq. (5.37) we find the following results for matrix elements derived from the 
SU(4) generating function: 



{Zi) 



C2F3 



2F3 



5 I 

2' 2 

4,5,6 

3 5 
2' 2 

3,4,5 



;16c2 



16c2 



(5.39) 





5 7 
5,6,7 


+ icSF3 


7 9 

^'^ ;16c2 
_6,7,8 


+ ^cSF3 


9 11 

^' 2 ;16c2 
7,8,9 


2F3 


3 5 

^'^ ;16c2 

3,4,5 





.40) 
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2-^3 


3 5 

3,5,6 


+ |c2 2F3 


"57 

2' 2 ;16c2 

_4,6,7 


+ icSF3 


7 9 

_5,7,8 




3 5 

^'^ ;16c2 
3,4,5 





We stress that these results are nothing more than observations based on series expansions. 
Despite some effort analogous expressions for iV > 4 have not been found. 



The generating functions, Ggu(Ar) (c, d) and Gu(Af)(c, c?), are not only of interest in 
Hamiltonian LGT. By differentiating Eq. (5.33) appropriately with respect to c and d and 
afterward setting c and d to zero, we obtain the number of singlets in a given product 
representation of S\]{N). This was discussed for the special case of SU(3) in Section 5.2. 
We now consider the general case in the calculation of Tk{n); the number of singlets in 
the SU(A;) product representation, 

(k (g) k) (g) • • • (g) (k (g) k) . (5.42) 

n 

As a group integral Tk{n) is given by 

Tk{n) = [ dU{\TTU\^r. (5.43) 

.ysu(fe) 

Integrals of this kind arc studied in combinatorics, in particular the study of increasing 
subsequences of permutations. An increasing subsequence is a sequences ii < i2 < ■ ■ ■ < im 
such that 7r(ii) < 7r(z2) < ■ ■ • < 7r(z^), where tt is a permutation of {1,2,..., k}. It has 
been shown that the number of permutations tt of {1, 2, . . . , A;} such that tt has no increas- 
ing subsequence of length greater than n is Tfe(n) [81]. In addition it is possible to prove 
that Tk (n) is the number of pairs of Young tableaux of size k and maximum height n via 
the Schensted correspondence [94, 81]. 



Making use of Eq. (5.33) we see that only the I = term contributes to Tk{n). Letting 
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X = cd we have: 



/ q2 \n r 



l<i,j<k 



c=d=0 



Hence the generating function for Tfe(n) is given by 



l<i,j<k 



(5.44) 



x=0 



E 

n=0 



n! 



|2 



det [ij_i(2v^)] 



l<jJ<A; ■ 



(5.45) 



a result first deduced by Gessel [77]. The first few T};{n) sequences are available as A072131, 
A072132, A072133 and A072167 in Sloane's on-line encyclopedia of integer sequences [95]. 

5.3.3 A More Complicated Integral 

We now move on to the more complicated integral 



i{c,d) = / dU e^p 
VSU(iV) 



c{TiU + TtU^) + dTiiU"') Vm G Z+. (5.46) 
This integral is of interest as a generating function for the calculation of integrals such as 



[ dUTk{U'^)e<^^+'^^^\ 



(5.47) 



For the simple case of SU(3) we can use the Mandelstam constraint, Tr([/^) = (TrC/)^ — 
2TrL'"''', to reduce such integrals to those obtainable from the generating function of the 
previous section. However for higher dimensional gauge groups not all trace variables can 
be written in terms of TrU and TrU^. For these gauge groups one must introduce the 
generating function, Hjn{c,d), to calculate integrals similar to Eq. (5.47). 



To calculate H„i{c,d) we start with the corresponding U(A'^) generating function and 
follow the procedure of Section 5.3.2 to obtain 



hmic,d) = [ dUexp\c{TiU + TiU^) + dTi{U"') 
Jv{N) L 



■^Sii-iN^ji-jNYiJ^ ^exp i(jfc-ifc)0fc + 2ccos^A; + rfe™'^'= .(5.48) 
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To proceed we need the following integral, 

/•27r J °o tfc /•27r t 

/ ^exp(mx + acosx + 6e*™^) = Y Tl ^e^('^+™'=)^+'^"°^^ 
^ 27r fn^^'io 27r 



= ^^/„+mfc(a). (5.49) 

A;=0 

Making use of Eq. (5.49) in Eq. (5.48) leads to the following expression for the ^{N) 
generating function, 

hm{c,d) = det[\m;j-i{c,d)\-^^-.^j^, (5.50) 

with 

Am;n(c, d) = ^ —In+mk{^c). (5.51) 
fe=0 

Extending to SU(A'^) following the prescription of Section 5.3.2 we arrive at the corre- 
sponding SU(A^) generating function. 



Hm{c,d) = ^ dei[\ra;i+j_i{c,d\^. .^j^. (5.52) 

l=—oo 

An example of an SU(Ar) integral derived from this generating function is the following: 

Jsv(N) <jd 



'SU(JV) 

" dd 



= — CX) 



d=0 

(5.53) 

d=0 



d °° 

— det[li+j-i{2c)+dli+j-i+m{^c)] 



Only two terms need to be kept in the fe-sum of Eq. (5.51) here because higher order 
powers of d vanish when the derivative with respect to d is taken and d set to zero. 



5.3.4 Another Integral 

Another integral that is calculable using the procedure of Section 5.3.2 is 

JSV{N) 
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Wc start with the corresponding U (N) integral and again use Weyl's parameterisation for 
U(iV) to give 



,(c,d) = / dUe 

JV{N) 



cTrC/^+dTrC/- 



(5.55) 

i{jk - ik)cpk + £6"^'^*= + de-^'^A . (5.56) 



To simplify this further we need the following integral: 



5m;n(c,d) = 



27r 



dx 



^ k\{k + n/m)\ 



k=0 



giLicjd) ifmln, 

m 

otherwise. 



k=0 1=0 

if m\n, 
otherwise. 



c^d^ dx 



ix[m(k—l)+n\ 



k\l\ Jo 27r 



(5.57) 



Making use of Eq. (5.57) in Eq. (5.56) leads to 



jmic,d) = det[gm;j-iic,d)]^^..^j^. 



(5.58) 



The corresponding SU{N) result is then 



Jm{c, d) = ^ det [gm;i+j-i{c, d)] 



l<i,j<N ■ 



l=oo 



(5.59) 



Due to the form of gm;n{c,d) the U(A'') generating function simplifies significantly for 
m> N. For this case we have, 



jm{c, d) — ^j£n...jjv^jl---iAr5m;ji-ii (c, d) ■ ■ ■ gm;jiy,-ijy,{c, d) 

= 9o{c,d)^, 



(5.60) 



since the only non-zero terms in the determinant occur when ji — ii = for I = 1, . . . ,N. 
We can make use of this result in the calculation of 



Am) 
JnN 



U(iV) 



dU |TV(?7"*) 



m\n\2 



(5.61) 
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a quantity which has been studied in the context of combinatorics [81] and random ma- 
trices [84] . Making use of our generating function we see that 



(5.62) 



c=d=0 



For the special case of m > iV we can make use of Eq. (5.60) to express as the sum 
of squares of multinomial coefficients as follows: 



c=d-- 



.0 V^^^c/; ) 



N 



c=d=0 



{cdY 



\dcdd ^ ^ (ii!)2---(/Ar!)2 

ilH l-'iv=fc 



c=d=0 



E 



|2 



ll,...,lf^>0 
'iH ^-ll\^='n 



(^i!)2---(/iv!)2- 



(5.63) 



With the help of Mathematica we can express with m > N in terms of sums of 
binomial coefficients for N = 2,3, 4: 



Am) 
Jn,2 



Am) 
Jn,3 



Am) 
JnA 



n\ In 



E 

fc=0 
n 

E 

fc=0 
n 

E 

fe=0 



n 



n\ ^ /'2k 



k \k 



n\V2n - 2k\ f2k 



k 



k J \ n — k 



(5.64) 



We have not found similar expressions for N > 5. 



A further exercise is to explore the large N limit of the generating function, jm{c,d). 
Let us expand jm{c, d), in the form of Eq. (5.55), in c and d about c = d = 0: 



Jr 



k% Mn) 



Tr{U"') Ti{U 



(5.65) 



A result due to Diaconis and Shahshahani [84], 



/ dU 

Jv(N) 



Tr(f/'") Tr(f/"») = dkiklm'' ViV > m/s. 



(5.66) 
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then allows us to access the large limit. In the ^ oo limit, for each term of the sum 
in Eq. (5.65) we have > mk and hence, using Eq. (5.66), we have. 

This result is consistent with the c = d and m = 1 case proved in Ref . [80] , 

lim det [V,(2c)],< = (5.68) 

iV — ^oo — — 

which is valid for all real c > 0. 



For the case of = mom, where mo € N, we also find an interesting simplification for 
jm(c, d). In this case the determinant appearing in Eq. (5.58) becomes a block determinant, 
which leads to the simplification. 



det[gm;j-i{c, d)] 



l<i,j,<mom 



det [gj-i{c, d)] 



l<i,j<mo 



(5.69) 



This result is equivalent to jm{c,d) 



N=mom 



= jiic,d) 



N=mo 



5.4 Conclusion 

In this chapter wc have calculated a number of SU(A^) group integrals. As generating 
functions they allow the analytic calculation of all expectation values that appear in vari- 
ational calculations of SU(A'') glueball masses in 2+1 dimensional Hamiltonian LGT. Our 
derivations have made use of techniques for the calculation of U(A/") integrals which have 
become standard tools in the fields of random matrices and combinatorics. The resulting 
generating functions generalise the analytic results that have been available for SU(2) for 
many years to the general case of SU(A^). 



The extension of the techniques presented here to 3+1 dimensions is highly non-trivial. 
The starting point for all derivations presented in this chapter was the fact that in 2+1 
dimensions the transformation from link to plaquette variables has unit Jacobian. This 
allows the matrix elements of functions of plaquette variables to be calculated analytically 
by differentiating a generating function appropriately. The viability of applying analytic 
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techniques in 3+1 dimensions is discussed in Chapter 8. 

Having derived analytic results for a selection of SU(iV) generating functions, in the 
next chapter we apply them in a variational calculation of SU(iV) glueball masses, with 
AT = 2, 3, 4 and 5, in 2+1 dimensions. The calculations make use of the simplest classically 
improved and tadpole improved Hamiltonians derived in Chapter 3 as well as the standard 
Kogut-Susskind Hamiltonian. 



Chapter 6 



SU(iV) Massgaps in 2+1 
Dimensions 



6.1 Introduction 



In this chapter we calculate the lowest lying glueball masses, or massgaps, for SU(A^) 
LGT in 2+1 dimensions, with = 2, 3, 4 and 5. We use the Kogut-Susskind, classically 
improved and tadpole improved Hamiltonians in their calculation and rely heavily upon 
the results of the previous chapter for the calculation of the required expectation values. 

The outline of this chapter is as follows. After defining our notation in Section 6.2 
we fix the variational vacuum wave function in Section 6.3. Following that we calculate 
lattice specific heats in Section 6.4 before studying SU(iV) glueball masses in Section 6.5. 
In Section 6.6 we present our results and compare them to similar calculations in both 
the Lagrangian and Hamiltonian formulations. Section 6.7 contains a discussion of further 
work and our conclusions. 
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6.2 Preliminaries 

Before fixing the variational ground state we introduce the following convenient notation. 
We define the general order a} improved lattice Hamiltonian for pure SU(-/V) gauge theory, 



7i:(k,«o) = ^I^Tr 



(1 - K)Ei{xf + -^Ei{x)\Ji{x)Ei{x + a%)V\{x) 



+^ E {{l + ^i^)Pij{x)-'^[Rij{x) + Rji{x)]], (6.1) 



x,i<j 



where the plaquette and rectangle operators are given in Eq. (3.64). The simplest Hamilto- 
nians derived in Chapter 3 can be expressed in terms of H as follows. The Kogut-Susskind 
(Eq. (3.37)) and C(a^) classically improved (Eq. (3.60)) Hamiltonians are given by H{0, 1) 
and ii,{l/6,l) respectively. The tadpole improved Hamiltonian (Eq. (3.63)) is given by 
H{l/6,uo), where uq is defined self-consistently as a function of /? = N/g^ as described in 
Section 4.3.4. 



With this new notation the vacuum energy density is given by 



Nr 



■in) 



(1 




+ 2(1 + 3k)/3, 



(6.2) 



where the expectation values, as usual, are taken with respect to the one-plaquette trial 
state defined in Eq. (4.22). The variational parameter, c, is fixed as a function of (3 by 
minimising the vacuum energy density. For the calculation of the expectation values we 
use the generating functions of Chapter 5 with all infinite sums truncated. Once the 
variational parameter is fixed the trial state is completely defined as a function of /?. 



6.3 Fixing the Variational Trial State 
6.3.1 Introduction 

In this section we fix the SU(iV) trial state for 2 < iV < 5, making use of the generating 
functions derived in Chapter 5. These generating functions allow the analytic calculation 
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of the plaqucttc and rectangle expectation values appearing in Eq. (6.2). The approach 
we take follows Section 4.3.2. For the Kogut-Susskind and classically improved cases, we 
simply minimise eo for a given value of (3. The value of c at which eo takes its minimum 
defines c as a function of (3. The tadpole improved calculation is more complicated because 
the mean plaquette, from which uq is calculated, depends on the variational state, |^o)- 
To complicate matters, the variational state in turn depends on the energy density and 
hence uq. For this case we make use of the iterative technique procedure introduced in 
Section 4.3.4. 

6.3.2 Results 

The results of the Kogut-Susskind, classically improved and tadpole improved SU(2), 
SU(3), SU(4) and SU(5) vacuum energy density calculations are shown in Fig. 6.1. The 

corresponding variational parameters, c(/3), are shown in Fig. 6.2. For SU(3) the gen- 
erating function of Eq. (5.14) is used to calculate the required plaquette and rectangle 
expectation values. The generating function of Eq. (5.33) is used for SU(4) and SU(5). 

The familiar strong and weak coupling behavior from variational calculations is ob- 
served in each case. The differing gradients for the improved and Kogut-Susskind varia- 
tional parameters for a given N in the weak coupling limit highlight the fact that when 
using an improved Hamiltonian one is using a different renormalisation scheme to the 
unimproved case. 

6.3.3 Dependence on Truncation 

In practice the A;-sum appearing in the SU(3) generating function (Eq. (5.14)) is trun- 
cated at a maximum value, kmax- The dependence of the variational parameter on various 
truncations of the fc-sum is shown in Fig. 6.3. We see that convergence is achieved up to 
/3 ~ 13 when keeping 20 terms. Further calculations show that when keeping 50 terms 
convergence up to /? 30 is achieved. 

The Z-sum appearing in the general SU(Ar) generating function of Eq. (5.33) is also 
truncated in practice. We replace the infinite sum over Z by a sum from — Zmax to Zmax- The 
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Figure 6.1: Analytic calculation of the 2+1 dimensional unimproved, improved and tadpole 
improved vacuum energy density in units of l/(aNp) for SU(2), SU(3), SU(4) and SU(5). 



dependence of the SU(3) and SU(4) variational parameters on /max is shown in Fig. 6.4. 
From the graphs we see that convergence is achieved quickly 8.8 /max incrGciscs foi both 
SU(3) and SU(4). The results for I 

max ^ 8 Siie barely distinguishable up to /3 — 80 with 

the scale used in the plots. 



6.4 Lattice Specific Heat 

In addition to the vacuum energy density we can also calculate the lattice specific heat 
defined in Eq. (4.20). The results for SU(2), SU(3), SU(4) and SU(5) are shown in Fig. 6.5. 
The SU(2) and SU(3) results are calculated with the aid of Eqs. (5.7) and (5.14) with the 
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Figure 6.2: Analytic calculation of the unimproved, improved and tadpole improved vari- 
ational parameter in 2+1 dimensions for SU(2), SU(3), SU(4) and SU(5). 



fc-sum of Eq. (5.14) truncated at kmax = 50. The SU(4) and SU(5) results are obtained 
using Eq. (5.33) to calculate the required matrix elements. For these cases the infinite 
Z-sum is truncated at ^max = 4, for which the generating function has converged on the 
range of couplings used. We recall from Chapter 4 that the location of the peak indicates 
the region of transition from strong to weak coupling. For an improved calculation one 
would expect the peak to be located at a smaller (3 (corresponding to a larger coupling) 
than for the equivalent unimproved calculation. We see that this is indeed the case for 
each example, with the tadpole improved Hamiltonian demonstrating the largest degree 
of improvement. 
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5 10 15 20 

Figure 6.3: Analytic calculation of the unimproved SU(3) variational parameter in 2+1 
dimensions, truncating the /c-sum of y{c,d) at k = kmax- 

6.5 Massgaps 
6.5.1 Introduction 

Having fixed the one-plaquette vacuum wave function, in this section we turn to investi- 
gating excited states. Our aim is to calculate the lowest lying energy eigenstates of the 
Hamiltonians described by Eq. (6.1) for SU(A'') with 2 < N <5. 

We follow Arisue [30] and expand the excited state in the basis consisting of all 
nxm rectangular Wilson loops {|n, = {j/)}^^'"' that fit in a given square whose 
side length Lmax defines the order of the calculation. Enumerating the possible overlaps 
between rectangular loops is relatively simple and so a basis consisting of rectangular loops 
is an ideal starting point. However, for an accurate picture of the glueball spectrum we will 
need to extend the rectangular basis to include additional smaller area loops. Without such 
small area nonrectangular loops, it is possible that some of the lowest mass states will not 
appear in the variational calculation described here. In order to ensure the orthogonality 
of l^o) and |0i) we parameterise the excited state as follows 

T 7"2 

^max max 

|0i) = XI Sn,m\n,m) = si\l), (6.3) 

n,m=l 1=1 
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Figure 6.4: Analytic calculation of the unimproved SU(3) and SU(4) variational parameters 
in 2+1 dimensions, truncating the /-sum of Eq. (5.33) at I = ilmax- 



with, 



10 



Y,mix)-{Wi{x))]\cpo). 



(6.4) 



Here {Wi{x)) is the expectation value of Wi{x) with respect to the ground state \(po) and 
the convenient label I = {n — l)Ljnax + has been defined to label the nx m rectangular 
state, \n, m). We define the particular form of Wi{x) to reflect the symmetry of the sector 
we wish to consider. For SU(A/") we take Wi{x) = Tr[wi{x) ± wj{x)] for the symmetric 
(0''"''") and antisymmetric (0 ) sectors. To avoid over-decorated equations, the particular 
Wi{x) in use is to be deduced from the context. Here wi{x) is the rectangular Wilson loop 
joining the lattice sites x, x + nai, x + nai + maj and x + naj, with 



n 



I - 1 



+ 1 and m = I — Lr, 



l-l 



(6.5) 



Here [k] denotes the integer part of k. In order to calculate excited state energies we 
minimise the massgap (the difference between the excited state and ground state energies) 
over the basis defined by a particular order L^ax- To do this we again follow Arisue [30] 
and define the matrices 



Niv 



1 



— {l\H-E^\l'), 



(6.6) 
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where Eq is the ground state energy, and 



Dw = ^m = Y.[{w!{x)w,m-{wi{x)r{wi,m]. (6.?) 



Extending the calculation to the general improved Hamiltonian H{k, uq) and making use 
of Eqs. (A. 15 — A.18) from Appendix A gives 



■f^ E E i (1 - ( [^f m^), w.m 



l,X X' 



+- 



(6i 



To minimise the massgap over a basis of a given order we make use of following diagonal- 
isation technique [96]. We first diagonalise Z), with 



(6.9) 



where V is diagonal. The n-th lowest eigenvalue of the modified Hamiltonian, 



H' = VS^NSV, 



(6.10) 



then gives an estimate for the massgap corresponding to the n-th lowest eigenvalue of the 
Hamiltonian, Am„. 



6.5.2 Classification of States 

States constructed from only gluon degrees of freedom can be classified in the continuum 
by their J^'^ quantum numbers. In Section 2.5 we discussed the assignment of P and C 
quantum numbers. In this section we discuss the important topic of building states with 
particular continuum spins on the lattice. Difficulties arise when the continuous rotation 
group of the continuum is broken down to the group of lattice rotations. The most serious 
difficulty to arise is an ambiguity in the assignment of continuum spins to states built from 
lattice operators. It is this ambiguity that we discuss in what follows. 

Let us start with an eigenstate, \(f)), of the Hamiltonian. In the continuum, follow- 



6.5 Massgaps 



89 



ing standard quantum mechanics, one can construct a state with spin J by taking the 
superposition 

f-2-K 

\ct>j) = / e'^'nel^), (6.11) 
Jo 

where TZo is the operator which rotates by an angle 9. This fohows from the fact that J is 
the generator of rotations. On a square lattice continuous rotations are not available. All 
rotations must be modulo 7r/2 and so the lattice analogue of Eq. (6.11) is 

3 

\<Pj) =J2e'^TZn^\<P). (6.12) 
n=0 

Here the state labelled {(pj) is no longer a pure spin J state. It also includes states with 
spins J± 4, J lb 8, . . . since these spins produce the same phases as J on the square lattice 
as shown in Eq. (6.12). To clarify this point continuum states with spin J it 4, J it 8, . . . 
will couple to the lattice state labelled J = 0. Thus we need to be careful when using 
continuum terminology to label lattice wave functions. 

To give a specific example, suppose we construct a wave function, \(p), on the lattice 
with lattice spin J = 0; a state built from Wilson loops which are unchanged by rotations 
of n7r/2 for all integers n. As explained above, this state is not a pure J = state; it also 
contains J = 4, 8, . . . states. Using a variational approach we can obtain estimates of the 
lowest energy eigenvalues of states with lattice spin J = 0. When the continuum limit is 
taken we obtain estimates of the lowest continuum energy eigenvalues for the states with 
spin 0, 4, 8, . . .. 

In the Lagrangian approach it is possible to suppress the unwanted spin Jib 4, Jib 8, . . . 
states in a given spin J calculation. By "smearing" links, one can confidently construct 
lattice states which do not couple with the unwanted higher spin continuum states, at 
least for the lowest energy eigenvalues [51]. The technique of "smearing" links has not to 
our knowledge been applied in the Hamiltonian approach. 

Another way to clarify ambiguities in spin assignment is to attempt to construct a 
J = 4 state on the lattice, devoid of J = contributions, and similarly a J = state. 
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devoid of unwanted J = 4 contributions. The construction of an exact J = 4 state is 
impossible on a square lattice due to the unavailability of 7r/4 rotations. One can however 
attempt to construct states that are approximately symmetric under rotations by 7r/4. In 
the Lagrangian approach in 2+1 dimensions, it has been demonstrated, for the simple case 
of SU(2), that such states can be chosen on a square lattice and that the approximate 7r/4 
symmetry becomes exact in the continuum limit [52] . This technique is readily applicable 
in the Hamiltonian approach but has not yet been attempted. 

Thus in our Hamiltonian calculation a lattice spin J state will correspond to a con- 
tinuum state with spin J, J ± 4, J ± 8, . . .. Using a variational approach we can obtain 
estimates of the lowest mass states in the continuum with these spin values. To improve 
the spin identification in the continuum beyond modulo 4 requires more work. Perhaps 
the multiple scaling regions visible in our massgap calculations highlight this ambiguity. 
Either way, it will prove interesting to compare the masses calculated here to that of Teper 
who has been careful to identify continuum spins correctly, at least for the lowest mass 
excitations. 

6.5.3 Calculating Matrix Elements 

Having described the minimisation process we now detail the calculation of the matrix 

elements Nut and Dul Our aim is to reduce Nui and Dut to polynomials of one plaquette 
matrix elements. This, again has been done for the case of SU(2) by Arisue [30]. Here 
we retrace his calculations for the general case of SU(A'') and extend them to incorporate 
improved Hamiltonians. We start with Dul 

Taking elementary plaquettes as our independent variables, it is easy to show that 
the only non-zero contributions to Dui occur when the rectangles I and V overlap. As 
an example of a contribution to D/;/, consider ADur; the case where iV/n;/ plaquettes are 
contained by both rectangles (these are the overlap plaquettes) and Ni plaquettes are 
contained by the rectangle I. In order to calculate such matrix elements we rely heavily on 
the orthogonality properties of SU{N) characters. We are interested in calculating SU(A'') 
integrals of the form 




(6.13) 
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where Up is a SU(A^) plaquctte variable and y is a product of any number of plaquettes 
not including Up. Here Xr{U) denotes the character corresponding to the representation 
r. For SU(2), r = 0, 1/2,3/2, ... and for SU(3), r = (A,/x) where A denotes the number 
of boxes in the first row of the Young tableau describing the representation and /x is the 
number of boxes in the second row. Similarly, for SJJ{N), r = (ri,r2, . . . ,rN-i). 



Performing a character expansion of the exponential in Eq. (6.13) gives: 

j dUpe^^^^hriUpV) = Y.j dUpCr'Xr'{Up)Xr{UpV). (6.14) 

r' 

This is simply a generalisation of a Fourier expansion. Here, the coefficient c^' is given by: 

Cr> = J dUpxAUp)e^^^-\ (6.15) 
Now, using the orthogonality relation, 

j dUpXAUpV)Xr{Up) = j-Sr'rXriV), (6.16) 

where dr is the dimension of the representation r, we obtain: 

J dUpe^^^-^XriUpV) = ^Xr{V) J dUpXr{Up)e^^^'^\ (6.17) 

This result allows us to integrate out a single plaquette from an extended Wilson loop in 

a given representation r. To complete the calculation we need to relate SU(A^) characters 
to traces of group elements. This can be done using Weyl's character formula [97]. For 
SU(iV), according to Bars [97], the dimensions and characters corresponding to the first 
few representations are given by: 



Xi{U) 


= TrC/ 


di{U) 


= N 


X2{U) 


= ^ [{TvUf + Tr[/2] 


d2{U) 


= In{n + i) 


XniU) 


= ^ [{TiUf - Tr?/^] 


du{U) 




X2l{U) 


= \ [{TiUf - TiU^] 


d2i{U) 


= ^{N -1)N{N + 1 


Xi^-^{U) 


= Tr?7^ 


diN-i{U) 


= N 



(6.18) 
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Here we have adopted the convention of dropping all zeros in the character labels. The 
Mandelstam constraints for the gauge group in question allows all characters to be ex- 
pressed in terms of a minimal set of trace variables. For example, for SU(3) we make use 
of the Mandelstam constraint, 

TriA^B) = Tr:ATr{AB) - TrA^^TrB + Tr{A^B), (6.19) 

where A G SU(3) and B is any 3x3 matrix, to express all characters in terms of TrU and 
TrU^. For example, for the case of SU(3), Eq. (6.18) simplifies to 



Xi{U) = TrU di = 3 

(6.20) 



Xn{U) = J [(TVC/)2 - IV(C/2)] = TVC/t = 3 



X2iU) = ^ [(TrC/)2 + TV([/2)] = (W)2 - TVC/t d2 = 6 

X2i{U) = ^[{TrUf-Tr{U^)]=TrUTrU^-l dsi = 8 

However, for general SU(A^) the Mandelstam constraints are difficult to calculate. In what 
follows we will need expressions for TrU, TrU^^, Tr{U^), {TrUf, and TrUTrU^i as Unear 
combinations of characters for SU{N). It is possible to do this without the use of the 
Mandelstam constraint. Such expressions are necessary in order to make use of Eq. (6.17) 
in the calculation of expectation values of trace variables. For TrU, TrU^ , Tr{U'^) and 
(TrU)^ the necessary expressions are easily obtained by rearranging Eq. (6.18), 

TrU = xi{U) 

{TrUf = X2{U)+Xn{U) 

TrU^ = X2{U)-xn{U) 

T>C/t = Xi^-^{U). (6.21) 



To express the remaining expression, TrUTrU^, linear combination of characters is 
not as easily done. For SU(3) one can simply rearrange the expression for X2i{U) in 
Eq. (6.20). For the general A'" case it is simplest to consider Young tableaux. In terms of 
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Young tableaux we have 



□ 



□ 



(6.22) 



Performing the product representation decomposition gives 



□ 



□ 
□ 

□ 



\ N-1 = 



□ 
□ 

□ 



□ 

> iV+ : 

□ 



} N-1 



(6.23) 



Converting back into the notation of characters and traces gives 



TvUTiU^ = l + X2iJv-2(?7). 



(6.24) 



Eq. (6.24) together with Eqs. (6.21) and (6.17) allow the analytic calculation of each 
contribution to Dui for all N. For the case of ADu' described earlier, we have 



ADii> = -Fz,{Ni + Ni, - 2Nini' 



{Nir^v) + Fz^z^ {Nmi') - 4Fz, {Ni)Fz, (iV;0,(6.25) 



where the character integrals are given by: 

n-l 



FzAn) 



(6.26) 



1 



1 



2 [A^(A^ + 1) 



n— 1 



{zl + + 



1 



2 IN{N - 1) 



n-l 



{Zl-Z^r, (6.27) 



2 [Af(iV + 1) 



n-l 



{zi+z^r 



1 



N{N - 1) 



n-l 



Fz,zM) = 1 + 



_(iV- 1)(7V + 1) 



n-l 



{Zl-Z^r, (6.28) 



(6.29) 



Here we have made use of the notation, Z„ := Tr (C/"), to denote the trace variables oc- 
cupying a single plaquette, U. The expectation values appearing in Eqs. (6.26 — 6.29) are 
easily expressed in terms of the generating functions of Eqs. (5.33) and (5.52). Differenti- 
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ating the generating functions appropriately gives 

{Zi) = 

{Zf±Z2) = 



1 






dc 


1 


9'^Gsv{N) 




dcdd 


1 


9'^Gsu(^N) 




dc? 



d=c 



d=c 



1 dH2 



d=c 



H2 dd 



(6.30) 



d=0 



In practice, we do not need to calculate all of these matrix elements. We see from Eq. (5.31) 
that {ZD and {ZiZi) are related by 



{ZD 



1 d^G'su(jv)(c,c) 
2Gsu(Ar) (c,c) dc^ 



{ZiZi). 



(6.31) 



This follows from the fact that a group integral does not depend on the choice of direction 
for the links. To be more precise, the result 



/ dUf{U) = [ dU^fiU^) = [ dUfiU^), 

JSV{N) JSV{N) JSIJ{N) 



(6.32) 



follows from the fact that dU and dU^ each define invariant Haar measures on SU(A'') 
which, by uniqueness, must be equal. 



We now move on to the calculation of Nur. It is easy to show that the only non-zero 
contributions occur when there is at least one common link and an overlap between the 
rectangles. The improvement term (the second term in Eq. (6.8)) only contributes when 
the two rectangles share at least two neighbouring links in a given direction. Consider the 

contribution ANni to Nui in which there arc Li common links and L2 common strings of 
two links in a given direction. Again we suppose Ni plaqucttcs arc enclosed by rectangle I 
and that there are A^^nr common plaquettes. Making use of Eq. (6.17) and Eqs. (A. 15 — 
A. 18) from Appendix A we obtain 



ANw = ^Fz,iNi+Ni,-2Ninv) 



,(6.33) 
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with 



L = (1 - k)Li + —L2. 



(6.34) 



'0 



For the case of SU(3) we can simpUfy this using Eq. (6.19) to 



mi 



„ = ^Fz,{Ni+Ni,-2Ni^l,) lFz2{Ni^i,)-2Fz,{Ni^i>)-3+^Fz^zA^ini') ■ (6.35) 



Having determined individual contributions to Du' and Nui, to complete their calcu- 
lation the possible overlaps between states I and I' of a given type must be counted. 

6.5.4 Choosing an Appropriate Vacuum State 

In Section 6.3 we calculated variational vacuum wave functions for pure SU(A'^) gauge 
theory for N = 2, 3, A and 5. Our motivation was to use these wave functions as inputs to 
calculations of SU{N) massgaps. We obtained wave functions with a variational parameter 
that was proportional to P in the large (3 limit and in the small /? limit. However, this 
behaviour is incompatible with the exact continuum vacuum wave function. For the one- 
plaquette trial state, given by Eq. (4.22), to be compatible with the exact continuum 
vacuum wave function in 2+1 dimensions, one must have c oc in the large (3 limit [98]. 
This result is independent of the dimension of the gauge group in question. To understand 
this result we must return to continuum non-Abelian gauge theory in 2+1 dimensions for 
which Greensite [99] argued that the infrared behaviour of the vacuum can be described 
by the wavefunction, 



The simplest lattice analogue of Eq. (6.36) is the one-plaquette trial state given by Eq. (4.22) 
Expanding the exponent of Eq. (4.22) about a using Eq. (3.34), we see that for the 
Eq. (4.22) to reduce to Eq. (6.36) in the continuum limit and for n to remain finite, 
the variational parameter c must be proportional to 




(6.36) 



For the case of SU(2) in 2+1 dimensions, in the scaling region it has been shown [98] 
that for compatibility with the exact SU(2) continuum vacuum wave function we must use 
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the Greensite vacuum wave function, 

/ 0.81 ± 0.02 ^ \ „„x 

|<^o) = exp E LJ ^^"^^^ 

It is this vacuum wave function that was used by Arisue [30] in the calculation that we 
generalise in this chapter. It would thus seem that using a variational wave function is 
not appropriate in the calculation of massgaps. However, this is not entirely the case as 
we will now show. 



Define the following quantities. 




and 



= |<^i (c(/3),di (/?),..., d,(/3))) = F(di (/?),..., d,(/?))|(/)o(c(/?))). (6.38) 

Here E is an energy functional, H a Hamiltonian, a trial vacuum state, defined in 
Eq. (4.22), dependent on the coupling only though the dependence of c on /?. The dk{0) 
are determined variationally by minimising — -E[0o] for each (3 and the function F 

is chosen so that \(f)i) and |^o) are orthogonal. 



Suppose that we calculate E[<j)i] — E[<j)Q\. We see from Eqs. (6.7) and (6.8) that this 
quantity depends on the coupling only through c, di, . . . ,dn and a multiplicative factor of 

aAm{c{l3), di(/3), . . . , d„(/3)) = E[^i] - £;[0o] (6.39) 

Suppose that this calculation exhibits the appropriate scaling behavior on some f3 interval 
S: 

^Am = /x(c(/3),di (/?),..., d„(/3)) = Amo V/3 e 5. (6.40) 

Here Amo is a constant. The variational parameters may be rescaled with respect to /3 
and the same constant result is obtained but for a different range of couplings T: 



fi{c{fm,di{fm,...,dn{fm) = ^rno yper. (6.41) 
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Here T is the interval for which b £ T ^ f(b) £ S (i.e. T is the image of S under 
f~^). If we choose such that c{f{f5)) gives the correct scaUng form (for example 

/(/3) K 1.62/?2/8 for SU(2) in 2+1 dimensions) we then have 

aAm(c(/(/3)),di(/(/3)), . . . ,d„(/(/3))) ^ E[4>^] - Eo V/? G T, (6.42) 

where Eq is the exact ground state energy. This follows from the fact that |0o(c(/(/?)))) 
has the correct scaling form and so in the continuum limit we have ii"|0o(c(/(/?)))) = 
(^0 + £)l^^o(c(/(/3)))). Here e is a small correction term that would become smaller if we 
improved our trial vacuum by adding more complicated Wilson loops in the exponent. 

Since |(/>i) is constructed orthogonal to the ground state we have 

aAm{c{f{(3)),di{f{P)),...,dn{fm)>Ei-Eo + e V/3 G T, (6.43) 

where Ei is the exact energy of the first excited state, and so 

g'^Amo >Ei-Eo + e V/3 e T. (6.44) 

We cannot therefore say that we have an upper bound on the massgap, only an approxi- 
mation. 

This argument demonstrates that if one obtains scaling with a given c(/9) = ci(/5) in 
the vacuum wave function, the same result will be obtained using a different function, 
c{(3) = C2(/3), although in a different region of coupling, provided that there exists a func- 
tion / such that Ci(/3) = C2(/(/?)). For this reason we expect that the calculation of 
massgaps using the variational trial state will agree with a calculation using the Greensite 
vacuum wave function. 

The argument above also demonstrates that the precise form of c{f3) is unimportant in 
the calculation of massgaps in this context. A numerical demonstration of this for the case 
of 2+1 dimensional SU(3) gauge theory is shown in Fig. 6.6. We see that each glueball 
mass calculation approaches the same constant scaling value (labelled by "Continuum 
limit") but at different values of /3. As it turns out, to calculate variational wave functions 
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for large dimension gauge groups is cumbersome. Numerical precision becomes a factor 
in the minimisation of the energy density. This problem is magnified in the calculation 
of tadpole improved results. We thus abandon the use of variational wave functions in 
the calculation of massgaps beyond SU(5). Instead we make use of the one-plaquette 
wave function of Eq. (4.22) and define a simple dependence c{(3), which most often will be 
c{(3) = (3. Calculations for iV > 5 are presented in Chapter 7. 

6.6 SU(2), SU(3), SU(4) and SU(5) Massgap Results 

In this section we present glueball mass results for SU(Ar) pure gauge theory in 2+1 di- 
mensions with iV = 2, 3, 4 and 5. For each SU(3) calculation we keep 80 terms in the 
fc-sum of Eq. (5.14) giving convergence up to /3 = 50. For N > 3 the truncation Zmax = 20 
is used. The generation of A^^^' and Dw and implementation of the minimisation process 
is accomplished with a Mathematica code. 

For the case of 2+1 dimensions we expect aAm/g^ to become constant in the scaling 
region. The convergence of the massgaps with Lj^ax is illustrated in Fig. 6.7. We notice 
that for AT > 2 only small improvements to the scaling behaviour are gained by extending 
the calculation beyond order 8 on the range of couplings shown. This suggests that a more 
complicated basis (including, for example, nonrectangular loops) is required to simulate 
SU(A'') excited states with N > 2 than for the case of SU(2). 

In Fig. 6.8 results for the lowest lying glueball mass, calculated with Kogut-Susskind, 
improved and tadpole Hamiltonians, are shown. We see that aAmf/g^ is approximated 
well by a constant, in very large scaling regions, for the lowest lying eigenstates for all N 
considered. The scaling behavior becomes significantly worse for the antisymmetric sector 
which is shown in Fig. 6.9 and for higher energy eigenvalues. This is because the simplistic 
form of our excited state wave function is not sufficient to reproduce the plaquette corre- 
lations required to simulate these higher order states. One would expect the simulation of 
higher order eigenstates to improve by including more complicated loops in our expansion 
basis or by using a more complicated ground state. The continuum limit excited states 
results for SU{N) are given in Tables 6.1 — 6.7. 
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For the unimproved SU(2) case, the masses of the lowest two cigcnstates agree closely 
with the calculations of Arisue [30] (respectively 2.056 ± 0.001 and 3.64 ± 0.03 in units of 
e^, which is related to g'^ by Eq. (3.28)) in which the Greensite vacuum wave function of 
Eq. (6.37) was used. This serves as a check on our counting in calculating the possible 
overlaps of excited states. Our calculation is in disagreement with that of Arisue at the 
third eigenstate, for which Arisue calculates a mass (5.15 ± 0.1)e^. Our fourth eigenstate 
is close in mass to Arisue's third and our third eigenstate does not appear in his results. 
The reasons for this are not clear. 

The lowest mass unimproved SU(2) results (in units of e^) lies between the previ- 
ous Lagrangian Monte Carlo calculation of Teper, 1.582 ±0.017 [51], and the Hamiltonian 
strong coupling expansion result, 2. 22 it 0.05 [100]. It should be made clear that the strong 
coupling expansion result quoted here assumes that asymptotic weak coupling behaviour 
is already established in the range, 1.4 < /3 < 2.2. Data obtained from an extrapolation of 
the strong coupling expansion are fit by a straight line in that region to give the quoted 
result. For completeness, we must also mention that the coupled cluster estimates agree 
on an increasing glueball mass up to /? 2 but become variable after that [101, 102]. 
Convergence with increasing order appears to be a serious problem in such calculations. 

The results for the SU(3) symmetric massgap (in units of e^) are to be compared 
to calculations by Luo and Chen 2.15 ± 0.06 [103], Samuel 1.84 ± 0.46 [104] and Teper 
2.40 ± 0.02 [51]. Our result of 3.26520 ± 0.00009 is considerably higher than all exist- 
ing comparable results. By including more complicated loops in the expansion basis one 
would expect to reduce this estimate. This is emphasised by the fact that when using only 
square basis states our result is considerably higher. To explain the discrepancy between 
our results and others it is important to note that since we use a basis of rectangles, we 
exclude the contribution of many nonrectangular small area diagrams that are included 
in the calculations of Teper and that of Luo and Chen. For this reason it may be the 
case that what we have interpreted as the lowest glueball mass in this chapter may, in 
fact, be a higher order excited state. Teper has calculated the masses of the three lowest 
mass glueballs for SU(3) in the 0++ sector [51]: 0++, 0++* and 0++**, with the respec- 
tive results, in units of e^: 2.40 ± 0.02, 3.606 ± 0.063 and 4.55 ± 0.11. It is interesting 
to note that our result is closer to Teper's first excited state. In the same study Teper 
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Uiiiiiipro\'e(i Iiiipro\'ed Tadjjok^ Iiiipro\'e(i 

Am^ 2.05691 ± 0.00002 2.0897 ± 0.0003 2.0965 ± 0.0006 

Amf 3.645 ±0.001 3.685 ± 0.001 3.6953 ± 0.0009 

Amf 4.5202 ± 0.0004 4.574 ±0.004 4.583 ± 0.004 

Amf 5.133 ±0.003 5.177 ±0.004 5.189 ±0.004 

Amf 5.867 ±0.006 5.932 ± 0.008 5.943 ± 0.008 



Table 6.1: Estimates of the lowest lying SU(2) glueball masses (in units of e^) computed 
with various Hamiltonians in 2+1 dimensions. The unimproved, improved and tadpole 
results are calculated in the respective scaling regions 13.5 < (3 < 30.0, 9.9 < /3 < 30.0 
and 9.25 < (3 < 30.0. 



also calculated glueball masses in the 0"^"^ sector for iV = 4, 5 and 6. The mass, in units 
of e^, of his 0++* state for SU(4) is 4.84 ± 0.12 and 5.99 ± 0.16 for SU(5). We notice 
that as N is increased the results presented here move closer to the mass of Teper's 0++* 
state, with the improved results being closer than the unimproved. In fact for SU(5), the 
results presented here, improved and unimproved, are consistent with Teper's 0^^* mass. 
This forces us to question the interpretation of the large (3 plateaux in Fig. 6.8 as scaling 
regions for the lowest mass glueballs. It is possible that the minima present in Fig. 6.8 
in the small (3 region are possible scaling regions. It is possible that our vacuum wave 
function and minimisation basis are insufficient to extend this scaling region over a wide 
range of couplings and that as our approximation breaks down we observe a level crossing 
effect. We examine the possibility of the small (3 minima being scaling regions in Chapter 7. 



The antisymmetric results presented here can also be compared with those of Te- 
per [51]. While each of our results is considerably higher than the masses of Teper's 
and * states, Teper's ** state is close in mass to the lowest mass state calculated 
here. Teper obtains the following masses, in units of e^, for the ** state: 5.42 ±0.16 for 
SU(3), 6.98 ± 0.26 for SU(4) and 9.18 ± 0.45 for SU(5). It is interesting to note that our 
corresponding lowest unimproved glueball masses are consistent with these results. Our 
improved results show better agreement with Teper's ** state for SU(3) and SU(5) 
than the corresponding unimproved results. The improved SU(4) results are also consis- 
tent with Teper's SU(4) ** mass although the agreement is closer for the unimproved 
result. 



When compared to equivalent unimproved calculations, the improved and tadpole im- 
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Uiiiiiir)ro\'cxi 


Iiiipro\'txl 


ladpok} Iiin)rov(xl 




3.265868 ± 0.000042 


3.32365 ± 0.00012 


3.32580 ± 0.00015 


Amf 


6.23903 ± 0.00065 


6.30391 ± 0.00083 


6.31192 ± 0.00084 


Amf 


7.5767 ± 0.0025 


7.6466 ± 0.0030 


7.6498 ± 0.0030 


Amf 


8.9462 ± 0.0029 


9.0118 ± 0.0044 


9.0206 ± 0.0045 


Amf 


10.0778 ± 0.0071 


10.1546 ± 0.0094 


10.1628 ± 0.0094 



Table 6.2: Estimates of the lowest lying symmetric SU(3) glueball masses (in units of e^) 
computed with various Hamiltonians in 2+1 dimensions. The results are calculated in the 
scaling region which minimises the standard error in each case. 





Unimproved 


Improved 


Tadpole Improved 


Amf 


4.59121 ± 0.00007 


4.6720 ± 0.0001 


4.6754 ± 0.0001 


Amf 


8.8122 ± 0.0012 


8.9276 ± 0.0016 


8.9284 ± 0.0017 


Amf 


10.5889 ± 0.0051 


10.7807 ±0.0051 


10.7794 ± 0.0051 


Amf 


12.5527 ± 0.0048 


12.6266 ±0.0081 


12.6138 ± 0.0080 


Amf 


14.052 ±0.012 


14.165 ±0.016 


14.157 ±0.016 



Table 6.3: Estimates of the lowest lying symmetric SU(4) massgaps (in units of e^) com- 
puted with various Hamiltonians in 2+1 dimensions. The results are calculated in the 
scaling region which minimises the standard error in each case. 



proved massgaps approach scaling faster as (3 is increased. This is evident in Figs. 6.7 
and 6.8 and is expected since, for an improved calculation one is closer to the continuum 
limit when working at a given coupling. However, for most improved calculations the 
scaling behaviour is marginally less precise than the equivalent unimproved calculation. A 
possible reason for this is that the one-plaquette trial state used here does not allow for 
direct contributions from the improvement term in the kinetic Hamiltonian. For this term 
to contribute directly one would need a trial state which includes Wilson loops extending 
at least two links in at least one direction. 



The improved SU(2) massgap can be compared to the coupled cluster calculation of 
Li et al [105]. Their result (in units of e^), Amf = 1.59, is again significantly lower than 
our result 2.0897 ± 0.0003. The difference is again attributable to the different choices 
of Wilson loops used in the simulation of states. While our calculation makes use of the 
simple one-plaquette ground state and a minimisation basis with only rectangular loops, 
the coupled cluster calculation of Li et al uses a more accurate ground state wave function 
consisting of an exponential of a sum of extended loops which are not necessarily rectan- 
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Uiiiiiipro\'c}cl 


Iiiiprowxl 


ladpc^k} IiiiiJixnxxi 




5.8903 ± 0.0001 


5.99434 ± 0.00009 


5.9983 ± 0.0002 


Amf 


11.2335 ±0.0036 


11.3696 ±0.0050 


11.3731 ± 0.0049 


Amf 


13.340 ±0.011 


13.658 ± 0.012 


13.663 ±0.011 


Amf 


15.881 ± 0.012 


15.890 ± 0.019 


15.890 ± 0.019 


Amf 


17.564 ± 0.025 


17.676 ± 0.035 


17.682 ± 0.035 



Table 6.4: Estimates of the lowest lying symmetric SU(5) massgaps (in units of e^) com- 
puted with various Hamiltonians in 2+1 dimensions. The results are calculated in the 
scaling region which minimises the standard error in each case. 





Unimproved 


Improved 


Tadpole Improved 


Amf 


5.32750 ± 0.00047 


5.39661 ± 0.00034 


5.39864 ± 0.00032 


Am2 


7.9389 ± 0.0021 


8.0142 ± 0.0028 


8.0145 ± 0.0028 


Am3 


8.9319 ± 0.0045 


9.0092 ± 0.0056 


9.0087 ± 0.0055 


Ami 


10.4711 ± 0.0058 


10.5514 ±0.0085 


10.5502 ± 0.0085 


Ami 


11.304 ±0.011 


11.384 ±0.015 


11.381 ±0.015 



Table 6.5: Estimates of the lowest lying antisymmetric SU(3) massgaps (in units of e^) 
computed with various Hamiltonians in 2±1 dimensions. The results are calculated in the 
scaling region which minimises the standard error in each case. 



gular. Without including additional small area Wilson loops we cannot be confident that 
the lowest mass state accessible with our minimisation basis is in fact the lowest mass 
state of the theory. Clearly there is scope for more work here. 



6.7 Conclusion 

In this chapter we have extended the analytic techniques of 2+1 dimensional Hamiltonian 
LGT, traditionally used for SU(2), to general SU(Ar). Impressive scaling is achieved over 
an extremely wide range of couplings for the lowest energy eigenstates in the symmetric 

and antisymmetric sectors. Our calculations use a one-plaquette trial state and a basis of 
rectangular states over which excited state energies are minimised. Such choices allow the 
use of analytic techniques in SU(A^) calculations. 

The results of this chapter give estimates of the lowest unimproved, improved and 
tadpole improved SU(Ar) glueball masses all of which are above current estimates. We 
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Uiiiiiir)ro\'cx"l 


Impro\'txl 


TMdi)(^k^ InH)i(^vcd 


Amf 


7.21479 ± 0.0012 


7.3310 ±0.0011 


7.33586 ± 0.00077 


Ami 


10.9117 ±0.0033 


11.0099 ±0.004909 


11.0617 ±0.0046 


Ami 


12.121 ±0.007 


12.2779 ± 0.0088 


12.292 ± 0.009 


Ami 


14.5012 ± 0.0089 


14.592 ± 0.014 


14.574 ± 0.015 


Ami 


15.4521 ± 0.0173 


15.545 ± 0.023 


15.555 ± 0.023 



Table 6.6: Estimates of the lowest lying antisymmetric SU(4) massgaps (in units of e^) 
computed with various Hamiltonians in 2+1 dimensions. The results are calculated in the 
scaling regions which minimise the standard error in each case. 





Unimproved 


Improved 


Tadpole Improved 


Ami 


9.067 ± 0.003 


9.239 ± 0.002 


9.248 ± 0.002 


Ami 


13.717 ±0.008 


13.89 ±0.01 


13.89 ± 0.01 


Ami 


15.054 ±0.015 


15.32 ±0.02 


15.327 ±0.019 


Ami 


18.08 ± 0.02 


18.12 ±0.03 


18.12 ±0.03 


Ami 


19.084 ± 0.036 


19.19 ±0.05 


19.20 ± 0.05 



Table 6.7: Estimates of the lowest lying antisymmetric SU(5) massgaps (in units of e^) 
computed with various Hamiltonians in 2+1 dimensions. The results are calculated in the 
scaling regions which minimise the standard error in each case. 



suspect that the reason for the discrepancy is a lack of small area nonrectangular states 
in our minimisation basis. A basis of rectangular states was used for simplicity. The 
inclusion of nonrectangular states is straightforward and only complicates the counting of 
overlaps between diagrams of a particular type. When not including sufficient small area 
diagrams it is possibly that the lowest mass states of the theory are not accessible over a 
large range of couplings. A further improvement to our calculation would involve the use 
of an improved trial vacuum state. Such a state would possibly include several extended 
Wilson loops in its exponent. Without the development of new techniques for performing 
the required integrals the use of such a vacuum wave function would require the use of 
Monte Carlo techniques for the calculation of expectation values. In this scenario many 
advantages of the Hamiltonian approach presented here would be lost. 

In the next chapter we extend the calculations presented here to SU(25) in an attempt 
to explore the mass spectrum in the large N limit of pure SU(A^) gauge theory. 
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Figure 6.5: The unimproved, improved and tadpole improved lattice specific heat in 2+1 
dimensions for SU(2), SU(3), SU(4) and SU(5). 
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Figure 6.6: An order 8 calculation of the lowest SU(3) glueball mass in 2+1 dimensions 
using various functional forms for c(/3). The horizontal line is the continuum limit result 
obtained in Section 6.6. 
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Figure 6.7: The unimproved 2+1 dimensional symmetric massgaps for SU(2), SU(3), SU(4) 
and SU(5). 
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Figure 6.8: The lowest lying 2+1 dimensional symmetric massgaps for SU(2), SU(3) (both 
with L^ax = 25), SU(4) (with Lmax = 16) and SU(5) (with Lmax = 12). 
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Figure 6.9: The lowest lying 2+1 dimensional antisymmetric massgaps for SU(3) (with 
Lmax = 25), SU(4) and SU(5) (both with Lmax = 12). 



Chapter 7 

Large AT Physics 

7.1 Introduction 

In this chapter we study the glueball mass spectrum for large N pure gauge theory in 2+1 
dimensions. Our primary motivation lies in recent developments in string theory where a 
direct connection between pure gauge theory in the large limit and certain string theo- 
ries has been conjectured. An independent test of this conjecture is provided by large N 
calculations of the glueball mass spectrum in LGT. Impressive Monte Carlo calculations 
by Teper [51] and Lucini and Teper [106, 2] have progressed to the stage where accurate 
values of glueball masses have been calculated up to iV = 6. This has allowed an extrapo- 
lation to the N ^ oo limit. With the analytic techniques of the Hamiltonian approach it 
is possible to extend these calculations to much larger A*", allowing a more reliable N ^ oo 
extrapolation. 

The outline of this chapter is as follows. We start in Section 7.2 with a historical 
background, where the original motivations for large N physics are discussed and the recent 
developments in string theory introduced. In Section 7.3 we consider the exploration of 
the large limit directly by considering the N ^ oo limit of the expectation values that 
appear in the variational calculation of massgaps. We then move on to the extrapolation 
technique in Section 7.4, where various glueball masses are calculated at finite values of 
AT up to iV = 25 in some cases. These results are then extrapolated to explore the large 
N limit of the 2+1 dimensional glueball spectrum. The N ^ oo limits obtained here are 
compared to the results of the analogous Monte Carlo studies in the Lagrangian approach. 
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7.2 Background 

The nonperturbative region of QCD has been examined predominantly with numerical 
techniques. While the numerical Monte Carlo simulations of Lagrangian LGT have made 
significant progress, there have been few developments in analytic techniques. One ana- 
lytic technique which has received considerable attention in recent years is the large N 
limit. 

In 1974 't Hooft proposed a study of S\]{N) gauge theories in the large N limit rather 
than the physically interesting N = 3 case [107]. It was hoped that the SU(iV) theory 
could be solved analytically and would be, in some sense, close to SU(3). Indeed it was 
shown that the large N theory simplifies drastically [107] but still captures at least some 
of the complexity of the SU(3) theory [108, 109]. Motivated by the fact that the QCD 
coupling constant, g^, is a poor expansion parameter because its value depends on the 
energy scale of the process under consideration, t' Hooft considered an expansion in an- 
other dimensionless parameter of SU(A^) QCD, 1/A^. Based on an ingenious organisation 
of Feynman diagrams, 't Hooft was able to show that in the large N limit, keeping the 
't Hooft coupling, g'^N, fixed, only planar diagrams remain. While the remaining planar 
theory can be solved exactly in two dimensions, the three and four dimensional theories 
have not yet been solved. 

Recent developments in string theories have attempted to address this problem. In 
1997 Maldacena conjectured a correspondence, in the duality sense, between supercon- 
formal field theories and string theory propagating in a non-trivial geometry [72]. The 
idea behind duality is that a single theory may have two (or more) descriptions with the 
property that when one is strongly coupled the other is weakly coupled. A technique for 
breaking the conformal invariance and supcrsymmctry constraints of Maldacena's original 
proposal was later developed by Wittcn [73]. This led to the hope that nonperturbative 
pure SU(iV) gauge theories could be described analytically in 3-1-1 dimensions by their 
string theory dual. To develop this hope into a concrete solution would require perturba- 
tive expansions within the dual string theory. Such developments would seem to be some 
way off. 
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In order to test the viability of using string theories to probe nonperturbativc QCD 
independent tests are needed. The glueball mass spectrum of QCD provides a perfect 
laboratory. Recent Monte Carlo calculations have provided stable estimates of glueball 
masses up to AT = 6 [51, 2] in 2+1 dimensions and up to AT = 5 [2] in 3+1 dimensions. 
Prom these estimates an AT — >^ oo limit can be extrapolated in each case. Convincing 
evidence of 0(l/iV^) finite N corrections have been demonstrated in 2+1 dimensions, 
verifying a specific prediction of 't Hooft's 1/N expansion in the quarkless case. Prom the 
string theory side, present calculations require a strong coupling limit {N,g'^N oo) to be 
taken. In this limit the string theory reduces to classical supcrgravity and the calculation 
of glueball masses, m(J^*^), is straightforward. In this limit the glueball is represented 
by a dilaton field whose mass can be extracted by solving the dilaton wave equation [110, 
111, 112]. Despite the extreme approximations required surprising quantitative agreement 
with the weak coupling results from LGT were obtained. Later it was realised that, in 
2+1 dimensions, lower mass states can be constructed from the graviton field [113]. With 
this approach the quantitative agreement with LGT is lost but qualitative agreement in 
the form of the prediction 

m(0++) < m(2++) < m(l-+) (7.1) 

still holds. The qualitative agreement between string theories and LGT persists in 3+1 
dimensions [114] with 

m(0++) < to(2++) < m(0-+), (7.2) 

for both theories. It appears that the identification of glueball states remains a prob- 
lem. A serious difficulty is the removal of unwanted states which have masses of the same 
magnitude as the glueballs. Recipes for removing these unwanted states have been pro- 
posed [115, 116] but new parameters need to be introduced. 

The exploration of large gauge theories outside of string theory has not been re- 
stricted to traditional Monte Carlo methods. Recent studies using the light-front Hamil- 
tonian of transverse LGT have shown agreement with the Monte Carlo calculations of 
Teper [117]. In this approach explicit calculations at A/" — oo are possible without the 
need for extrapolation. Continuum calculations, without the use of string theory, have also 
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commenced. An impressive series of papers [118, 119, 120] has led to specific predictions 
for the string tension for all in 2+1 dimensions [120]. These predictions lie within 3% 
of lattice calculations up to iV = 6 beyond which LGT results are not available. A recent 
review is available in Ref. [121]. 



7.3 Direct Calculation of Asymptotics 

In this section we consider the direct calculation of the asymptotics of the generating func- 
tions defined in Chapter 5. In this way it is hoped that analytic forms for the asymptotics 
of the matrix elements on which the glueball masses depend will be obtained. We start 
with the generating function of Eq. (5.33). 



To calculate the asymptotic form of Eq. (5.33) in the large A^ limit we make use of 
a differential equation technique and some elementary linear algebra [122]. Consider an 
N X N matrix, Mni{x), depending on parameters x and I, for which we wish to calculate 
the determinant. Provided an N x N matrix Tni{x) exists such that 

^M^{x)=Tm{x)Mm{x), (7.3) 

the determinant of Mni{x) satisfies the simple first order differential equation 

ddeiMm 



dx 



- (x) = TV [Tmi {x)\ det Mmi {x). (7.4) 



The problem of finding the determinant reduces to one of finding a suitable Tni{x) and 
solving the first order differential equation in Eq. (7.4). 



Consider the N x N matrix Mpji{x) whose (i,j)-th entry is defined as follows 

[Mm{x)],^ = I±i^j-i{2V^). (7.5) 

Making use of Mathematica we calculate a suitable Tni{x) such that Mni{x) satisfies 
Eq. (7.3) for various N and I. In this way a general form for Tf^i{x) may be guessed. Once 
a form is guessed the proof that it satisfies Eq. (7.3) for all N and I is simply a matter of 
substitution. Calculations for AT = 2, . . . , 10 and I = 0, . . . , 6 suggest a general form for 
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Tni (x) which has a trace given by, 



X 



N 



TiTmix) 1 (^ + 1)1(^ + 2)! + (iV + l)!(7V + 3)! 
TrT^lix) = ^ + JL. + ... 1^0, (7.6) 

with the higher order terms vanishing in the N ^ oo hmit. To find an asymptotic form 
for det Mjsfi we must however do better than this and find the asymptotic form for the 
matrix Ti^i{x), not simply its trace. At this stage we have not found such an asymptotic 
form. Without such a result we can only conjecture an asymptotic form for detMjv«(x). 
However all is not lost. The success or failure of any asymptotic form that we conjecture 
can be determined numerically. 



Substituting Eq. (7.6) in Eq. (7.4) and solving the resulting differential equation leads 
to the asymptotic form, 



exp - 



Nx 



Here Cni does not depend on x. We determine Cat; by noticing that 



(7.7) 



Cm = limx-^l'l/^exp (-^^) detMmix). 



(7.8) 



Again we use Mathematica and calculate C^i for various A'^ and /. Some numerical exper- 
imentation leads to the conjecture that, in the large A^ limit, 

112!...(|i|-l)! ™|/2 



We denote the fractional error of this expansion by eiyi{x). To verify this result numerically 
we plot em for various N and I as a function of x. The typical examples of / = 2 and 
I = —4 with N = 1, . . . , 8 are shown in Fig. 7.1. Fast convergence over a large range of x 
ans N is increased is clear. Similar results are obtained for other values of L In each case 
it is clear that convergence hastens for N > Making use of Eqs. (5.33) and (7.9) we 
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Figure 7.1: A plot of ej\[i{x) versus x, for 1 = 2 and I = —4 demonstrating the eonvergence 
of the exact result detMjV2(a3) to the asymptotic form in the large A'" limit, Eq. (7.9). 



obtain the large N asymptotics of the SU(iV) generating function, 

1!2!...(U| -1)! 



Gsv{N) 



cd 



cj N\{N + 1)1... {N + 



(cd) 



N\l\/2 



m 



(7.10) 



We now move on to the second generating function of Chapter 5, that of Eq. (5.52). For 
this case we define 



PNi{x,y) = det [Ii+j-i{2x) + yli+j-i+2{'2x)] 



l<i,j<N ■ 



(7.11) 



For the purpose of calculating the matrix element {Z2) in the large N limit, we are inter- 
ested in the large N asymptotic form of 



d_ 
dy 



det Pm{x,y) 



(7.12) 



y=o 



Again we perform numerical experiments using Mathematica and conjecture that there 
exists a SNi{x,y) such that 



—Pmix,y) = Smix,y)PNiix,y), 



(7.13) 
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with 



TiSm{x,y) 

We then have, 

— det Pm{x,y) 
dy 



„2N 



y=o 



{N-1){N + 1) 



N\l\x^ 



(iV + |Z|-l)(iV + |Z|)(iV + |Z| + l) 



I ^ 0. (7.14) 



Tr [S mix, y)] det PNiix,y) 



y=o 
N\l\x'^ 



-1 



(TV + \l\ - l)iN + |/|)(7V + |Z| + 1) 1 J-^ (iV + m)\ 



n 



^- ^^^1^1 exp ( ^^llf^ ) . (7.15) 



N + 



Here we have recognised that Pni{x,0) = Mj^i{x) and used Eq. (7.9). Making use of 
Eq. (5.53) we can sum Eq. (7.15) over I to obtain the asymptotic form of the following 
SU(Ar) integral, 



„2N 



sv(N) ^ ' (iV-l)!(iV + l) 



2Nl(? 



l-i 



^(iV + /-l)(iV + 0(iV + / + l) jn (iV + m)! 



N + 



(Ar + l)(iV + 2)iV! 



0. 



(7.16) 



Compiling the above results we expect the following simple asymptotic forms for the 
matrix elements on which the glueball masses depend: 

{Z2) ~ 

{ZD ~ 

(ZiZi) ~ l + c^ 

{Zi) ~ c. (7.17) 

These limits are verified numerically in direct calculations of the matrix elements as A'^ is 
increased but only for small c. For finite N, glueball mass scaling occurs at large c, outside 
the range of c for which the asymptotic forms listed in Eq. (7.17) hold. The problem can 
be explained as follows. In the above we have taken the naive iV — 00 limit. The limit in 
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which we are interested is the one in which N ^ oo with g'^N held constant. In practice 
c is given a dependence on g"^ and hence the requirement that g^N remain fixed forces 
a functional dependence of c on N. For instance if we have c oc N/g"^ then c oc AT^ 
in the large N limit if g'^N is held constant. The large N limit of the matrix elements 
therefore cannot be analysed until the functional dependence of c on /3 is fixed. It would 
therefore seem that an accurate description of the expectation values arising the glueball 
mass calculation in the large N limit, requires many more terms than calculated here. 

7.4 Extrapolation 

7.4.1 Introduction 

It this section we calculate 2+1 dimensional SU(Ar) massgaps for various N with the 
method of Section 6.5, with the intention of exploring the large N limit of the results 
obtained. With this method we can extend the calculation of the lowest energy massgaps 
to SU(25) on a desktop computer. This should be compared to current results in the 
Lagrangian approach which are limited, at present, to SU(6). With large values of 
available a more reliable extrapolation to the N ^ oo limit is possible. 

For the minimisation process, we use the same basis of rectangular states as in Sec- 
tion 6.5. With this basis the states and J , with J = or 2, are accessible. The 
J = 2 states require the use of a minimisation basis that contains no square states so that 
the excited states are invariant under rotations by vr but not tt/2. The dependence of the 
massgaps on the coupling does not change significantly as N is increased. We do however 
observe the appearance of what could possibly be additional low (3 scaling regions as N is 
increased. 

7.4.2 Convergence with Imax 

In Section 6.3.3 the dependence of the variational parameter on the truncation of the ^-sum 
in Eq. (5.33) was considered. The calculation of the variational parameter depends only on 
the plaquette expectation value. The dependence of the massgaps on the truncation should 
also be explored as their calculation incorporates on additional matrix elements. We find 
that fast convergence is achieved as /max is increased for the O"*""^ massgaps. A typical 
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Figure 7.2: Example spin massgaps in units of ja demonstrating different convergence 
properties when truncating the Z-sum of Eq. (5.33) at I = iZmax- 



example is shown in Fig. 7.2(a), where the massgaps are almost indistinguishable up to 
/? = 80 on the scale of the plot. For the convergence is not as fast. A typical example 
is shown in Fig. 7.2(b). We find that in order to obtain convergence of the antisymmetric 
massgaps up to ^ 350 for SU(25) we need a truncation of Imax = 30. 



7.4.3 The Small /3 Minima 

The massgaps for > 2, in units of Ng'^ /a do not differ significantly as functions of ■ 
The lowest lying 0"*"+ glueball masses have the characteristic form shown in Fig. 7.2(a) for 
each N. We observe a minimum at small f3 and a scaling plateau at large f3. 



In this section we examine the minima occurring in the lowest O"^"*" and glue- 
ball masses and consider the possibility that they may correspond to a scaling region. 
It appears that for SU(3) this region has been interpreted as a scaling region by some 
authors [123, 124, 76]. From the scaling arguments of Section 6.5.4, by changing the 
functional dependence of the parameter appearing in the vacuum state the scaling region 
can be modified. In this way it is possible to choose the variational parameter's depen- 
dence on (3 such that the region in which the massgap takes its minimum is stretched 
over a large j3 interval. In the same way, the scaling plateau occurring for large f3 may be 
compressed. With this in mind, it is unclear which scaling region we should take to be 



118 



Large N Physics 



the correct one. For this reason we analyse both potential scaling regions. We start, in 
this section, with the small /? minima and consider the large /? plateaux in the next section. 



Fig. 7.3 shows the lowest lying SU(A'') O"*""*" glueball masses in 2+1 dimensions in units 
of g^N/a as functions of for various N G [3,25]. It is apparent that the minima 
depend only weakly on N. We see that the minima appear to approach a finite limit from 
below as N ^ oo. On the scale of Fig. 7.3 the minima corresponding to A?^ = 15 and 
N = 25 are barely distinguishable. 



Let us denote the minima in Fig. 7.3 by oAM*^/ {Ng"^) and consider them as a functions 
of N. Fitting aAM''/(Ng'^) to the model 

qAM^ 72 

lv^ = ^^ + iv^ ^'-'^^ 

for N > 8 gives best fit parameters 

71 = 0.83262 ± 0.00022 

72 = -0.97 ±0.13 

73 = 0.990 ±0.035. (7.19) 

We fit on the data iV > 8 to minimise pollution from next to leading order corrections. If 
we assert that the power of the leading order correction for finite N must be an integer 
power of 1/A^, then from Eq. (7.19) that order must be 2 in agreement with the predictions 
of large N QCD. Fitting to a model with 1/iV^ corrections for N > 8 gives 

M,(IV^)= 0.83256 ±0.00007 -5^155^5:55!?. (r.20) 

We can go further and attempt to find the next to leading order corrections by fitting to 
the model 

^ = ^' + N2+J^- (7-21) 
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The best fit parameters when fit on the whole data set are 

71 = 0.83287 ± 0.00011 (7.22) 

72 = -1.116 ±0.028 (7.23) 

73 = 1.715 ±0.085 (7.24) 

74 = 1.627 ±0.050 (7.25) 

which is in disagreement with integer power next to leading order corrections. A good fit 
however is achieved to the model with 1/iV^ next to leading order corrections: 

;..(Ar=) = 0.83233 ± 0.00018 - M^"^^"-"") ^-772^ a090 ^ 

The precise locations of the minima, aAM'^/{Ng'^), are plotted as a function of in 
Fig. 7.5 along with the fitted models of Eqs. (7.20) and (7.26). 



Prom Eq. (7.20), our best fit result for the ^ oo limit of a/^M" / {N g'^) is 0.83256 ± 
0.00007. The error here is purely statistical. We should, as always, expect a significant 
systematic error attributable to our choice of ground state and minimisation basis. This 
result may be compared to the Monte Carlo calculation of Lucini and Teper who obtain 
the N ^ CO limit of the lowest 0''"'' glueball mass as 0.8116 ± 0.0036 [106]. This is com- 
parable to the result presented here. The result of Lucini and Teper is based on a linear 
extrapolation to the N —>■ oo limit of 2 < A" < 6 data. Their result differs significantly 
from the one presented here at the leading order finite TV corrections. The correction term 
obtained by Lucini and Teper, — (0.090 ± 0.028)/A/"^, has the same sign but is significantly 
smaller than the one presented here. When our data is fit on the range 3 < A^ < 6 (we 
don't obtain a small (3 minimum for SU(2)) our slope is halved but is still significantly 
larger than that of Lucini and Teper. It should be pointed out that Lucini and Teper's 
calculation was performed in the Lagrangian approach in which the coupling is the so 
called Euclidean coupling gE- The Euclidean coupling and the Hamiltonian coupling, g^, 
are equal up to order g'^ corrections. The precise relation between the couplings is only 
known for small g^ [61], a case which does not apply for the small (3 minima. It appears 
that our simple calculation induces a level crossing whereby the lowest mass glueball state 
switches to a higher mass state beyond the small /? minimum. Although we do not have an 
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explanation for this, presumably, by including additional states in the minimisation basis 
or implementing a more complicated vacuum state, the level crossing would no longer 
appear and the small /? minima would extend into large /? scaling regions. This should be 
checked in further studies. A glueball mass extracted from a large (3 scaling region can 
be confidently compared to a corresponding Lagrangian calculation, for in that region of 
couplings the ratio of to is unity. 



We obtain small /3 minima for the massgaps but only for iV > 5. Plots of these 
massgaps with N > 5 are shown in Fig. 7.4. The behaviour is significantly different to 
that observed for the symmetric massgap. By A/" = 25 the minima do not appear close 
to convergence. Indeed if convergence is occurring at all, it is significantly slower than 
was observed for the 0^^ state. To explore this further, in Fig. 7.6 we plot the minima 
of Fig. 7.4, which wc again denote by aAM'^ /{Ng'^), as a function of 1/^/N. We observe 
linear behaviour in the large N limit. Fitting the N > 9 data to the model 

71 + TTZTT^, (7-27) 



Ng"^ " Ar73/2' 



gives best fit parameters 



71 = 0.41 ± 0.01 

72 = 1.25 ±0.02 

73 = 0.98 ±0.04, (7.28) 

which is consistent with 73 = 1. Fitting the > 9 data to a model with leading order 
finite N corrections starting at 1/ y/N gives 

, , 1.255 ±0.003 

uUN) = 0.41390 ± 0.00007 + = . (7.29 



As was done for the symmetric case we can attempt to go further and find the form of the 
next to leading order finite N corrections. Fitting the N >5 data to the model 



(7.30) 
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gives best fit parameters 



71 



0.410 ± 0.001 



(7.31) 



72 



1.2718 ± 0.0045 



(7.32) 



73 



-76 ± 23 



(7.33) 



74 



9.5 ±0.4 



(7.34) 



indicating vanishingly small next to leading order corrections for N > 5. The models vi 
and f2 are plotted against the N > 5 data in Fig. 7.6. 

The small (3 minima appear to converge to a non-zero N ^ oo limit. The conver- 
gence is significantly slower than for the small /? minima in the symmetric case. It is also 
clear that the leading order finite N corrections are not the expected 0{1/N^) but rather 
0{l/\/lV). The N —>■ oo limit of 0.41 it 0.01 is approximately half the corresponding 0"'""'" 
result. Such a state does not appear in the calculation of Teper. 

It is possible that the small /? minima are spurious. Dimensional analysis gives the 
expected scaling form but does not allow us to decide which scaling region is preferable, 
or if one is a lattice artifact. However, based on lattice calculations to date, we expect the 
0++ glueball to be lighter than the in the continuum limit. Spurious scaling regions 
have been observed in other lattice calculations [47] . It is clear that additional analysis is 
needed here. Again, an important step will be to include additional nonrectangular states 
in the minimisation basis. 

7.4.4 The Large /3 Plateaux 

Having considered the small /? minima as possible scaling regions we now move on to 
the large /? plateaux. This scaling region appears for all N and for all states considered 
and so its interpretation as a genuine scaling region is less dubious. Furthermore, the 
glueball mass results extracted from these scaling regions may be confidently compared to 
the corresponding Lagrangian calculations since the ratio of the Hamiltonian to Euclidean 
coupling is unity up to small 0{g^) corrections. 

We start with the 0+''" state for which the best scaling behaviour is obtained. We 
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calculate the five lowest lying massgaps corresponding to the five lowest glueball masses 
accessible with our choice of ground state and minimisation basis. We consider values of N 
in the range 3 < A'^ < 25. For each N considered we find a large P scaling plateau for each 
of the five lowest mass states. The lowest lying massgap is shown in Fig. 7.7 for a range 
of N. We observe that in units of Ng^/a the massgaps do not depend strongly on N and 
that in the scaling region they appear to approach a finite limit. Similar observations can 
be made for the four higher mass states obtained, although the scaling behaviour is less 
precise. We show the second lowest glueball mass in Fig. 7.8 as an example. A continuum 
limit for each massgap is obtained in the scaling region by fitting to a constant. For each 
fit we use a region of at least 10 data points which minimises the standard error. The 
continuum limit results for the five lowest lying massgaps, denoted by aAM'^/{Ng'^), are 
shown as functions of in Figs. 7.7 and 7.10. Also shown in the plots are fits to models 
with leading order large N corrections starting 1/A/"^, 

-r=p^+§2+w^- (7-35) 



Here the subscript i labels the i-th lowest glueball mass. The values of the best fit param- 
eters are given in Table 7.1. 
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Figure 7.4: L^ax = 6 SU(A'') lowest mass antisymmetric massgaps in 2+1 dimensions in 
units of g^N/a as functions of The Z-sum of Eqs. (5.33) and (5.52) truncated at 
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Pi 


Qi 
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1.23526 ±0.00023 


-1.540 ±0.018 


1.97 ±0.16 


2 


2.36924 ± 0.00072 


-2.478 ± 0.047 


-1.62 ±0.42 


3 


2.88446 ± 0.00071 


-3.435 ± 0.017 




4 


3.35422 ± 0.00047 


-3.476 ± 0.012 


Qa 


5 


3.7667 ± 0.0013 


-4.114 ±0.092 


1.88 ± 0.83 



"Set to zero to obtain a stable fit. 

Tabic 7.1: The best fit parameters for the five lowest energy O^"*" massgaps when fitting 
Eq. (7.35) to the available data. 



These results should be compared with those of Lucini and Teper who obtain large N 
glueball masses in units of Ng^/a in the 0++ sector with linear fits given by 

0.8116(36) - 
1.227(9) - 

0++** ; 1.65(4) - (7.36) 

The fit for tlie O'*'^** state is obtained using 4 < < 6 data. The remaining fits are 
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Figure 7.5: The continuum limit lowest mass symmetric 2+1 dimensional massgaps in 
units of Ng'^/a as a function of taken from the small /3 minima. The dashed lines 

are fits to Eqs. (7.20) and (7.26). 



obtained using 2 < iV < 6 data. We find that the lowest glueball mass extracted from 
our large /? plateaux is consistent, in the iV — ^ oo limit, with the state which Lucini and 
Teper label 0"*"+**. The slopes of the fits are of the same sign but differ significantly in 
magnitude. The 0'''"^** state of Lucini and Teper does not appear in our data in the form 
of a large f3 scaling plateau. There is however a hint of an approach to scaling in the 
vicinity of their result in our second lowest massgap data. This effect is only visible in our 
data for > 13 and occurs for small /? as shown in Fig. 7.11. The values of P for which 
this effect appears are quite close to those for which the small beta minima are observed 
in the lowest mass eigenstate. Similar effects are observed in the data but the effect 
is much less convincing with the data currently available. 



Let us consider this small /3 region as a possible scaling region. We fit a constant to 
the available A'^ > 15 data on a range of at least 4 data points chosen so that the standard 
error is minimised. These scaling values are plotted as a function of l/N"^ in Fig. 7.12. 
Also shown is the best fit linear model 

4+ = 1.7605 ± 0.0032 - (7.37) 
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Figure 7.6: The continuum limit lowest mass antisymmetric 2+1 dimensional massgaps in 
units of Ng'^/a as a function of 1/y/N taken from the small (3 minima. The dashed lines 
are fits to Eqs. (7.29) and (7.30). 



This produces an A?" — *■ oo limit which is close to, but inconsistent with, the result obtained 
by Lucini and Teper for their O"*""*"** state. 



We now move on to the states. For these states the scaling behaviour is not 
as precise as that obtained for the 0++ states. Again we calculate the five lowest lying 
massgaps corresponding to the five lowest mass glueballs accessible with our choice of 
minimisation basis and ground state. For each N considered up to 25 we find a large (3 
scaling plateau for each of the five states. The lowest lying massgap is shown in Fig. 7.13 for 
a range of N. The second lowest energy massgap is shown in Fig. 7.14. Considerably less 
data has been obtained for the states due to the large Zmax required for convergence as 
discussed in Section 7.4.2. Despite this, the results qualitatively replicate those of the 0^+ 
states. We observe that in units of Ng^ /a the massgaps do not depend strongly on N and 
that in the scaling region they appear to approach a finite limit. Similar observations can 
be made for the four higher mass states, although the scaling behaviour is less precise. The 
continuum limit values extracted are plotted as functions of 1/A/"^ in Figs. 7.15 and 7.16. 
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1.8896 ±0.0011 


-1.829 ± 0.068 


6.95 ± 0.62 


2 


2.8930 ± 0.0044 


-2.96 ±0.23 


6.6 ±2.0 


3 


3.20871 ± 0.0047 


-3.35 ±0.22 


8.8 ± 1.9 


4 


3.83647 ± 0.0024 


-3.41 ± 0.05 




5 


4.0759 ± 0.0019 


-2.45 ±0.12 


-7.4 ± 1.1 



"Set to zero to obtain a stable fit. 

Table 7.2: The best fit parameters for the five lowest energy massgaps when fitting 
Eq. (7.38) to the available data. 



The dashed lines show fits of the continuum limit values to the model 



with the parameters for each excited state given in Table 7.2. 



These results should be compared again with those of Lucini and Teper who obtain 
large glueball masses in units of Ng^/a in the sector with linear fits given by 

0" : 1.176(14) + "•"<^''> 



0—* : 1.535(28) 



Ar2 
0.35(35) 



0—** : 1.77(13) + (7.39) 

The result for ** is from Rcf. [51]. All fits were obtained using 3 < < 6 data. The 
N ^ oo limit of our lowest lying massgap extracted from the large P plateaux is consistent 
with Teper's ** result. 



Having considered spin states we now move on to spin 2, the only other spin accessi- 
ble when using a basis of rectangular states. For the case of spin 2 the scaling behaviour 
of the massgaps is significantly worse than for spin 0. More troublesome is the fact that 
the convergence of the 2++ massgaps with increasing /max is slower than the case of . 
The situation for the 2 state is markedly better with the convergence of massgaps with 
increasing Zmax being no different to that of 0++. For this reason we concentrate on the 
2 sector here. 
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3.20379 ± 0.00006 


-2.9571 ± 0.0032 


5.73912 ± 0.02585 


2 


4.07376 ± 0.00066 


-2.189 ±0.050 


-9.46 ± 0.43 


3 


4.96259 ± 0.00090 


-4.947 ±0.045 


1.99 ±0.35 


4 


5.26717 ± 0.00040 


-6.221 ± 0.031 


15.92 ±0.26 


5 


5.744 ± 0.024 


-5.13 ±0.56 


Qa 



"Set to zero to obtain a stable fit. 

Table 7.3: The best fit parameters for the five lowest energy 2 massgaps when fitting 
Eq. (7.40) to the available data. 



The 2 states produce less precise scaling behaviour than the spin states. However 
large P plateaux appear for each of the five lowest lying massgaps. We use these regions 
to estimate their respective continuum limits. For these states small f3 minima do not 
appear. The lowest lying 2 massgap is shown in Fig. 7.17 for a range of N. Once again 
we observe that in units of Ng'^/a the massgaps do not depend strongly on N and that 
in the scaling region they appear to approach a finite limit. Similar observations can be 
made for the four higher energy states obtained, although the scaling behaviour worsens 
as the mass of the state increases. The extracted continuum limit values are plotted as 
functions of in Figs. 7.18 and 7.19. The dashed lines show fits to the model 

'r=p^ + ^ + ir^■ (7.40) 

The best fit parameters for each excited state are given in Table 7.3. 



Once again these results should be compared with those of Lucini and Teper who 
obtain large N glueball masses in units of Ng^/a in the 2 sector with linear fits derived 
from 3 < iV < 6 data given by 

2- : 1.87(12) - (7.41) 

The 2 * result is from Ref. [51]. The N ^ oo results presented here are significantly 
higher. This time no correspondence between our results and those of Lucini and Teper 
can be obtained. 
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Figure 7.7: Lowest mass symmetric 2+1 dimensional massgaps in units of Ng^ /a as a 
function oil/g^. 



Wc finish this section by presenting tlic mass spectra obtained for tlie 0^+, and 
2 sectors. The results are summarised by the plots in Figs. 7.20, 7.21 and 7.22. 



7.4.5 An Empirical Observation 

In this subsection we present an empirical observation that could possibly be useful in the 
construction of simple glueball models. Since Lucini and Teper include many small area 
Wilson loops in their construction of excited states, we assume that they calculate the 
lowest few glueball states without omission. We assume the calculations presented here 
give higher mass glueball states. Presumably, the correct enumeration of these states is 
not possible without the inclusion of nonrectangular states in the minimisation basis. For 
example, the third lowest mass state in a given sector calculated here could be the true 
seventh lowest mass state. However, when our results are combined with those of Lucini 
and Teper an interesting empirical observation can be made. We can choose a labelling of 
the 0"*""^, and 2 excited states such that the large N limit of their masses lie on a 
straight line as shown in Fig. 7.23. We have included the results of Lucini and Teper for 
, *, 2 and 2 * which do not correspond to any of the states calculated here. 
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The straight Unes are fits to the model, 

mn(J^^)=7i(2n + 72) (7.42) 

where 71 and 72 are parameters and 72 is restricted to integer values. For the J^^ states 
considered we obtain the following best fit models 

m„(0++) = (0.256 ± 0.002) (2n + l) 
mn(0~) = (0.151±0.002)(2n + 5) 

mn(2~) = (0.1495 ±0.0008)(2n + 7). (7.43) 

Here n > 1 labels the n-th lowest mass state. We note that the fit improves in accuracy 
as n is increased. It is interesting to note that the slopes of the and 2 models are 
consistent suggesting that the constant of proportionality does not depend on J. Another 
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Figure 7.9: Continuum limit lowest mass 0"''^ SU(A^) massgap in units of Ng^/a as a 
function of 1/N'^. The dashed line is the fit to the quadratic model of Eq. (7.35). 



interesting observation can be made by recasting the models in the form 

m„(0++) = (0.256 ± 0.002) (2n + + 1) 
m„(0— ) = (0.151 ±0.002)(2n + 4 + 1) 

m,„(2— ) = (0.1495 ±0.0008)(2n + 6 + 1). (7.44) 

We notice a similarity with the two dimensional harmonic oscillator spectrum, 

En0^2n + J + 1, (7.45) 

when we take into consideration the ambiguity modulo 4 of spin identification on the 
lattice. With this in mind we propose a simple model for the J^^ spectrum 

^n(j''^) = Cpc [2n + 7(J^^) + 1] , (7.46) 

where Qpc is a spin independent parameter and 7(J^*-^) is an integer for which 7(J^*^) = 
J mod 4. Prom Eq. (7.43) we have C— ~ 0.15 and C++ = 0.256 ± 0.002. To check this 
simple model we can attempt to predict the lowest lying states obtained by Lucini and 
Teper in the sectors that have not been considered in this chapter. We start with the 2++ 
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sector in which Tcpcr and Lucini obtain 1.359(12) and 1.822(62) for the N ^ oo hmit 
of the 2++ and 2++* glueball masses in units of Ng^ /a. The predictions of Eq. (7.46), 
with 7(2++) = 2 and C++ = 0.256, are 1.28 and 1.79 for the 1st and 2nd lowest glueball 
masses in units of Ng^/a. As we would expect the prediction is better for the higher mass 
state. We can also consider the 1+"'" sector which has an equivalent spectrum to 1 ^ due 
to the phenomenon of parity doubling [51]. Lucini and Teper obtain 1.98(8) for the mass 
of the 1++ glueball in units of Ng'^/a. The model of Eq. (7.46), with 7(1"^+) = 5 and 
("++ = 0.256, gives 2.048. In the 1 sector the agreement is not as good, with Lucini and 
Teper obtaining 1.85(13) for the mass of the 1 glueball and the model of Eq. (7.46), 
with 7(1"") = 9 and C— = 0.15, giving 1.812. 

It will be interesting to see if the model suggested here stands up to further calculations 
with an extended minimisation basis and in other J^'~^ sectors. The analytic techniques 
presented here are at a great advantage to the standard Monte Carlo techniques of La- 
grangian LGT for the purpose of testing glueball models. High order excited states are 
readily accessible; in a variational calculation with s states in the minimisation basis, s 
glueball states are accessible. This is in contrast to competing Lagrangian calculations in 
which only 3 states are currently accessible in some sectors. 

7.4.6 Discussion 

In this section we have calculated variational mass spectra for pure SU(A'') gauge theory 
with a simple basis of rectangular states in 2+1 dimensions. Such a basis is easy to work 
with computationally and is therefore a good starting point. However to accurately explore 
the pure gauge mass spectra additional states need to be included. Perhaps a more suit- 
able basis would be the set of all states with an area less than some maximum value which 
would define the order of the calculation. The next stage beyond this would be to include 
higher representation states in the form of Wilson loops in which some links are covered 
more than once. In this way multiple glueball states could presumably be explored. The 
inclusion of additional small area states in the basis, in particular states which are not 
symmetric about reflections in coordinate axes, would also allow the calculation of spin 1 
glueball masses and also massgaps in the P = —C sector. 

It is important to point out that the inclusion of additional states in our minimisation 
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basis does not present a significant challenge. The only complication would be in counting 
the possible overlaps of particular states. While more difficult than the counting required 
here, the process could presumably be automated using symbolic programming and tech- 
niques from graph theory. The analytic techniques used here would still be applicable 
until higher representation states were included in the minimisation basis. At that stage 
further character integrals would be required to handle the additional integrals arising in 
the calculation. 

A final extension of the method presented here would be to make use of an improved 
ground state; a ground state which includes additional Wilson loops in the exponent. How- 
ever, the advantage of using analytic techniques would then be lost, unless new technology 
was developed for tackling the resulting integrals. One would need to calculate the re- 
quired matrix elements via Monte Carlo simulation on small lattices, losing the advantage 
of working on an infinite lattice with analytic expressions that we have here. 

7.5 Conclusion 

In this chapter we have applied the analytic techniques developed in Chapter 5, in a study 
of the large N glueball mass spectrum in 2-1-1 dimensions. 

A direct attempt at calculating asymptotic expansions of the generating functions was 
carried out in Section 7.3. While some low order large N expansions were obtained, many 
terms appear to be needed for a direct N ^ oo calculation of the matrix elements appear- 
ing in glueball mass calculations at small couplings. 

Having tried a direct approach to calculating matrix elements in the N ^ oo limit 
without success, an extrapolation technique was used. In Section 7.4 we calculated glue- 
ball masses at finite N, with N as large as 25 in some cases. This allowed accurate 
N ^ oo extrapolations to be made. Evidence of leading order finite N corrections 

to the glueball masses was obtained, confirming a specific prediction of large N gauge 
theory. The interpretation of the possible scaling regions was discussed and agreement 
with the Lagrangian study of Lucini and Teper [106] was obtained in some cases in the 
AT — oo limit. The discrepancies are attributable to our use of only rectangular states in 
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the minimisation basis. Without the inclusion of smah area, nonrectangular states, it is 
hkely that some low energy states are inaccessible. Further work will fill in the incomplete 
spectra presented here. Interesting empirical observations were made when the results 
presented here were combined with those of Lucini and Teper; the enumeration of excited 
states can be chosen so that the glueball mass spectrum has the structure of a two dimen- 
sional harmonic oscillator. To develop this observation into a model would require a more 
complete calculation of the mass spectrum with nonrectangular states in the minimisation 
basis and possibly a more complicated vacuum trial state. 

Clearly much work remains to be done before an accurate picture of the 2+1 dimen- 
sional pure SU(A/") gauge theory spectrum is achieved within a Hamiltonian variational 
approach. We have demonstrated however that the analytic techniques of Hamiltonian 
LGT can be used for massgap calculations with N as large as 25 on a desktop com- 
puter. This is significantly closer to the N ^ oo limit than is currently possible with 
the Lagrangian approach using supercomputers. The inclusion of additional states in the 
minimisation basis, while not presenting a major challenge, will allow a thorough study of 
the 2+1 dimensional pure SU(A'^) mass spectrum at least up to A'^ = 25. 

In the next chapter we consider the extension of the analytic techniques used with 
success in 2+1 dimensions here, to the physically interesting case of 3+1 dimensions. 
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Figure 7.10: Continuum limit 0"'""'' SU(iV) massgaps in units of Ng^ /a as functions of 
1/N'^. The dashed lines are fits given in Eq. (7.35). 
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Figure 7.11: Second lowest mass symmetric 2+1 dimensional massgaps in units of Ng^ /a 
as a function of l/{Ng^). 
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Figure 7.12: Continuum limit extrapolations (in units of Ng^/a as a function of l/N"^) 
derived from the low /3 scaling region that appears for iV > 13 in the second lowest mass 
0++ SU(iV) massgap. The dashed line is the fit to the quadratic model of Eq. (7.37). 
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Figure 7.13: Lowest mass 2+1 dimensional massgaps in units of Ng'^/a as a function 
of l/{Ng^). 
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Figure 7.14: Second lowest mass 2+1 dimensional massgaps in units of Ng^ /a as a 
function of ■ 
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Figure 7.15: The 2+1 dimensional continuum limit lowest SU(iV) glueball mass in 
units of Ng'^/a as a function of The dashed line is the fit to the quadratic model of 

Eq. (7.38). 
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Figure 7.16: Continuum limit SU(iV) massgaps in units of Ng^ /a as functions of 
1/N'^. The dashed lines are fits given in Eq. (7.38). 
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Figure 7.17: The lowest lying 2+1 dimensional 2 SU(7V) massgaps in units of Ng'^/a 
as a function of 1/^^. 
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Figure 7.18: The 2+1 dimensional continuum limit lowest 2 SU(A^) glueball mass in 
units of Ng^/a as a function of The dashed line is the fit to the quadratic model of 

Eq. (7.40). 
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Figure 7.19: The 2+1 dimensional continuum limit 2 SU(A^) massgaps in units of Ng^ /a 
as functions of 1/iV^. The dashed lines are fits given in Eq. (7.40). 
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Figure 7.20: An estimate of the continuum limit mass spectrum for 0"*"^ SU(Ar) glueballs 
in units of Ng^/a as a function of The Hues are fits to the models of Eq. (7.35). 
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Figure 7.21: An estimate of the continuum limit mass spectrum for SU(A^) glueballs 
in units of Ng^ /a as a function of l/N"^. The lines are fits to the models of Eq. (7.38). 
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Figure 7.22: An estimate of the continuum limit mass spectrum for 2 SU(A'') glueballs 
in units of Ng^/a as a function of The hues are fits to the models of Eq. (7.40). 




Chapter 8 

3+1 Dimensions 



8.1 Outline 

In this chapter we explore the viabihty of extending the analytic techniques used with 
success in Chapters 6 and 7 to the calculation of glueball masses in the pure gauge sector 
in 3+1 dimensions. The primary difficulty lies in the calculation of expectation values in 
3+1 dimensions. In Section 8.2 we briefly review what has been achieved in Hamiltonian 
LGT in 3+1 dimensions. We discuss the difficulties faced in 3+1 dimensions and possible 
solutions in Section 8.3. In Section 8.4 we consider the problem of Gauss' law constraints. 
This is a topic that has been discussed in the context of Hamiltonian LGT most recently 
by Ligterink, Walet and Bishop [24] and concerns the constraint equations that appear 
when non-abclian gauge theories are canonically quantised. In Section 8.5 we move on to 
the calculation of variational glueball masses on a single cube using the analytic variational 
technique discussed in Chapters 6 and 7. We finish in Section 8.6 with a discussion of the 
viability of pursuing analytic techniques for pure SU(iV) LGT in 3+1 dimensions based 
on the results of Section 8.5. 

8.2 Introduction 

Prom a renormalisation point of view the key difference between 2+1 and 3+1 dimensional 
gauge theory lies in the units of the coupling constant. As was seen in Chapter 3, 2+1 
dimensional gauge theory has a coupling constant, e^, with the dimensions of mass and so 
the coupling constant explicitly sets a mass scale for calculations on the lattice. In con- 
trast the 3+1 dimensional coupling constant is dimensionless. This makes the extraction 
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of continuum physics from lattice calculations more subtle in 3+1 dimensions than in 2+1, 
as discussed in Section 2.6. 

On a practical level, in Hamiltonian calculations there is a more serious problem faced 
in moving from 2+1 to 3+1 dimensions. The analytic techniques that we have used with 
success in Chapters 6 and 7 are no longer applicable. These techniques rely heavily upon 
the fact that in 2+1 dimensions a change of variables from links to plaqucttcs has unit 
Jacobian. The form for the equivalent Jacobian in 3+1 dimensions is considerably more 
complicated. The most comprehensive study of the change of variables from links to pla- 
quettes is due to Batrouni [125, 126]. We discuss this change of variables in more detail 
in Section 8.3. 

The extension of the techniques used in Chapters 6 and 7 to 3+1 dimensions is not 
straightforward. There are a number of immediate problems. Firstly, since plaquettes are 
not independent variables in 3+1 dimensions one can not automatically work in the infinite 
volume limit. In a precise study one would need to calculate identical quantities on dif- 
ferent sized lattices and extrapolate to the infinite volume limit. Secondly, in the context 
of analytic calculations, even on small lattices the integrals involved in the calculation of 
basic matrix elements are considerably more complicated than those encountered in 2+1 
dimensions. Such matrix elements could in principle be carried out analytically on small 
lattices but since the number of integration variables increases quickly with the volume of 
the lattice a calculation on even a 5^ lattice would seem exceedingly difficult. How quickly 
the infinite volume limit is reached will therefore determine the worth of pursuing analytic 
Hamiltonian methods in 3+1 dimensions. Finally, there is the complication of Gauss' law 
which we discuss in Section 8.4. 

The only Hamiltonian techniques to have been applied with any success to the case 
of SU(3) gauge theory in 3+1 dimensions have been strong coupling expansions, the t- 
expansion and exponential wave function methods. Each of which we now summarise. 

Strong coupling perturbative techniques were used in the early days of LOT in an 
attempt to bridge the gap between the strong and weak coupling limits. Strong coupling 
expansions of the Callan-Symanzik (3 function were calculated to 0{g~'^'^) [127, 16] and 
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showed signs of interpolating smoothly between the strong and weak coupling limits. This 
suggested that the /? function was a smooth function of the coupling with its only zero at 
g = 0, providing a strong argument at the time for the continuum limit of LGT confining 
quarks. Corresponding strong coupling expressions for glueball masses did not share the 
same success. Despite strong coupling calculations to 0{g~^^) [128] scaling was not ob- 
served in 0++, 1'' or 2++ glueball masses. 

The t-expansion was introduced by Horn and Weinstein [129] as an analytic method 
suitable for the study of LGT in the Hamiltonian formulation. It has been applied in the 
calculation of glueball masses in 3+1 dimensions for SU(2) [70] and SU(3) [130, 131, 132] 
LGT in the pure gauge sector. More recently it has been used in an attempt to calculate 
the lowest hadron masses [133]. For each case, however, asymptotic scaling of masses was 
not directly observed. Extrapolation techniques such as Pade approximants were required 
to probe the weak coupling region. The extrapolated values mass ratio results agreed with 
Monte Carlo estimates of the time. 

The coupled cluster method and related exponential wave function techniques have 
received by far the most attention in Hamiltonian LGT in recent years. Essentially these 
techniques aim to solve the Kogut-Susskind eigenvalue equation by making a suitable 
ansatz for the wave function. The coupled cluster method was originally constructed with 
applications in nuclear physics in mind [22, 23] but has since found the majority of its 
applications in molecular physics [134]. Its application in the context of Hamiltonian LGT 
is described in Refs. [135] and [25]. The truncated eigenvalue method, developed by Guo, 
Chen and Li [136], is another exponential wave function technique to have found applica- 
tion in Hamiltonian LGT. 

A number of groups have made considerable progress in the application of exponential 
wave function techniques to Hamiltonian LGT in 3+1 dimensions. While most studies 
have explored gauge groups other than SU(3) in less than three dimensions, studies of 
SU(3) glueballs in 3+1 dimensions have commenced. Results tangent to the expected 
scaling form, indicating an approach to scaling, have been obtained for SU(3) pure gauge 
theory in a coupled cluster calculation by Leonard [76]. However convergence with in- 
creasing orders appears to be slow. Finite order truncation errors appear to be under 
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more control in the truncated eigenvalue method. The first calculations of 3+1 dimen- 
sional SU(3) glueball masses with this method [137] gave a 1"*"" to O"*""*" mass ratio that 
was consistent with the Monte Carlo results of the time. The agreement was not as 
good for the to 0++ mass ratio. A convincing demonstration of asymptotic scaling 
has not yet been produced in a calculation of glueball masses for SU(3) in 3+1 dimensions. 

More promising 3+1 dimensional results have been obtained for higher dimensional 
gauge groups. Ironically these results have been obtained with much simpler methods 
than either of the exponential wave function methods described above. Asymptotic scaling 
of the lowest 0++ glueball mass has been demonstrated by Chin, Long and Robson in a 
variational calculation on a small volume (6^ sites) lattice for SU(5) and SU(6) [138]. 
This calculation follows the variational technique described in Chapter 6 but uses Monte 
Carlo rather than analytic techniques to calculate the required expectation values. It uses 
only plaquette states in the minimisation basis rather than the large basis of rectangular 
states used in Chapters 6 and 7. Naively one would hope that the same method could be 
explored on similar sized lattices for larger N, with additional states in the minimisation 
basis, using the analytic techniques of Chapters 6 and 7. We take this as our motivation 
for the studies presented in this chapter. It will become clear however that attempting a 
similar calculation to Chin, Long and Robson on a single cube using analytic techniques 
presents a significant challenge. 

8.3 The Move to 3+1 Dimensions 

In LGT one usually works in 2+1 dimensions to test a technique with the intention of 
later extending it to the physically relevant 3+1 dimensions. This may be justified in 
the Lagrangian approach where the time coordinate is treated on the same footing as 
the spatial coordinates and adding another dimension equates to nothing more than an 
increased load on computer memory. However, in Hamiltonian LGT two serious technical 
differences exists between 2+1 and 3+1 dimensions. The first is Gauss' law and the sec- 
ond is related to constructing a Jacobian for transforming from link to plaquette variables. 

To understand the first difference one needs to recall that Hamiltonian LGT is for- 
mulated in the temporal gauge which sets Aq = 0. If one starts with the Yang-Mills 
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Lagrangian and performs the standard equal-time quantisation, one runs into problems 
because the time derivative of Aq does not appear in the Lagrangian. The variational 
principle gives equations of motion for the space-like components resulting in the standard 
Yang-Mills Hamiltonian. For the time-like component we obtain a set of algebraic con- 
straint equations which are the analogue of Gauss' law. There is one constraint equation 
for each colour component of Aq, — 1 for SU{N), at each lattice site. As has been 
pointed out in the context of Hamiltonian LGT by Ligterink, Walet and Bishop [24], it is 
only when one works with a set of variables whose number of degrees of freedom matches 
the number of unconstrained degrees of freedom in the theory, that one can avoid the tech- 
nicalities of constraint equations. For this reason, in 2-1-1 dimensions we really are quite 
lucky. The number of plaquette variables on a square two dimensional lattice is precisely 
equal to the number of unconstrained variables [24] . For a three dimensional cubic lattice 
this is not the case. One could envisage constructing a polyhedral lattice in d dimensions 
such that the number of faces (plaquettes) would equal the number of unconstrained vari- 
ables. However such lattices appear to be prohibited by Euler's equation which relates the 
number of vertices, edges and faces of polyhedra. 

The other serious technical difference between 2-1-1 and 3-1-1 dimensions, that of con- 
structing a Jacobian for transforming from link to plaquette variables, is relevant to both 
the Hamiltonian and Lagrangian formulations of LGT. It only becomes important when 
the method of choice relies on plaquette variables for its implementation. The transfor- 
mation from link variables to plaquette variables is required in the approach taken in 
Chapters 6 and 7 to make use of the analytic results available for certain group integrals. 
Such a transformation does not need to be made if one is happy to use Monte Carlo tech- 
niques to handle the integrals as is the case in Ref. [138]. 

The most complete treatment of the transformation from link to plaquette variables 
is due to Batrouni [126, 125] who worked in the Lagrangian formulation. His approach 
was based on the continuum work of Halpern [139] who constructed a field-strength for- 
mulation of gauge theory. In Halpern's construction the definition of the field strength 
Ffj^iy is inverted to give an expression for as a function of the field strength. To do this 
requires the choice of a suitable gauge. The Jacobian of the transformation is precisely 
the Bianchi identity; a constraint equation on For the Abelian case the Bianchi 
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identity is equivalent to the requirement that the total magnetic flux leaving a volume is 
zero. Batrouni developed an equivalent field strength formulation on the lattice. On the 
lattice, the link operators, Ui, correspond to the vector potentials, A^, and the plaquette 
variables correspond to the field strengths, F^^. Batrouni demonstrated that the Jacobian 
of the transformation from link variables to plaquette variables was the lattice analogue 
of the Bianchi identity. For Abelian gauge theories the lattice Bianchi identity can be 
separated into factors which depend only on the plaquette variables of elementary cubes, 
with one factor for each cube of the lattice. For non-abclian theories the Bianchi identity 
has only been found to separate in this way for special types of lattices, the largest volume 
example being an infinite tower of cubes. For the infinite lattice the Bianchi identity is 
a complicated nonlinear inseparable function of distant plaquette variables. Interestingly 
it is the only source of correlations between plaquette variables in LGT. Mean plaquette 
methods have been developed to deal with the added complications of the non-abelian 
Bianchi identity but have not progressed far [140]. 

To summarise, Hamiltonian LGT in 3+1 dimensions faces some serious problems. Care 
must be taken in choosing appropriate variables to work with if constraints on the lattice 
electric fields are to be avoided. Additionally, if one wishes to use analytic techniques to 
calculate matrix elements the complications of the Bianchi identity restricts the calcula- 
tions to small lattices. 

8.4 Constraint Equations 

In this section we focus on one of the problems faced in moving from 2+1 to 3+1 di- 
mensions, that of Gauss' law. Prom the discussion of Section 8.3 it would seem that one 
needs either to construct an appropriate set of variables with precisely the correct number 
of unconstrained degrees of freedom or be faced with the problem of building constraint 
equations on the lattice electric fields. As has been pointed out by Ligterink, Walet and 
Bishop [24] there is an alternative. The problem of satisfying Gauss' law can be solved 
by working with wave functions that are annihilated by the generator of Gauss' law. We 
discuss this matter in that which follows. 

In the Hamiltonian formulation of LGT we work in the temporal gauge which sets 
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Aq = 0. In terms of link variables this is equivalent to setting all time-like links to the 
identity. This however is only a partial gauge fixing; setting = does not completely 
eliminate the arbitrariness in the definition of the vector potential, it merely reduces the 
set of possible gauge transformations to purely space-like ones. To see this we revert to 
continuum gauge theory. 



In the continuum one encounters difficulties when attempting to canonically quantise 
SU(Ar) pure gauge theory [141, 142]. The difficulties stem from the fact that it is impossible 
to define a momentum conjugate to the time-like component of the vector potential Aq. 
This follows from the non-appearance of OqAq in the Yang-Mills Lagrangian. The simplest 

way around the problem is to fix to the temporal gauge Aq = before quantising. The 
result is a Hamiltonian whose equations of motion do not mirror those of the Lagrangian; 
the non-abelian generalisation of Gauss' Law, 



is absent. One does not need to impose Gauss' law as a constraint on the electric fields, 
= Defining G = V ■ E, with the usual canonical commutation relations, it can 
be shown that commutes with the Hamiltonian. Consequently if one chooses to work 
with states such that 



then one recovers Gauss' law ensuring that the physics of the Lagrangian formulation is 
recovered. 



On the lattice a similar problem arises [24] . One must therefore ensure that any state 

used is annihilated by the generator of Gauss' law, in order to recover physical results with 
certainty. The generator of Gauss' law on the lattice can be written as [12, 16], 



DjF^^ = V -£ = 0, 



(8.1) 



IV) = 0, 



(8.2) 



i 



(8.3) 
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In this notation the lattice electric fields satisfy the commutation relations [16] 

[St{x),U^{y)] = X''Ujiy)6ijd^y (8.4) 
Sli{x + ia),Uj{y)\ = -Uj{y)X''Sij5^y. (8.5) 

It should be pointed out that in this notation we have uj{x) = U-i{x + ia). For physical 
states we must therefore have Q^{x)\iIj) = for each lattice site, x, and all a = 1, ... , N'^—1. 
It should be checked that this is the case for the one-plaquette exponential trial state. 
Consider first a state consisting of a single plaquctte acting on the strong coupling vacuum, 



\pij{x)) = Tr Ui{x)Uj{x + ia)U-i{x + ja)C/_j (x) |0) 



(8.6) 



Since |0) is annihilated by the electric field, using Eqs. (8.4) and (8.5) we immediately 
have 



g»(y)|p,,(x)) = 0. 



(8.7) 



for all lattice sites, y, not lying on the corners of the plaquette, pij(x). Consider now sites 
lying on the corners of the plaquette in question. In particular consider y = x. Making 
use of the commutation relations of Eqs. (8.4) and (8.5), 



g'^{x)\mj{x)) 



= Tr 

-Tr 
= 0. 



X''Ui{x)Uj{x + ta)U-i{x + ja)U-j{x) |0) 
Ui{x)Uj{x + ia)U-i{x + ja)U-j{x)X"-] |0) 



The same applies for other sites on the plaquette. This argument is not specific to plaque- 
ttes. Gauss' law is found to be satisfied locally by any closed Wilson loop on the lattice, 
traced over colour indices, acting on the strong coupling vacuum. It is easy to extend this 
result to products of such loops acting on |0). To see this we consider how Gauss' law 
applies to the product of two plaquettes 

\pij{x)pij{x + ai)) = pij{x)pij{x + ai)\0). (8.9) 

Only at the sites x + ai and x + a{i + j) does this case differ from the single plaquette 
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example . Let us consider Gauss' law at a;+ = x + ai. Once again, using the fact that the 
lattice electric field annihilates the strong coupling vacuum, and Eqs. (8.4) and (8.5) we 
have 

g%x+)\pij{x)pij{x+)) = Y,[£^{x+) + £lk{x+),Pij{x)pij{x+)]\(}) 

k 
k 

+Pij{x) [£^{x+)+£lk{x+),pij{x+)] |0)} 
= 0. (8.10) 

This result is easily extended to arbitrary products of closed Wilson loops, traced over 
colour indices, acting on |0). Consequently, any function of such loops acting on |0) also 
satisfies Gauss' law provided it admits a Taylor scries expansion. The one plaquette trial 
state thus obeys Gauss' law and is therefore suitable for use in simulating physical states. 

8.5 The One Cube Universe 

In this section we adapt the analytic techniques used with success in Chapters 6 and 7 for 
use in 3+1 dimensions. As a starting point we consider the case of a single cube. This 
will serve as a test case to assess the feasibility of a larger volume study. Our aim is to use 
a small basis of states to calculate 3+1 dimensional SU(A^) glueball masses variationally, 
and check if an approach to the correct scaling form is observed. We do not expect to 
achieve scaling with such a small lattice. However, an approach to scaling would warrant 
an extended study on larger lattices. We start with a general description of the approach. 
The starting point is the choice of trial state. As in 2+1 dimensions we choose to work 
with the one-plaquette trial state of Eq. (4.22), with the variational parameter again being 
fixed by minimising the vacuum energy density of Eq. (6.2). In 3+1 dimensions however 
the expression for the plaquette expectation value in terms of the variational parameter 
is significantly more complicated than for the case of 2+1 dimensions. Having fixed the 
variational parameter we then construct a small basis of states, with each state fitting on 
a single cube, and minimise the glueball mass over this basis. The process of minimising 
the massgap follows precisely Chapter 6. The key difficulty of working in 3+1 dimensions 
is the calculation of the required integrals. In the next section we explain how an analytic 
approximation to these integrals can be obtained. 
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8.5.1 SU(A7^) Integrals in 3+1 Dimensions 



The analytic techniques used in Chapter 6 rely on the fact that the transformation from 
link to plaquette variables has unit Jacobian in 2+1 dimensions. Batrouni has calculated 
the Jacobian for arbitrary numbers of dimensions [125, 126]. For a lattice consisting of a 
single cube, the result of Ref. [125] is 

= Einx^(^^)' (8-11) 

r i=l 

where Pi,. . . ,Pq are the six plaquette variables on the single cube and the sum is over 
all characters Xr of SU{N). The second line is simply a character expansion of the first 
line. In the variational study of glueball masses, we need to calculate the integrals of 
overlapping trace variables on a single cube. It is always possible to reduce these integrals 
to integrals involving non-overlapping trace variables using the orthogonality properties of 
the characters. For example, consider the expectation value, on a single cube, of a twice 
covered bent rectangle, with respect to the one-plaquette trial state of Eq. (4.22), 



X J dP,dPeXr{P5)Xr{P6) ^P^Pe) + XuiP^Pe)] e'^^(^5+^«+^^^6^). (8.12) 

To proceed with this integral we need to perform character expansions. The orthogonality 
of the characters can then be used to write characters over two plaquettes in terms of one- 
plaqucttc characters. To demonstrate how this is done we consider the character expansion 
of one of the integrals in Eq. (8.12), 

J dPxr{P)X2{PP')e''^'^''^'"^ = Y^^r'J dPxAP)X2{PP'). (8.13) 

r' 

Here the c^' are given by 

cr' = j dPxAP)Xr{P)e''^^''^'"\ (8.14) 
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Making use of the orthogonality of characters, given by Eq. (6.16), we obtain 



J dPxr{P)X2{PP')e''^'^''^'"^ = ^X2{P)X2{PP'). (8.15) 

Making use of this and proceeding similarly for the analogous integral involving Xii{.PbP&)i 
we can reduce Eq. (8.12) to 



(00 



do 



/ rfQx2(Q)x.(Q)e^^('2+^')]' + ^ [/ dQxii(Q)x.(Q)e^^(«+'2') 



•. (8.16) 



The plaquette matrix element provides a more straightforward example; 

(<^o| Q I'/'o) [/ dPxr{P)e'^^^''^'"^f J dQTVQxr(Q)e'=^'-(«+«^). (8.17) 

To proceed further we need expressions for the integrals over characters that appear in 
each matrix element of interest. We define these character integrals generically by 



Cr,r,...r„ (c) = J dPXr, {P)Xr, {P) ' ' ' Xr„ (P)e'=^(^+^^) . 



(8.18) 



All SU{N) integrals encountered in the calculation of glueball masses in 3+1 dimensions 
can be expressed in terms of character integrals. It is possible to calculate them generally, 
however three points need to be considered. Firstly, if we are to use the machinery of 
Chapter 5, we need to know how to express a given character in terms of trace variables. 
Secondly, as the dimension of the gauge group increases the number of independent trace 
variables increases rapidly. In order to contain the number of integrals required for a 
given order of approximation, it is necessary to express all high order trace variables in 
terms of a basis of lower order ones. Finally, the collection of SU(A'^) integrals presented 
in Chapter 5 must be extended. We now address each of these points in turn. 



The problem of expressing the characters of SU(7V) in terms of trace variables was 
solved by Bars [97] who showed that for SU(A'"), 

Xn{U)= Y: ^[E'k,-n]]J—^{TrWp. (8.19) 
ki,...,k„ \i=i / j=i 
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General characters can then be expressed in terms of Xn{U) using standard techniques 
from group theory [92] , 



lrir2---rN- 



,([/) =det [xr-,+i-,(f/)]i<,,.<jv_i- (8-20) 



For a given character, Eq. (8.20) produces a multinomial of traces of different powers of 
U. However, not all of these trace variables are independent. The so called Mandelstam 
constraints define the relationship between dependent trace variables. For SU(2) and 
SU(3) we have 

TrU = Tri/t V?7 G SU(2) and 

^ ' (8.21) 

Ti{U^) = (Tri7)2 - 2Tr(i7t) V?7 G SU(3). 

Similar relations can be obtained to express all trace variables in terms of TrU for SU(2), 
and for SU(3) in terms of TrU and TrU^ . To reduce the calculation of character integrals 
to a manageable size, we require expressions for high power trace variables in terms of a 
minimal set of lower order trace variables. To do this we proceed as follows. We start with 

an alternative expression of det U = 1, satisfied for all SU(A'') matrices U, 

^ili2---iN^iljl^i2j2 ' ' ' ^iNjN ~ ^jlj2---jN- (8.22) 

Here the colour indices of the group elements have been made explicit and all repeated 
indices are summed over. Multiplying Levi-Civita symbols produces sums of products 
of delta functions, the precise form of which depends on how many pairs of indices are 
contracted. A standard result from differential geometry will be useful here, 

£ai...a„£6i...6„ = det (<5ai6, ) i<^_^<„ ■ (8-23) 

Multiplying both sides of Eq. (8.22) by Uj^j^^^ and contracting over repeated indices 
gives 

^ili2---iN-ljN+l^hjl^i2j2 ' ' ' = ^jlj2---i]V^jjvijV+l 

^lll2---lN-liN^hi2--iN^hjl^i2j2 ' ' ' ^iN-dN-l ~ ^lll2---lN-liN^jlj2—jN^jj^iff- (8.24) 



In the last line we have renamed dummy indices and multiplied both sides by a Levi-Civita 
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symbol. Making use of Eq. (8.23) in Eq. (8.24) produces a product of delta functions on 
each side of the equation. By introducing appropriate combinations of delta functions and 
contracting over repeated indices we can construct trace identities as we please. We find 
that identities useful in our context are produced by the introduction of delta functions 
and an additional matrix, A, which is not necessarily an element of SU(Ar), into Eq. (8.24) 
as follows: 

N-l 

^kl---kN-liN^h---iN^iljl^jlkl ^kmjrrJ^imjm 

N-l 

= ^ki---kN-iiN^ji---jN^jikiUjj^i^ Y\_ ^kmjm- (8.25) 

m=2 

Here we have again renamed dummy indices. Some examples of SU(A'^) identities obtained 
from Eq. (8.25) valid for all N x N matrices A are as follows: 

-TrUTi{UA) + Ti{U^A) = -TrATiU^ + Tr(C/U) V f/ e SU(3) 
-{TrUfTr{UA) + Tr(C/A)Tr(C/^) + 2Tj:UTi{U^A) - 2Tr{U^A) 

= -2TiATiU^ + 2Tr{U^ A) V i7 G SU(4).(8.26) 

We require formulae for Tr(C7") in terms of lower order trace variables for SU(A^). Such for- 
mulae, known commonly as Mandelstam constraints, are obtained by setting A = 1/^-^+^^ 
in Eq. (8.25). In Appendix B we list the Mandelstam constraints obtained in this way up 
to SU(8). We see that for SU(A/^) it is possible to express all characters in terms of A'" — 1 
trace variables. In what follows we will choose to express the general SU(A/") characters in 
terms of the set of trace variables, {TrU^ , TrU, TrU'^, . . . , TrU'^~^}. It may well be possible 
to reduce the size of this set. This would indeed improve the efficiency of our technique. 
However, no effort has been made to do this at this stage. 

The procedure for calculating the general character integral, Cr-i^...r„ic), is then as 
follows. We first express the characters Xri,---,Xrn terms of trace variables using 
Eqs. (8.19) and (8.20). We then simplify these expressions using the Mandelstam con- 
straints of Appendix B. The final step is to perform the integrals over trace variables, 
which is an increasingly non-trivial task as one increases the dimension of the gauge group. 
For instance, the general character integral for SU(A') involves integrals over powers of the 
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trace variables TrJ7^, TrU, TrlJ^, . . . TiU^ ^. To proceed we need the following integral, 

^o^''o!Hc)= / d?7(TrC/^ifi(TrJ7'^2)P2.. /XrC/5'=)P'=e^^'"(^+^^). (8.27) 
Jsmm 



For SU(iV), the cases of interest to us here are described by, j''^^-'^^ j^_2(c). To calculate 
this integral we need to extend the work of Chapter 5 and consider the generating function, 

G,,..,,(c,7i, . . . ,7.) = / dC/e^^(^+^^)+2L.7.TV(f/.). (8.28) 

JSV{N) 

Following the procedure of Chapter 5 we obtain, 

oo 

G5i...q^(c,7i,...,7fc) = J2 det[A;+j_i,gi,...,qfc(c,7i,...,7fc)]^<.,^.<^, (8.29) 



l=—oo 

with, 

Am,gi,...,gfc(c,7l,---,7fc) = ^ ~i ^Ini+siq,+-+Skqki^c). (8.30) 

•^"^ Si ! • • • Shi 

Sl,S2, — ,Sk=0 

We then obtain an expression for Tq^.'.'.q^ {c) by differentiating Gq^,,,q^ appropriately; 

(8.31) 



QPl-\ hPfc 

'^qi.'.'.qk ('^) ^ Q^l . . . QryPk Gqi -qk (c, 7l ) • • • ) 7fe) 



7l=--=7fc=0 



8.5.2 The Variational Ground State 



In this section we fix the variational ground state following the usual procedure of minimis- 
ing the unimproved vacuum energy density given by Eq. (6.2) with k = 0. This equation 
is independent of the number of dimensions. The dimensionality of the lattice arises at 
the stage of calculating the plaquette expectation value. 

Making use of Eq. (8.17), the variational parameter can be fixed as a function of 
(3 for the one-cube lattice. In practice, the character sum in Eq. (8.17) and the infi- 
nite ^-sum in Eq. (8.29) need to be truncated. We truncate the infinite /-sum at Umax 
and instead of summing over all SU(A'') characters, we sum over only those characters, 
r = (ri,r2, . . . ,rjv-i), with Vmax > ''i > ^2 > • • • > tn-i- With this truncation scheme, 
memory constraints restrict calculations of the variational SU(A'') ground state to N <7, 
when working with rmax = 2 and Imax = 2 on a desktop computer. 
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Figure 8.1: The one cube SU(iV) variational parameter as a function of (3 for various N 
showing the dependence on the character sum truncation. 



The dependence of the variational parameter on TjYidx with Ifnax — 2 is shown for 
various gauge groups in Fig. 8.1. As 

fffiax incrBasGS the variational parameter appears 
to converge for each N considered. Moreover the convergence appears to improve as the 
dimension of the gauge group is increased. With the exception of SU(3), the Vmax = 1 
and 2 results are indistinguishable on the range < ^ < 0.7, where ^ = l/{Ng'^). As this 
is the range of interest to us later in the chapter we restrict further calculations to Vmax = 1- 

The SU(A'^) variational parameters on the one-cube lattice with r^ax = 1 and Imax = 2 
are shown for various N in Fig. 8.2. The results do not differ greatly. The corresponding 
variational energy densities are shown in Fig. 8.3. 



8.5.3 Expressions for the Glueball Mass 

We follow precisely the method described in Section 6.5 for the calculation of 3+1 dimen- 
sional glueball masses on a single cube. Here, however, we choose a different basis of states 
to minimise over. Instead of rectangular loops, we use states which fit on a single cube. 
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Figure 8.2: The one cube SU(A^) variational parameter in units of A/" as a function of 
l/{Ng'^) for various N. 



We start with a basis of two states, the plaquettes and the bent rectangles, 

B = {|1),|2)}, 



with 



and 



i) = J2i^iix)-{F,ix))]\ct>o) 



Fi{x) 
F2{x) 



(8.32) 



(8.33) 



(8.34) 



In 3+1 dimensions, the "+" sign corresponds to the O"*""*" state and the "— " sign corresponds 
to the 1"^~ state [128]. In order to calculate the glueball masses, we need expressions for 
the matrix elements N^, and Df^, of Eqs. (6.6) and (6.7) respectively. Here the superscript, 
C, denotes the charge conjugation eigenvalue, C = ±1, of the state in question. Making 
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1 2 3 4 5 

1 

Figure 8.3: The one cube SU(A^) energy density in units of iV as a function of for 
various N. 



use of Eq. (6.17) and carefully counting the number of possible overlaps between different 
loops we arrive at 




((Q)) 



(8.35) 



Here we have introduced the notation 



((O)) = (^o|0|^o). 



(8.36) 
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The combinatorics which lead to the coefficients in the matrix elements in Eq. (8.35) are 
a result of coTinting the possible overlaps within a single cube. The remaining matrix 
elements N22 and D22 can be calculated similarly, resulting in 




Having calculated N^^, and D^l, , we follow the minimisation procedure described in Chap- 
ter 6 to arrive at the glueball mass, AM^'-^. 

8.5.4 Results 

In this section we present calculations of the (symmetric) and 1^ (antisymmetric) 
glueball masses on the one-cube lattice for SU(iV) with 4 < A*" < 7. We first define the 
rescaled glueball mass, /U^*", corresponding to AM^'-^ as follows: 

= log(aAM-r"'-) - §i log (^1^) + (8.38) 

Here, as usual, P and C denote the parity and charge conjugation eigenvalues of the state 
in question. Prom the discussion Section 2.6 asymptotic scaling of a glueball mass is ob- 
served if the corresponding rescaled glueball mass becomes constant for some range of 
couplings. 

The results for the rescaled symmetric glueball mass are shown in Fig. 8.4. Fig. 8.4(a) 
shows the rescaled glueball mass calculated with only plaquettes in the minimisation basis. 
Fig. 8.4(b) shows the same quantity calculated with the minimisation basis of Eq. (8.32). 
The aim of this exploratory study is to observe whether or not a move toward scaling is 
apparent as the number of states in the minimisation basis is increased. This is clearly not 
the case for the symmetric glueball mass on the range of couplings explored here. There 
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Figure 8.4: The one cube SU(A^) 0"'""'" rescaled glueball mass as a function of l/{Ng^) for 
various N obtained with different minimisation bases. 

is no sign of jU"'""'" becoming constant on any range of couplings within < ^ < 0.55. 



The prospects are marginally better for the antisymmetric case. The results for the 
rescaled antisymmetric glueball mass are displayed in Fig. 8.5. A move toward scaling is 
observed, most clearly for the N = 7 case, near ^ = 0.38. This matches the region of 
couplings for which Chin, Long and Robson observed scaling of the 0++ glueball mass, for 
N = 5 and 6, on a 6^ lattice using only plaquettes in their minimisation basis [138]. 

8.5.5 String Tension 

If the exploratory study presented here was to be extended to larger lattices, at some stage 
it would be useful to compute the string tension, a. It is common for calculations in La- 
grangian LGT to express results for masses in units of s/a. This allows masses calculated 
on the lattice to be expressed in MeV since the string tension can be calculated in MeV 
from the decay of heavy quarkonia for example. In the Hamiltonian formulation, precise 
calculations of the string tension have only been performed in the strong coupling regime. 
Variational estimates, at least for SU(2) [1], have not exhibited asymptotic scaling when 
making use of the one-plaquette ground state of Eq. (4.22). It is possible that the situa- 
tion may be improved for higher dimensional gauge groups but this has not yet been tested. 
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Figure 8.5: The one cube SU(A'^) l"*" rescaled glueball mass as a function of l/{Ng^) for 
various N obtained with different minimisation bases. 



In this section wc calculate the symmetric and antisymmetric glueball masses in units 
of i/o". For the string tension, since reliable variational results are not available, we use 
the strong coupling expansions of Kogut and Shigemitsu [143]. These are available for 
SU(5) and SU(6), as well as SU(2) and SU(3) which we do not consider here. 



Since the square root of the string tension has units of mass, the ratio of a mass to 
the string tension is constant in a scaling region. Again, with the crude model presented 
here, wc do not expect to observe scaling. We seek only an indication of an approach to 
scaling. Such behaviour would warrant further study. 



The results for the symmetric states are shown in Fig. 8.6. Since we perform calcula- 
tions with two states in the minimisation basis, the two lowest mass states are accessible. 
In Fig. 8.6(a) the results for the lowest mass state are shown. Calculations with different 
minimisation bases are shown, with the "SU(A^)-plaquettes" label indicating that only 
plaquettc states are used. For each case the glueball mass has a local minimum in the 
range 0.33 < ^ < 0.37. The various results do not differ greatly. There is no sign of 
improved scaling behaviour when either bent rectangles are included or TV is increased. 
The effect of including bent rectangles is to lower the local minimum. 
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The second lowest glueball mass is shown in Fig. 8.6(b). Again, no improvement in 
scahng behaviour is seen as N is increased. In Fig. 8.6 the horizontal lines indicate the 3+1 
dimensional SU(5) calculations (and error bars) of Lucini and Teper [2]. Interestingly, both 
masses calculated here with a simplistic model, have minima that lie within the error bars. 

We now move on to the antisymmetric states. The results for the antisymmetric glue- 
ball masses in units of ^/a arc shown in Fig. 8.7. The results arc more promising than 
the symmetric case. We see an improved approach to scaling when bent rectangles are 
included in the minimisation basis for the lowest glueball mass. For this case the glueball 
mass shows promising signs of becoming constant in the ranges 0.32 < ^ < 0.39 for SU(6) 
and 0.3 < ^ < 0.37 for SU(5). Only marginal improvement in scaling behaviour is evident 
as N is increased from 5 to 6. For each of these cases no data is available for comparison 
to our knowledge. Such promising signs are not apparent in the second lowest glueball 
mass in this sector, as seen in Fig. 8.7(b). 

It would be interesting to perform analogous calculations for SU(7), for which the most 
promising results were displayed in Section 8.5.4. However strong coupling expansions 
for the SU(7) string tension are not available in Ref. [143] and have not been published 
elsewhere to our knowledge. 

8.6 Future Work 

In this chapter we have studied the variational 0++ and l+~ glueball masses on a single 
cube in 3+1 dimensions. Our intention was to determine the viability of using analytic 
Hamiltonian methods in 3+1 dimensional glueball mass calculations. With such a crude 
small volume approximation and a minimisation basis containing only two states, the 
observation of asymptotic scaling was not expected. Promising signs of an approach to 
asymptotic scaling were displayed by the 1"' glueball mass, as N was increased, at cou- 
plings close to the scaling window observed by a comparable, but larger volume, study 
by Chin, Long and Robson [138]. Interesting results were also observed for calculations 
of the glueball masses in units of the string tension. With no variational results for the 
string tension available, the strong coupling results of Kogut and Shigemitsu [143] were 
used. Interestingly, the mass of both SU(5) 0++ glueball states calculated had minima 
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(a) Lowest energy state (b) Second lowest energy state 



Figure 8.6: The one cube SU(5) and SU(6) 0"'""'" glueball masses in units of as functions 
of l/{Ng'^). The horizontal hnes indicate the result and error bars of the SU(5) 0++ 
calculation of Lucini and Teper [2] . 



which were consistent with the Lagrangian calculation of Lucini and Teper [2]. Better 
scaling behaviour was exhibited by the SU(5) and SU(6) 1"' states, although no alterna- 
tive calculations are available for comparison to our knowledge. For the lowest mass 1"*"" 
states calculated, the scaling behaviour was improved by increasing the number of states 
in the minimisation basis and the dimension of the gauge group. The promising results 
observed in this chapter warrant further study on larger lattices, with additional states in 
the minimisation basis. 



To extend this calculation to larger lattices and minimisation bases, two challenges will 
be faced. Firstly, the correct implementation of the Bianchi identity, within our plaquette 
based analytic approach, forces the number of integration variables to grow quickly with 
the volume. New approaches to handling the general character integrals may need to be 
developed in order to avoid the inevitable memory restrictions. Secondly, great care will 
need to be taken in the counting of overlaps between different states. This process could, 
in principle, be automated using techniques from graph theory and symbolic programming. 

It would also be of interest to extend the calculation presented here to larger N. To 
do this a more efficient technique for handling the general character integrals would need 



8.6 Future Work 



165 



4.5 



4- 





1 

SU(6) 


1 1 






-— SU 5 








oU(D 


-plaquettes 






su(5: 


-plaquettes 










/ 












1 


1 1 





0.2 



0.25 



0.3 

1 



0.35 



(a) Lowest energy state 



0.4 




(b) Second lowest energy state 



Figure 8.7: The one cube SU(5) and SU(6) 1+ glueball masses in units of ^/a as functions 
of l/iNg-"). 



to be developed to minimise the memory demands of the calculation. The following is 
just one possibility. It is likely, that by expressing the products of characters appearing in 
the character integrals as sums of characters, one could reduce the memory requirements 

of the calculation drastically. The general character integrals could then be expressed in 
terms of integrals of the form. 



Gric) 



SV{N) 



(8.39) 



where r = (ri,r2, . . . rjv-i) labels a representation of SU(A^). Having completed the cal- 
culations presented in this chapter, it was discovered that Gr{c) can be handled using the 
techniques of Chapter 5, with the result. 



Gr{c) = det[/^.+i+j_i(2c)]^<._ 



,j<N ■ 



(8.40) 



l=—oo 



The final stage of this improvement would be to find a convenient way to express the 
products of characters appearing in the calculation as a sum of characters. This would be 
possible with the symbolic manipulation of Young tableaux. 
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Chapter 9 



Conclusion 



This thesis has been concerned with the study of two topics: improvement and the use 
of analytic techniques in Hamiltonian LGT. In this chapter we summarise our key results 
and discuss the directions that further work should take. 

Chapters 3 and 4 were concerned primarily with the issue of improvement. In Chap- 
ter 3 a technique allowing the straightforward construction of lattice Hamiltonians was 
developed. Using this direct approach a number of improved Hamiltonians were derived. 
The key advantage of our approach is that it allows the straightforward construction of 
improved kinetic Hamiltonians with a finite number of terms, a task that has proved 
nontrivial in the past. In Chapter 4 two simple tests of improvement, a strong coupling 
calculation of the U(l) Callan-Symanzik /3 function in 3+1 dimensions and a mean field 
calculation of the SU(2) lattice specific heat in 2+1 dimensions, were presented. For each 
example considered, the tadpole improved Hamiltonians (when they were used) achieved 
the greatest level of improvement followed by the classically improved Hamiltonians and 
then the Kogut-Susskind Hamiltonians. 

In Chapter 5, techniques which allow the analytic calculation of SU(A^) matrix ele- 
ments in 2+1 dimensions were developed. This chapter generalised the methods that have 
been employed in variational SU(2) LGT calculations for almost 20 years, to the general 
case of SU(Ar). A close connection with the fields of random matrices and combinatorics 
was also discussed. 
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Conclusion 



In Chapter 6, the analytic techniques of Chapter 5 were apphed in the calculation of 
N = 2,3,4, and 5 glueball masses in 2+1 dimensions. Kogut-Susskind, improved and 
tadpole improved Hamiltonians, all of which were derived in Chapter 3, were used and 
their results compared. For each case considered, the lowest mass results obtained were 
considerably higher than the lowest glueball masses obtained by comparable Hamiltonian 
and Lagrangian calculations. However, the lowest mass states obtained in the 0++ and 
sectors were close in mass to the respective excited states, O"^"^* and **, of Teper 
and Lucini [51, 106]. In the 0"'"''' sector this agreement improved as N was increased, and 
for each N the agreement was closer for the improved calculations. This situation was 
replicated in the sector if the SU(4) results were excluded. For SU(4), the lowest 
unimproved glueball mass was closer to the mass of Teper and Lucini's ** state than 
the improved results. The problem of enumerating states can be attributed to the fact 
that the calculations presented in Chapter 6 used a basis of rectangular states. The com- 
peting calculations incorporate additional nonrectangular states. Without the inclusion 
of small area, nonrectangular states in the calculation, it is possible that some low mass 
states are inaccessible. Further calculations, which include nonrectangular states in the 
minimisation basis, are needed to clarify this issue. Including such additional states will 
also allow spin 1 and P = —C states to be studied, and permit a more complete study of 
the spin and 2 glueball mass spectra. 

In Chapter 7 we presented the first steps of a study of the SU(A^) glueball mass spec- 
trum in the large N limit using the techniques of Chapter 6. Our motivation was the 
Maldacena conjecture which, in its simplest form, predicts a precise correspondence be- 
tween certain string theories and gauge theories in the large A'' limit. Using analytic 
Hamiltonian techniques, calculations with N as large as 25 were possible on a desktop 
computer, in contrast to comparable Lagrangian calculations which have reached N = 6. 
With such large N glueball masses accessible, reliable extrapolations to the N ^ oo were 
possible. The results provided evidence for leading order finite N corrections to the 

glueball masses, verifying a specific prediction of large N gauge theory. 

In Chapter 8 we moved on to the physically interesting case of 3-1-1 dimensions. There 
we studied the viability of adapting the analytic techniques used with success in Chapters 6 
and 7 to 3-1-1 dimensions. A simple variational calculation of 0++ and l"* glueball masses 
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was carried out on a single cube, with only two states in the minimisation basis. With such 
a simplistic approximation asymptotic scaling was not expected. However, it was hoped 
that with the use of analytic Hamiltonian techniques, large values of N could be reached. 
This was not the case, with memory constraints restricting calculations to A?^ < 7 on a one 
cube lattice. Our motivation for attempting this calculation was a variational Hamiltonian 
study by Chin, Long and Robson [138] which demonstrated asymptotic scaling of the 0+"*" 
glueball mass on a 6^ lattice for A?^ = 5 and 6, with improved scaling observed as N was in- 
creased. Promising signs of an approach to asymptotic scaling were observed in Chapter 8 
for the 1"' state, in the same region of couplings as scaling was seen by Chin, Long and 
Robson. No such signs were observed for the O"*""*" state. The extension to larger lattices 
thus seems warranted. Such extensions will be hampered by memory constraints unless 
new techniques for handling the general character integrals are developed. Another im- 
portant improvement, which does not pose a serious problem, will be the inclusion of more 
states in the minimisation basis. The applicability of analytic variational techniques in 
3-1-1 dimensions will thus be determined by how quickly the infinite volume limit is reached. 

While we have derived various improved Hamiltonians and put them to use in glueball 
mass calculations in 2-1-1 dimensions, their use in 3-1-1 dimensional non-abelian LOT is 
untested. The next stages of a 3-1-1 dimensional calculation, within the variational analytic 
approach presented here, would incorporate improved Hamiltonians and move beyond a 
single cube. It is possible that in such a calculation the use of an improved Hamiltonian 
could mean the difference between observing and not observing scaling on the range of 
couplings available. 

The results of this thesis demonstrate that the analytic techniques of Hamiltonian LGT 
can compete with the traditional Monte Carlo techniques of Lagrangian LGT, at least in 
the calculation of glueball masses in 2-1-1 dimensions. The use of the analytic techniques 
presented here are by no means restricted to variational calculations. It is possible that 
they will find use in coupled cluster calculations, or in the moment calculations of the 
plaquette expansion method [26]. The move to 3-1-1 dimensions is formidable but the 
calculations of Chapter 8 suggest that analytic techniques may by feasible, at least for 
large A'^, provided the infinite volume limit is reached quickly. 



170 Conclusion 



Bibliography 



[1] H. Arisue, M. Kato, and T. Fujiwara, Prog. Theor. Phys. 70, 229 (1983). 

[2] B. Lucini and M. Teper, JHEP 06, 050 (2001), hep-lat/0103027. 

[3] D. J. Gross and F. Wilczek, Phys. Rev. Lett. 30, 1343 (1973). 

[4] H. D. Politzer, Phys. Rev. Lett. 30, 1346 (1973). 

[5] V. N. Gribov and L. N. Lipatov, Yad. Fiz. 15, 781 (1972). 

[6] K. Hagiwara, K. Hikasa, K. Nakamura, M. Tanabashi, M. Aguilar-Benitez, C. Am- 
sler, R. Barnett, P. Burchat, C. Carone, C. Caso, et al., Phys. Rev. D66, 010001 
(2002), URL http://pdg.lbl.gov. 

[7] A. Chodos, R. L. JafFe, K. Johnson, C. B. Thorn, and V. F. Weisskopf, Phys. Rev. 
D9, 3471 (1974). 

[8] J. M. Cornwall and A. Soni, Phys. Lett. B120, 431 (1983). 

[9] N. Isgur and J. Paton, Phys. Rev. D31, 2910 (1985). 

[10] R. W. Johnson and M. J. Teper, Phys. Rev. D66, 036006 (2002), hep-ph/0012287. 

[11] K. G. Wilson, Phys. Rev. DIO, 2445 (1974). 

[12] J. B. Kogut and L. Susskind, Phys. Rev. Dll, 395 (1975). 

[13] M. Beneke (2002), hep-lat/0201011. 

[14] S. M. Ryan, Nucl. Phys. Proc. Suppl. 106, 86 (2002), hep-lat/0111010. 

[15] R. Gupta (1997), hep-lat/9807028. 

[16] J. B. Kogut, Phys. Rept. 67, 67 (1980). 



172 



BIBLIOGRAPHY 



[17] G. S. Bali et al. (TXL), Phys. Rev. D62, 054503 (2000), hep-lat/0003012. 

[18] S. Aoki et al. (CP-PACS), Phys. Rev. Lett. 84, 238 (2000), hep-lat/9904012. 

[19] T. Yoshie (CP-PACS), Nucl. Phys. Proc. Suppl. 109 A, 25 (2002). 

[20] E. B. Gregory, S.-H. Guo, H. Kroger, and X.-Q. Luo, Phys. Rev. D62, 054508 (2000), 
hep-lat/9912054. 

[21] B. Bringoltz and B. Svetitsky (2002), hep-lat/0209005. 

[22] F. Coester, Nucl. Phys. 7, 421 (1958). 

[23] F. Coester and H. G. Kiimmel, Nucl. Phys. 17, 477 (1960). 

[24] N. E. Ligterink, N. R. Walet, and R. F. Bishop, Annals Phys. 284, 215 (2000), 
hep-lat/0001028. 

[25] B. H. J. McKellar, C. R. Leonard, and L. C. L. Hollenberg, Int. J. Mod. Phys. B14, 
2023 (2000). 

[26] L. C. L. Hollenberg, Phys. Rev. D47, 1640 (1993). 
[27] L. C. L. Hollenberg, Phys. Rev. D50, 6917 (1994). 

[28] J. A. L. Mcintosh and L. C. L. Hollenberg, Nucl. Phys. Proc. Suppl. 106, 257 (2002), 
hep-lat/0111045. 

[29] J. A. L. Mcintosh and L. C. L. Hollenberg, Phys. Lett. B538, 207 (2002), 
hep-lat/0111061. 

[30] H. Arisue, Prog. Theor. Phys. 84, 951 (1990). 

[31] H. Pritzsch and M. Gell-Mann, in Proceedings of the XVI International Conference 
on High Energy Physics, edited by J. D. Jackson and A. Roberts (1972), vol. 2, p. 
135. 

[32] H. Fritzsch, M. Gell-Mann, and H. Leutwyler, Phys. Lett. B47, 365 (1973). 
[33] G. S. Bah (2001), hep-ph/0110254. 

[34] C. Morningstar, AIP Conf. Proc. 619, 231 (2002), nucl-th/0110074. 



BIBLIOGRAPHY 



173 



[35] S. Godfrey and J. Napolitano, Rev. Mod. Phys. 71, 1411 (1999), hep-ph/9811410. 

[36] C. Amsler, Phys. Lett. B541, 22 (2002), hep-ph/0206104. 

[37] C. Amsler and F. E. Close, Phys. Rev. D53, 295 (1996), hep-ph/9507326. 

[38] W.-J. Lee and D. Weingarten, Phys. Rev. D61, 014015 (2000), hep-lat/9910008. 

[39] J. Sexton, A. Vaccarino, and D. Weingarten, Phys. Rev. Lett. 75, 4563 (1995), 
hep-lat/9510022. 

[40] M. Acciarri et al. (L3), Phys. Lett. B501, 173 (2001), hep-ex/0011037. 
[41] K. Benslama et al. (CLEO), Phys. Rev. D66, 077101 (2002), hep-ex/0204019. 
[42] K. Symanzik, Nucl. Phys. B226, 187 (1983). 
[43] K. Symanzik, Nucl. Phys. B226, 205 (1983). 

[44] F. Jegerlehner, R. D. Kenway, G. Martinelli, C. Michael, O. Pene, B. Petersson, 
R. Petronzio, C. T. Sachrajda, and K. Schilling (2000), ECFA-99-200, CERN-2000- 
002. 

[45] G. P. Lepage, in Perturbative and Nonperturbative Aspects of Quantum Field Theory 
: Proceedings of the 35. Internationale Universitdtswochen fur Kern- und Teilchen- 
physik, Schladming, Austria, March 2-9, 1996, edited by H. Latal and W. Schweiger 
(1996), p. 1. 

[46] X.-Q. Luo, S.-H. Guo, H. Kroger, and D. Schiitte, Phys. Rev. D59, 034503 (1999), 
hep-lat/9804029. 

[47] A. Wichmann, Ph.D. thesis, ITKP Bonn (2001). 

[48] S. A. Chin, C. Long, and D. Robson, Phys. Rev. D37, 3001 (1988). 

[49] C. Long, D. Robson, and S. A. Chin, Phys. Rev. D37, 3006 (1988). 

[50] B. Berg and A. Billoire, Nucl. Phys. B221, 109 (1983). 

[51] M. J. Teper, Phys. Rev. D59, 014512 (1999), hep-lat/9804008. 

[52] R. Johnson and M. Teper, Nucl. Phys. Proc. Suppl. 73, 267 (1999), 
hep-lat/9808012. 



174 



BIBLIOGRAPHY 



[53] R. W. Johnson, Phys. Rev. D66, 074502 (2002), hep-lat/0206005. 

[54] D. Q. Liu and J. M. Wu, Mod. Phys. Lett. A17, 1419 (2002), hep-lat/0105019. 

[55] W. E. Caswell, Phys. Rev. Lett. 33, 244 (1974). 

[56] D. R. T. Jones, Nucl. Phys. B75, 531 (1974). 

[57] T. van Ritbergen, J. A. M. Vermaseren, and S. A. Larin, Phys. Lett. B400, 379 
(1997), hep-ph/9701390. 

[58] R. F. Dashen and D. J. Gross, Phys. Rev. D23, 2340 (1981). 

[59] A. Hasenfratz and P. Hasenfratz, Nucl. Phys. B193, 210 (1981). 

[60] M. Liischer and P. Weisz, Nucl. Phys. B452, 234 (1995), hep-lat/9505011. 

[61] C. J. Hamer, Phys. Rev. D53, 7316 (1996). 

[62] M. J. Tepcr, in Confinement, Duality, and Nonperturbative Aspects of QCD, edited 
by P. van Baal (1998), hep-lat/9711011. 

[63] M. Creutz, Quarks, Gluons and Lattices (Cambridge University Press, 1983), (Cam- 
bridge Monographs On Mathematical Physics). 

[64] J. B. Kogut, Rev. Mod. Phys. 55, 775 (1983). 

[65] M. Liischer and P. Weisz, Phys. Lett. B158, 250 (1985). 

[66] X.-Y. Fang, S.-H. Guo, and J.-M. Liu, Mod. Phys. Lett. A15, 737 (2000). 

[67] G. P. Lepage and P. B. Mackenzie, Phys. Rev. D48, 2250 (1993), hep-lat/9209022. 

[68] M. G. Alford, W. Dimm, G. P. Lepage, G. Hockney, and P. B. Mackenzie, Nucl. 
Phys. Proc. Suppl. 42, 787 (1995), hep-lat/9412035. 

[69] M. G. Alford, W. Dimm, G. P. Lepage, G. Hockney, and P. B. Mackenzie, Phys. 
Lett. B361, 87 (1995), hep-lat/9507010. 

[70] D. Horn, M. Karliner, and M. Weinstein, Phys. Rev. D31, 2589 (1985). 

[71] G. G. Batrouni and M. B. Halpern, Phys. Rev. D30, 1782 (1984). 

[72] J. M. Maldacena, Adv. Theor. Math. Phys. 2, 231 (1998), hep-th/9711200. 



BIBLIOGRAPHY 



175 



[73] E. Witten, Adv. Theor. Math. Phys. 2, 253 (1998), hep-th/9802150. 

[74] I. Montvay and G. Miinster, Quantum Fields on a Lattice (Cambridge University 
Press, 1994), chap. 3.2.3, (Cambridge Monographs On Mathematical Physics). 

[75] K. E. Eriksson, N. Svartholm, and B. S. Skagerstam, J. Math. Phys. 22, 2276 (1981). 

[76] C. R. Leonard, Ph.D. thesis, The University of Melbourne (2001). 

[77] I. M. Gessel, J. Combin. Theory Ser. A 53, 257 (1990). 

[78] R. P. Stanley, Enumerative Combinatorics (Cambridge University Press, 1999), 
vol. 2, p. 453. 

[79] B. E. Sagan, The Symmetric Group (Springer- Verlag, 2001), chap. 2.11, 2nd ed. 

[80] J. Baik and E. M. Rains, Duke Math. J. 2, issue 1:1 (2001), math. CO/9905083. 

[81] E. M. Rains, Electron. J. Combin. 5(1), R12 (1998). 

[82] H. Widom, A. R. Its, and C. A. Tracy, Physica D152-153, 199 (2001). 

[83] A. Regev, Adv. in Math. 41, 115 (1981). 

[84] P. Diaconis and M. Shahshahani, J. Appl. Probab. 31 A, 49 (1994). 

[85] H. Widom and C. A. Tracy, Comm. Math. Phys. 207, 665 (1999). 

[86] M. Creutz, J. Math. Phys. 19, 2043 (1978). 

[87] R. De Pietri, Class. Quant. Grav. 14, 53 (1997), gr-qc/9605064. 

[88] C. Rovelh and L. Smolin, Phys. Rev. D52, 5743 (1995), gr-qc/9505006. 

[89] K. Ezawa, Phys. Rept. 286, 271 (1997), gr-qc/9601050. 

[90] J. B. Kogut, M. Snow, and M. Stone, Nucl. Phys. B200, 211 (1982). 

[91] R. Brower, P. Rossi, and C.-I. Tan, Nucl. Phys. B190, 699 (1981). 

[92] H. Weyl, The Classical Groups, Their Invariants and Representations (Princeton 
University Press, 1946). 

[93] I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products (Aca- 
demic Press, 1994). 



176 



BIBLIOGRAPHY 



[94] P. van Mocrbeke, in Random Matrices and Their Applications (Cambridge Univer- 
sity Press, 2001), MSRI-publication #40, math. CO/0010135. 

[95] N. J. A. Sloane, The On-line Encyclopedia of Integer Sequences, 
http : //www. research, att . coin/~nj as/ sequences/. 

[96] S. P. Tonkin, Nucl. Phys. B292, 573 (1987). 

[97] I. Bars, J. Math. Phys. 21, 2678 (1980). 

[98] J. Greensite, Phys. Lett. B191, 431 (1987). 

[99] J. P. Greensite, Nucl. Phys. B158, 469 (1979). 

[100] C. J. Hamer, J. Oitmaa, and W.-H. Zheng, Phys. Rev. D45, 4652 (1992). 

[101] C. H. Llewellyn Smith and N. J. Watson, Phys. Lett. B302, 463 (1993), 
hep-lat/9212025. 

[102] Q.-Z. Chen, S.-H. Guo, X.-Y. Fang, and W.-H. Zheng, Phys. Rev. D50, 3564 (1994). 

[103] X.-Q. Luo and Q.-Z. Chen, Mod. Phys. Lett. All, 2435 (1996), hep-ph/9604395. 

[104] S. Samuel, Phys. Rev. D55, 4189 (1997), hep-ph/9604405. 

[105] J.-M. Li, S. Guo, and X.-Q. Luo, Commun. Theor. Phys. 34, 301 (2000). 

[106] B. Lucini and M. Teper (2002), hep-lat/0206027. 

[107] G. 't Hooft, Nucl. Phys. B72, 461 (1974). 

[108] E. Witten, Nucl. Phys. B160, 57 (1979). 

[109] T. Eguchi and H. Kawai, Phys. Rev. Lett. 48, 1063 (1982). 

[110] C. Csaki, H. Ooguri, Y. Oz, and J. Terning, JHEP 01, 017 (1999), hep-th/9806021. 

[Ill] R. de Mello Koch, A. Jevicki, M. Mihailescu, and J. P. Nunes, Phys. Rev. D58, 
105009 (1998), hep-th/9806125. 

[112] M. Zyskin, Phys. Lett. B439, 373 (1998), hep-th/9806128. 

[113] R. C. Brower, S. D. Mathur, and C.-L Tan, Nucl. Phys. B574, 219 (2000), 
hep-th/9908196. 



BIBLIOGRAPHY 



177 



[114] R. C. Browcr, S. D. Mathur, and C.-I. Tan, Nucl. Phys. B587, 249 (2000), 
hep-th/0003115. 

[115] J. G. Russo, Nucl. Phys. B543, 183 (1999), hep-th/9808117. 

[116] C. Csaki, J. Russo, K. Sfetsos, and J. Terning, Phys. Rev. D60, 044001 (1999), 
hep-th/9902067. 

[117] S. Dalley and B. van de Sande, Phys. Rev. D63, 076004 (2001), hep-lat/0010082. 

[118] D. Karabah and V. P. Nair, Phys. Lett. B379, 141 (1996), hep-th/9602155. 

[119] D. Karabah, C.-J. Kim, and V. P. Nair, Nucl. Phys. B524, 661 (1998), 
hep-th/9705087. 

[120] D. Karabah, C.-J. Kim, and V. P. Nair, Phys. Lett. B434, 103 (1998), 
hep-th/9804132. 

[121] V. P. Nair, Nucl. Phys. Proc. Suppl. 108, 194 (2002). 

[122] C. Krattenhaler, Seminaire Lotharingien Combin. 42, B42q (1999), The Andrews 
Festschrift. 

[123] X.-Y. Fang, P. Hui, Q.-Z. Chen, and D. Schutte, Phys. Rev. D65, 114505 (2002). 

[124] Q.-Z. Chen, X.-Q. Luo, S.-H. Quo, and X.-Y. Fang, Phys. Lett. B348, 560 (1995), 
hep-ph/9502235. 

[125] G. G. Batrouni, Nucl. Phys. B208, 467 (1982). 

[126] G. G. Batrouni, Ph.D. thesis, Lawrence Berkeley Laboratory, University of California 
(1983), lBL-15993. 

[127] J. B. Kogut, R. B. Pearson, and J. Shigemitsu, Phys. Rev. Lett. 43, 484 (1979). 

[128] C. J. Earner, Phys. Lett. B224, 339 (1989). 

[129] D. Horn and M. Weinstein, Phys. Rev. D30, 1256 (1984). 

[130] C. P. van den Doel and D. Horn, Phys. Rev. D33, 3011 (1986). 

[131] C. P. van den Doel and D. Horn, Phys. Rev. D35, 2824 (1987). 



178 



BIBLIOGRAPHY 



[132] D. Horn and G. Lana, Phys. Rev. D44, 2864 (1991). 

[133] D. Horn and D. Schreiber, Phys. Rev. D47, 2081 (1993). 

[134] R. F. Bishop and H. G. Kiimmel, Physics Today 40(3), 52 (1987). 

[135] D. Schiitte, W.-H. Zheng, and C. J. Hamer, Phys. Rev. D55, 2974 (1997), 
hep-lat/9603026. 

[136] S.-H. Guo, Q.-Z. Chen, and L. Li, Phys. Rev. D49, 507 (1994). 

[137] L. Hu, X.-Q. Luo, Q.-Z. Chen, X.-Y. Fang, and S.-H. Guo, Commun. Theor. Phys. 
28, 327 (1997), hep-ph/9609435. 

[138] S. A. Chin, C. Long, and D. Robson, Phys. Rev. Lett. 57, 2779 (1986). 

[139] M. B. Halpern, Phys. Rev. D19, 517 (1979). 

[140] G. G. Batrouni, Nucl. Phys. B208, 12 (1982). 

[141] L. D. Fadeev and A. A. Slavnov, Gauge Fields: Introduction to Quantum Theory 
(The Benjamin/Cummings Publishing Company, Inc., 1980), chap. 3.2. 

[142] R. Jackiw, in Lectures on QCD, edited by F. Lenz, H. Griesshammer, and D. Stoll 
(Springer- Varlag, 1997), pp. 90-127. 

[143] J. B. Kogut and J. Shigemitsu, Phys. Rev. Lett. 45, 410 (1980), erratum- 
ibid.45:1217,1980. 



Appendix A 



Commutation Relations 



In improved Hamiltonian LGT calculations, one encounters matrix elements of the form: 
(</.o|^Tr[f,(a;)fi(a;)]|(/)o) and ((/-ol ^ TV [£:,(a;)C/,(a;)£:i(a; + ai)C//(a;)] |0o). (A.l) 

The first of these is easily handled. One simply writes the electric field operators as 
Si{x) = £^{x)X" and uses the normalisation condition, Eq. (3.5), to give 

(<^o|^TV[f,(a;)f,(x)]|0o) = ^ ^^(./'olff (a;)fr(a;)|</>o). (A.2) 

Let the trial state, \(t)o), have the form l^o) = e'^|0), where 5 is a function of link variables 
and = S. Then 

j2{M£nx)£n^Mo) = [£"(^)^^'] [£n^),e'] io). (a.s) 

Making use of the Baker-Hausdorff formula, Eq. (3.19), we can derive the following result, 

e^£t{x)e-^ = £t{x)-[£t{x),S] + \[[Sr{x),SlS] + ... 

= £rix)-[£r{x),S]. (A.4) 

The last line follows from the fact that 5 is a function of link variables and the commutation 
relation of Eq. (8.4). Some rearrangement of Eq. (A.4) leads to the useful result, 

[£t{x), e^] = [£t{x), S] = [£t{x), S] , (A.5) 
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Commutation Relations 



which we apply in Eq. (A. 3) to obtain, 

^(</>o|^r(x)^r(a;)|</>o) = -Y,{M[£nx),S][Snx),S]\<Po) 

= lY.(^o\[£nx)A£nx),s]]\cpo). (A.6) 

x,i 

In the last line we have used integration by parts [1]. 



The second of the matrix elements in Eq. (A.l) is more difficult to handle. For this 
case it is convenient to expand Si{x + ai) in the basis \f{x + ai): 



£i{x + ai) = Af (a; + ai)£^{x + ai), 



(A.7) 



where, 



Xf{x + ai) = ul{x)\'^Ui{x). 



(A.8) 



We now derive the commutation relations between S^{x) and Uj{y). We start with the 
analogous relation to Eq. (3.5), 



Tr 



Xt{x)Xf{x) =-S^p 



(A.9) 



which follows trivially from Eq. (A.8). This relation allows us to invert Eq. (A.7) and 
write in terms of f " 



2Tr 
2Tr 



£i{x)Xf{x) 



Xt{x)\^ £^{x). 



(A.10) 



We can form a similar relation between A" and A" as follows. Since Xf{x) G S'U{N), we 
can expand in the Gell-Mann basis as follows, 



Af(a;) = cr(a.)A^ 



(A.11) 



Here cf^{x) are constants which we now determine. Multiplying Eq. (A. 11) throughout 
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by A'^, tracing and making use of Eq. (3.5) gives 



cf{x) = 2TV te(a;)A^ 



(A.12) 



From Eq. (A. 10) we see that £^{x) = c^^ {x)£^ (x). The commutation relations foUow 
immediately, 



St{x),U,{y)\ = cfix)[s^ix),Ujiy)]=cf{x)\f^d,j5^yU,{x) 

= 6ijS^y~Xfix)Ui{x). (A.13) 



These commutation relations enable the simplification the second of the matrix elements 
in Eq. (A.l) using the basis defined in Eq. (A. 7). Making use of Eqs. (A.IO) and (A.13) 
leads to the following simplification: 



^o\Y,T^T[£iix)Uiix)£i{x+ai)uj{x)\\(l)o) = ^ ^(0o| [ff (oj), [^f (x+ai), |<^o).(A.14) 



In order to calculate such matrix elements the following results are useful 



[£nx),mx)}AB 



- at) 



£^{x + ai),Ui{x)Ui{x + 
£f{x + ai), {Uiix^- ai))cD 



AT^-l 



Ui{x)Ui{x + ai), (A. 15) 



^SBc{Uiix)}AB'{Uiix + ai)}B'D - ^{Ui{x)}AB{Ui{x + ai)}cD, 

£t{x), {U^i{x)}AB\ [£t{x + ai), {Ui{x + ai)}cD 



(A.16) 



^{uj{x)}cBmx + ai)}AD + ^{U;{x)}ABmx + ai)}cD, (A.17) 

£^{x + ai),{U^i{x + ai)}cD = 



[£nx),mx)}AB 

-l{Ui{x)}cB{ulix + ai)}AD + -^{Uiix)}AB{ul{x + ai)}cD. 



(A.18) 



Here {X}ab denote the colour indices of the matrix X and an implicit sum over repeated 
colour indices is understood. These results follow simply from the commutation relations 
of Eqs. (8.4) and (A.13), and the following SU(7V) formula: 



^AB^CD = 2^AdSbC - -^SabScD- 



(A.19) 
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Appendix B 



Mandelstam Constraints 



TrU"- = TrU"--^ - TrU^TrU""'^ + Trf/Trf/""^ yU G SU(3) (B.l) 
W = -TrC/"-^ + Tr?7^TrC/""3 - ^{TiUfTrU"'-^ + ^TrU'^TrU"--^ 

+Tri7TV?7"-^ V?7 G SU(4) (B.2) 

3 2^ ^ 2 

+TrC/TV?7"-^ V?7 G SU(5) (B.3) 

TVC/" = -TVf/"-6 + TrJ/^TrC/^-s - — (TrC/)^TVC/"-^ + -(Trf/)^ TrC/^TVC/""^ 

24 4 

8^ 3 4 6^ ^ 

--Tr[/TrC/2Tr?7"-''^ + -TrC/^TrC/"-^ - -Tri/^Trt/"-^ + -TrC/^TrC/"-^ 
2 3 2 2 

+Trf/TVC/'*-i V[/ G SU(6) (B.4) 

120^ 12^ 
+ ^Tr;7(Tr[/2)2Tr[/"-^ + ^{TiUfTiU^TrU''-^ - ^Tri/^Trf/^Trt/""^ 

--Tr?7Tr?7''TVC/"-^ + -Tr?7^TrC/"-^ - —(TrU)'^TrU''~^ + -(TrUfTrU^TrU''-'^ 

4 5 24^ ' 4^ ^ 

-^(Tr[/2)2Tr?7"-^ - ^TrUTrU^TiU''-'^ + jTtU'^TtU''-^ + ^{TrUfTrU''~^ 

^ o ^ ^ 

+Tri7Tr?7"~^ V?7 G SU(7) (B.5) 
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Mandelstam Constraints 



720^ ' 48^ ' 

-^(Tri7)^(Tr?7^)^Tr?7"~^ + ^(Tr?7^)^Tr?7"""^ - ^(Tr?7)^Tri7^Tr?7"~^ 

- — (TrC/)3TrC/2Trf/"-^ + ■^Trt/(Trf/2)2TrC/"-^ + -(Tr?7)2Tr?7-^Tr[/''-^ 
12 8 6 

--Tri/^Tri/^Tr?/"-^ - -Tri7Tri7''TrC/"-5 + -Tri/^Tr?/''-^ - — (Tri7)^Tri7"-'^ 
6 4 5 24 

4^ ^ 8^ ^ 3 4 

+ ^(Trf/)3TrC/''-^ - ;^Trf/TrC/2TrC/'*-3 + ^Trf/^TrC/'^-^ - ;^(TrC/)2Tri7"-2 

O ^ O ^ 

+^Tr?72Tri7""2 + Tri7Tr?7""^ V?7 G SU(8) (B.6) 



1.4 



1.3 - 



1.2 - 



1.1 - 



1.0 - 



0.9 - 



0.8 - 



0.7 - 



0.6 



A A A A ^ ^ 

□ □□□□□□□□□□□□□□□□□ 

□ 



j^XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXx 



+ + + + 





+ 


su 


[3) 


X 

+ 


X 


su 


4) 


X 




su 


5J 




□ 


su 


7) 


_x 

X 


■ 


su 


10 


+ + 


o 


su 


15 




• 


su 


20 




A 


su 


25 



I I I 

0.00 0.50 1.00 1.50 2.00 

1 

iVg2 





2.90 


1 1 1 1 1 




2.85 
2.80 


- ^^+(iV2) - 




2.75 


N 

N 


^ M 

<^ 


2.70 
2.65 


S 

S 

\ 

\ 

S 

\ 

N 

~~ \ ~ 

\ 




2.60 


\ 

\ 

\ 

~" \ ~ 

\ 

\ 




2.55 


\ 

\ 

\ 

\ 




2.50 


\ 

^ - 

\ 

\ 




2.45 


1 1 1 1 1 



0.00 0.02 0.04 0.06 0.08 0.10 0.12 



1 



3.0 



2.5 - 



SU(15) 
SU 16 
SU 18 
SU 20 
SU 25 



? 2.0 - 



X + ■ 

X X ^ + + 



.■■ * ^ + * 



1.5 ^wpf^ 



_■ « + 



1.0 



0.00 



0.02 



0.04 



0.06 



1 

iVg2 



0.08 



0.10 



Figure 1 




Figure 1 




